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ABSTRACT 

Sound  production  and  modification  in  an  Argon  discharge  plasma 
was  investigated  theoretically  and  experimentally.   Sound  pro- 
duction by  ions  in  the  cathode  region  was  shown  to  be  feasible  and 
consistent  with  earlier  experimental  results.   A  numerical  analysis 
of  the  theoretical  dispersion  relation  of  Ingard  and  Schulz  showed 
that  plasma  amplification  effects  will  not  compensate  for  acoustic 
losses.   The  theory  of  the  mobility-limited  thermionic  diode 
operated  as  a  microphone  was  expanded.   Such  microphones  were  used 
in  an  experimental  investigation  of  sound  production  and  modifi- 
cation by  a  discharge  plasma.   Possible  plasma  effects  on  sound 
propagation  were  observed.   The  characteristics  of  sound  pro- 
duction by  the  plasma  were  found  to  be  in  reasonable  agreement 
with  theory.   A  departure  from  theory  at  low  frequencies  was  inter- 
preted as  arising  from  the  effects  of  particle  diffusion,  an  effect 
not  previously  observed  in  acoustics. 
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I.   INTRODUCTION 

Plasma  acoustics  deals  with  the  interactions  between  charged 
particles  and  sound  fields.   The  problems  encountered  are  different 
from  those  of  ordinary  acoustics  since  charged  particles  experience 
bulk  forces  from  external  fields  whereas  neutral  particles  do  not. 
When  neutral  particles  are  present,  they  are  indirectly  influenced 
by  the  external  fields  through  their  interactions  with  the  charged 
particles.   The  complete  mechanism  of  interaction  is  very  complex, 
involving  ionization  and  recombination  as  well  as  elastic  and  non- 
ionizing inelastic  collisions. 

This  paper  focuses  on  the  problem  by  studying  the  effect  of 
the  interactions  on  the  neutral  particles.   Two  aspects  are 
treated.   One  is  the  creation  of  sound  by  the  current  modulation 
of  a  discharge.   The  other  is  the  modification  of  the  sound  when 
it  propagates  through  a  discharge.   The  two  aspects  are  closely 
r  elated. 

The  creation  of  sound  in  a  gas  requires  a  localized  time-de- 
pendent pressure  variation.   In  a  weakly  ionized  gas,  the  partial 
pressures  of  the  ion  and  electron  fluid  components  are  negligibly 
small  with  respect  to  the  partial  pressure  of  the  neutral  gas.   For 
this  reason,  the  acoustic  wave  made  in  a  discharge  plasma  is  con- 
cerned only  with  the  neutral  component  except  when  the  effect  of 


"Sound"  and  "acoustic  signal"  or  "acoustic  wave"  are  used 
synonomously  and  are  defined  to  mean  any  fluctuation  in  the 
pressure  of  the  neutral  particles. 


ions  and  electrons  on  the  neutral  component  is  of  interest.   The 
pressure,  expressed  as  p  =  kNT,  can  be  changed  by  either  a  change 
in  temperature,  T,  or  particle  density,  N. 

In  the  positive  column  of  a  glow  discharge,  energy  is  con- 
tinually fed  into  the  plasma  from  the  electric  field,  which  draws 
its  energy  from  an  external  power  supply.   The  energy  is  received 
by  the  electrons  and  ions  in  the  plasma  as  kinetic  energy.   Since, 
in  the  aggregate,  these  particles  are  traveling  at  their  terminal 
velocity,  referred  to  as  drift  velocity,  the  increased  kinetic 
energy  of  the  individual  electrons  and  ions  is  dissapated  through 
collisions  with  the  neutral  particles.   The  average  rate  of  energy 
drawn  from  the  field,  per  unit  volume,  can  be  written 

H  =  H  +  H.  =  eN  u  e'^  +  eN-M-.E"" 
e    1      e  e        x  i 

where  \^      and  M- .  are  the  electron  and  ion  mobilities.   Since  the 
e       1 

electron  mobility  is  much  larger  than  the  ion  mobility  and  the 
other  terms  are  equal,  the  ion  term  *is  negligible.   Thus  only  the 
electrons  need  to  be  considered  in  the  energy  transfer  process. 

In  the  equilibrium  situation,  the  energy  received  by  the 
electrons  must  be  balanced  by  losses.   The  energy  distribution  of 
the  electrons  closely  approximates  a  Maxwell  Boltzman  distribution, 
The  difference  is  due  to  the  effect  of  the  electric  field  and  the 
energy-dependent  cross  sections  of  the  neutral  particles  with 
which  the  electrons  interact.   The  electrons  transfer  energy  to 
the  neutrals  through  elastic  and  inelastic  collisions.   The  energy 
transferred  in  inelastic  collisions  is  usually  transferred  out  of 
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the  system  of  neutral  particles  through  radiation.   The  energy 
received  by  the  neutral  particles  in  elastic  collisions  is  thermal 
energy  which,  at  equilibrium  is  balanced  by  thermal  losses  to  the 
container  and  neutral  gas  outside  of  the  ionized  region.   This 
energy  transfer  between  the  hot  electrons  and  cold  neutrals  is 
comparable  to  mixing  two  gases  at  different  temperatures  except 
that  the  hot  gaseous  component  stays  hot.   The  elastic  collision 
losses  are. normally  much  larger  than  the  inelastic  losses  which 
can  be  neglected  for  tliis  analysis. 

Thus  the  processes  discussed  suggest  a  mixture  of  two  fluids 
of  different  temperature  which  retain  their  identity  in  the  mixture, 
The  hot  fluid  receives  energy  from  an  external  source  to  balance 
the  energy  transferred  within  the  mixture  to  the  colder  fluid.   The 
cold  fluid  is  hotter  than  its  surroundings,  to  which  it  loses 
energy . 

The  characteristic  times  for  the  three  different  processes; 
energy  gain  by  the  electrons,  transfer  to  the  neutrals  and  dis- 
sapation  to  the  surroundings,  are  different.   The  latter  process 
is  very  slow.   The  characteristic  time  for  energy  transfer  to  the 

neutrals  is  of  the  order  of  the  electron-neutral  collision  time, 

-8 
typically  ^  10    seconds.   Thus,  if  the  energy  transfer  to  the 

electrons  is  caused  to  vary  at  acoustic  frequencies,  by  varying 
N   or  E,  the  time  delay  can  be  neglected  and  the  energy  transfer 
to  the  neutrals  can  be  considered  in  phase  with  the  energy  transfer 
to  the  electrons.   Since  the  time  for  energy  transfer  to  the  sur- 
roundings is  much  greater,  the  temperature  change,  and  therefore 
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the  pressure,  of  the  neutrals  will  be  essentially  in  phase  with 
the  energy  transferred  in.   This  means  that  a  change  in  either  E, 
the  electric  field  or  N  ,  the  electron  density,  can  generate  a 
time  dependent  pressure  component  in  the  neutral  gas  which  can  act 
as  an  acoustic  source. 

Amplification  of  a  wave  occurs  when  the  wave  causes  a  change 
in  the  medium  which  leads  to  a  change  in  the  wave  in  phase  with 
itself,  i.e.  "positive  feedback"  occurs.   This  can  occur  in  a 
discharge  plasma  through  the  interaction  of  the  pressure  wave  in 
the  neutral  gas  with  the  electron  energy  source  term  H  .   The 
neutral  gas  density  changes  in  phase  with  the  acoustic  pressure. 
Since  both  the  mobility  and  electron  density  depend  on  the  neutral 
density.  +be  f feedback  can  be  introduced  through  thtrm  even  v;ith  no 
change  in  the  electric  field. 

In  the  simplest  analysis,  the  electron  density  varies  in  phase 

with  the  neutral  density.   If  n   and  n   are  the  perturbations  in 

en  ^ 

the   electron   and   neutral    gas    densities    N      and  N    ,    then 

en' 

n 

n     77-  N  . 
e    N    e 
n 

This  indicates  that  the  energy  transfer  to  the  neutral  gas  varies 
in  phase  with  the  acoustic  pressure  variation.   Since  this  would 
produce  a  pressure  variation  in  phase  with  the  original  pressure 
variation,  it  constitutes  an  amplification  mechanism.   For  net 
growth  to  occur,  the  acoustic  losses  due  to  thermal  conductivity 
and  viscosity  must  be  less  than  the  acoustic  source  gain. 


12 


A  more  sophisticated  analysis  indicates  that  phase  differences 
occur  between  the  electron,  ion  and  neutral  density  variations  with 
resultant  generation  of  a  space  charge  electric  field.   The  phases 
of  all  of  these  quantities  influence  the  phase  of  the  feedback  and 
determine  whether  the  plasma  effects  increase  or  decrease  the  net 
attenuation  of  the  acoustic  wave. 

A.   CHRONOLOGY  OF  THE  EXPERIMENT 

Interest  in  plasma  acoustics  here  at  the  Naval  Postgraduate 
School  arose  out  of  the  study  of  moving  striations.   Earlier  work 
performed  here  was  that  of  Partlow  [l],  Carreta  and  Moore  [2], 
Wilson  [3],  Cranda^ll  [4]  and  Melville  [5].   Melville's  work 
represented  a  shift  away  from  the  study  of  moving  striations.   This 
was  prompted  by  a  theoretical  paper  by  Ingard  in  I966  [6],  which 
predicted  amplification  of  an  acoustic  wave  in  a  discharge  plasma. 
A  more  comprehensive  theory  by  Ingard  and  Schulz  followed  in 
1967  [7]. 

The  experimental  investigation  described  in  tliis  work  began 
in  1968  as  a  natural  continuation  of  the  earlier  work  performed 
here.   The  earlier  work  here  and  elsewhere  had  demonstrated  that 
sound  could  be  generated  by  a  plasma.   This  coupled  with  the 
prediction  of  effects  on  the  sound  as  it  passed  through  a  discharge 
led  the  author  to  the  consideration  of  a  combined  experiment.   In 
the  proposed  experiment,  a  discharge  would  be  used  to  create  sound, 
and  it  would  be  observed  before  and  after  it  propagated  through  the 
discharge.   With  the  proper  techniques,  the  experiment  could  provide 
information  about  both  the  creation  of  sound  and  its  modification 

in  the  discharge. 
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A   numerical    analysis    of   the    Ingar d-Schulz    dispersion   relation 
contained    in    [7]    was    performed  as    part    of    this    investigation    in 
which   the    variables    were   neutral    density,    neutral    temperature, 
electron    density,    electron    temperature  and   angular    frequency    of 
the   acoustic    signal.      The   attenuation  was    minimized  with    no  con- 
straints   placed   on    the    values    of    the   variables    except    for    upper 
and   lower    limits.      The    results    showed   that    the  minimum   attenuation 
occurred  at    high    neutr£il    densities,    high   electron    densities    and 

temperatures,    low   neutral    temperatures    and    low   frequencies.      The 

.  .    .   2 

maximum  relative  ampliiication   was  also  investigated.   The  con- 
clusion drawn  from  a  study  of  the  behavior  of  these  two  quantties 
was  that  the  experiment  should  test  the  theory  in  the  medium 

observed  by  Yatsui,  et.al.  [8,9]  indicated  the  need  to  design  the 
experiment  so  that  the  polarity  of  the  discharge  could  be  reversed. 
Also,  it  was  felt  that  thermocouples  were  needed  to  determine  the 
neutral  gas  temperature  in  and  out  of  the  discharge. 


2 
"Relative  amplif ic:at ion"( Ap)  is  the  change  in  attenuation 

due  to  the  plasma  amplification  referred  to  the  classical  atten- 
uation.  It  is  defined  as  ♦ 

=  A 


where  p  is  the  attenuation  coefficient,  I  Signifies  ionized  and  u 
signifies  unionized.   A  negative  value  means  that  an  amplification 
is  operative  in  the  plJisma  case. 
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The  experiment  was  designed  to  use  two  mobility-limited 
thermionic  diodes  as  microphones  to  observe  the  signal  at  two 
points  with  the  signal  being  produced  by  a  remote  source  discharge. 
Such  an  arrangement  would  permit  a  direct  determination  of  the 
attenuation.   The  theory  of  the  diode  microphones  was  developed 
beyond  the  earlier  elementary . form  so  that  an  extensive  cali- 
bration was  not  needed  to  use  it  as  an  analytical  tool.   Technical 
problems  were  solved  so  that  the  tube  could  be  constructed  with 
a  smooth  surface,  free  from  side  arms,  bumps  and  bubbles  which 
affect  an  acoustic  signal,  in  spite  of  the  numerous  tube  elements 
contained.   A  circuit  was  designed  to  control  the  discharge  cur- 
rent to  produce  the  sound  in  the  desired  form. 

"ulscd  currents  v,'crc  uccd  to  produce  the  sound  v/ith  the  sonrc*^ 
discharge.   The  current  and  sound  waveforms  were  Courier  trans- 
formed to  determine  their  spectra.   By  a  process  of  comparison, 
extrapolation  and  comparison,  the  frequency  dependence  of  the 
acoustic  attenuation  and  the  source  transfer  function  were  obtained, 
The  measured  values  of  the  experimental  parameters  were  used  to 
determine  the  predicted  values  of  attenuation  from  the  Ingard- 
Schulz  equations.   These  equations  predicted  that  the  plasma  would 
not  affect  the  attenuation  under  the  conditions  used.   The  experi- 
mental value  for  the  attenuation  co-efficient  agreed  with  the  pre- 
dicted value  within  ex^jer  imental  accuracy  although  the  results 
indicated  the  possibility  of  a  plasma  effect. 
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The  results  of  Yatsui  et.al.  [8,9]  were  then  reevaluated. 
It  was  noted  that  while  the  relative  change  in  the  attenuation 
coefficient  was  small,  since  the  attenuation  was  very  large,  the 
amplification  effect  was  large.   In  view  of  these  results,  it 
was  recognized  that  the  quantity  to  optimize  is  the  difference 
in  attenuation  co-efficients  between  the  ionized  and  un-ionized 
cases . 

Earlier  a  theoretical  study  of  the  frequency  characteristics 
of  a  discharge  sound  source  was  undertaken  by  Schulz  [ lo] .   He 
concluded  that  the  energy  transfer  from  the  electrons  was  the 
dominant  source  mechanism.   The  predicted  frequency  dependence 
is  sin  (V n  a/c)/(V  TT  a/c) .   When  the  source  transfer  function  was 
pypf^r  imc^nta  1  1\;  determined,  it  was  found  to  be  consistent  with 
theory  at  high  frequencies  but  not  at  low  frequencies.   The  low 
frequency  effect  was  attributed  to  diffusion  effects. 

In  the  following  two  sections  of  this  chapter,  earlier  experi- 
ments relating  to  the  creation  and  m'odif icat ion  of  sound  are 
discussed.   Following,  Chapter  III  discusses  the  theoretical  treat- 
ment of  sound  creation  and  modification.   This  order  of  presentation 
was  chosen  because  it  was  felt  to  offer  a  better  insight  into  the 
aspects  which  the  theory  should  treat. 

B.   PREVIOUS  EXPERIMENTAL  WORK  ON  ACOUSTIC  INTERACTIONS 

Before  considering  specific  papers  of  direct  interest  to  this 
experiment,  a  quick  general  review  of  previous  experimental  investi- 
gations of  the  acoustic  interaction  is  in  order.   Interest  in  the 
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interaction  between  charged  particles  and  neutral  acoustic  waves 
is  recent.   A  relationship  between  moving  striations  and  acoustic 
waves  was  suggested  by  Watanabe  and  Oleson  [ll]  in  1955. 
Wojaczeck  [l2]  reported  observations  of  a  sound  interaction  in 
i960  during  experiments  in  Xenon.   He  noted  on  effect  in  his  main 
discharge  which  occurred  when  an  auxiliary  discharge  was  pulsed, 
and  he  identified  the  effect  as  due  to  sound  waves  generated  by 
the  pulsed  discharge. 

Since  that  time,  many  workers  have  reported  observations  of 
sound  waves  generated  by  discharges.   The  source  and  type  of 
disturbance  in  the  discharges,  which  were  observed  to  produce 
sound,  have  been  sufficiently  varied  to  conclude  that  it  is  a 
general  phenomenon.   Sound  has  been  observed  l1-3J  whose  source 
was  interpreted  to  be  the  current  fluctuations  associated  with 
moving  striations.   Sound  production  by  artificially  introduced 
current  fluctuations  has  been  reported  [4,13-15]  for  a  variety  of 
conditions.   Modulation  of  the  electron  temperature  in  a  discharge 
by  an  R.F.  field  was  observed  by  Fitaire  and  Mantei  [16]  to  produce 
sound.   Multiple  pulse  emission  from  an  HCN  laser,  reported  by 
Turner,  et.el.,  [l7]j  was  interpreted  by  Born  [18]  as  due  to  sound 
production  by  the  discharge. 

The  variety  of  techniques  used  to  study  the  sound  produced  by 
pulsed  discharges  is  almost  as  interesting  as  the  observations 
reported  [4,5,12,19-25].   The  pulses  ranged  in  duration  from  one  to 
fifty  microseconds  [23,18],  in  pressures  ranging  from  5  x  10    torr[l9J 
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to  20  torr  [l8j.   In  some  experiments,  the  sound  was  detected  by 
observing  its  effect  on  some  property  of  the  discharge  which 
produced  the  sound.   An  example  of  this  is  the  work  reported  by 
Berlande,  et.al.  [24],  in  which  they  observed  perturbations  in 
the  light  of  the  afterglow  of  the  plasma  which  produced  the 
sound.   Nygaard  [2l]  utilized,  the  light  perturbation  effect  by 
using  one  discharge  to  produce  the  sound  while  he  observed  the 
afterglow  of  a  second  discharge.   Goldstein  et.al.  [23]  used  a 
mobility-limited  thermionic  diode  as  a  microphone.   These  and 
other  techniques  were  used  in  some  cases  because  of  the  problems 
associated  with  conventional  acoustic  detection  techniques  under 
the  conditions  involved.   In  other  cases,  they  wece  used  because 
•♦-'Vi/->  sc^'r^'^  a*^'^  z."'"s  *^^^'t''^o+T_Q»^  ivas  a*^  n*^'^ ^i-r^ r^^^+Qr?  r'^su"^'*"  — ^  a*^ 
exper  iment . 

Only  a  few  papers  have  been  published  concerning  the  effect 
of  a  plasma  on  sound  propagating  through  it  [5>8,9,26],  as 
opposed  to  the  effect  on  the  discharge.   Because  of  their  relevance 
and  sparsity,  these  papers  will  be  discussed  in  detail  in  the  next 
section. 

C.   SPECIFIC  EXPERIMENTS  OF  INTEREST 

The  two  aspects  of  this  experiment  are  the  creation  of  sound 
and  its  modification  in  a  gas  discharge.   Earlier  experiments  have 
examined  both  of  these  and  as  the  experimental  results  have  been 
reported,  theories  have  been  developed  to  explain  the  experimental 
results.   The  theoretical  aspects  of  the  experiment  will  be  dis- 
cussed in  the  next  chapter. 
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1 .       Sound  Production 

A  review    of   the    existing   experimental    evidence   concerning 
sound  production    can   point    out    features    which   must   be    explained 
by   any   theoretical   model    of   the  process. 

Wilson   [3]    determined   from  phase   measurements    of   the   sound 
produced   by    moving    striations,    that    the   only    sound   source   position 
consistent    through    the    experiment   was   at    or    near    the   cathode.       By 
comparing    the    current    and   sound   waveforms,    he   concluded   that    the 
sound   was    produced   by    fluctuations    in   heating   processes    occurring 
near    the   cathode    due    to   current    fluctuations.      He   suggested   that 
the  mechanism   must    involve   an    interaction    between    the    fluctuations 
in    ion    density   and   temperature   and   the   neutrals. 


iNygaar  u   anu   Mej-A.-^    l^-lj    are   niojie    SijeCxXxc    da     cO    Ciitr    iOc««,ivjn 


*.  —  -..1, 


of  the  source.   They  state  that  the  acoustic  disturbance  from  their 
H.V.  pulse  discharge  originates  in  the  negative  glow  region.   They 
postulate  a  high  pressure  driver  region  near  the  cathode  maintained 
by  electron  heating  and  momentum  transfer  from  the  ions. 

Hayess  [20]  also  interpreted  the  sound  production  as  arising 
from  localized  heating  of  the  neutral  gas.   However  his  experiment 
brings  in  a  different  aspect.   He  used  a  high  voltage  pulse  applied 
to  an  electrode  near  a  discharge  already  running,  so  that  the  other 
electrode  of  the  circuit  had  a  much  greater  charge  carrier  con- 
centration than  the  exciting  electrode.   He  found  that  the  sound 
source  was  always  at  the  exciting  electrode,  regardless  of  the 
polarity  of  the  voltage  pulse.   By  varying  the  rise  time  of  the 
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voltage  pulse,  Hayess  established  that  it  had  particular  importance 

in  the  sound  excitation  process.   Pulses  with  rise  times  of 

5  X  10   V/sec  produced  sound,  those  with  rise  times  of  1.4  x  10  V/sec 

did  not.   In  addition,  sine  waves  with  an  average  rise  time  of 

7 
2  X  10   v/sec  did  not  lead  to  sound  production.   Hayess  also 

established  that  the  presence  of  charge  carriers  at  the  exciting 

electrode  reduced  the  sound  produced.   From  consideration  of  these 

factors,  Hayess  concluded  that  energy  gained  in  the  high  field  by 

the  charge  carriers  is  transferred  to  the  neutrals.   He  further 

concluded  that  the  ions  were  responsible  for  the  inert ial  phenomena, 

i.e.  dependence  on  rise  time. 

In  contrast  to  the  above  experiments,  the  observation  of 
smjn-i  production  by  modulation  of  the  electron  temperature  in  a 
discharge,  reported  by  Fitaire  and  Mantei  Ll6],  may  be  deemed  as 
experimental  confirmation  of  the  roll  of  the  electron  temperature 
in  the  acoustic  source  term  of  Ingard  and  Schulz's  equations  L?!. 
In  interpreting  these  experimental  results,  it  must  be  borne  in 
mind  that  there  is  a  distinction  between  the  energy  of  electrons 
in  a  low  field  region,  such  as  the  positive  column  and  a  high 
field  region  such  as  the  cathode  fall.   In  the  latter,  the  motion 
is  highly  directional  and  can  not  be  considered  thermal. 

Upon  considering  the  above  experimental  results,  the 
following  conclusions  may  be  drawn.   There  are  two  mechanisms 
through  which  sound  is  produced  in  a  discharge.   One  mechanism 
occurs  near  the  cathode  and  probably  involves  ions  in  a  high  field, 
the  cathode  fall.   Nygard's  contention  that  the  source  is  the 
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negative  glow  region  is  in  contradiction  with  the  results  of 
Hayess  concerning  the  high  field.   The  other  mechanism  occurs  in 
the  positive  column  and  involves  the  electron  temperature.   The 
former  mechanism  appears  to  dominate  the  production  process. 
Further  comment  shall  be  reserved  for  a  discussion  of  sound  pro- 
duction in  Chapter  II. 

2  .   Sound  Modification 

The  experimental  evidence  concerning  the  effect  of  the 
plasma  on  sound  is  not  so  conclusive.   Wojaczek  Ll2j  observed  the 
perturbation,  produced  by  the  sound,  in  the  light  emitted  from  a 
Xenon  gas  discharge.   On  the  basis  of  his  calculations,  he  con- 
cluded that  the  attenuation  was  much  greater  than  expected.   He 
attributed  this  to  possible  in  terac  t  ion  Vvith  the  discharge. 

Later  Hayess  [19]  pointed  out  that  the  attenuation 
Wojaczek  predicted  was  in  error  because  the  dominant  loss 
mechanism,  wall  effects,  was  neglected.   Hayess  recomputed  the 
predicted  values.   With  the  wall  losses  taken  into  consideration, 
the  experimental  values  agreed  with  the  predicted  values  within 
the  experimental  uncertainty.   Wojaczek 's  experiment  was  conducted 
at  pressure  from  .07  torr  to  1.5  torr  with  currents  from  1.5  to 
0  amperes. 

Working  in  Argon,  Yatsui  et.al.  [8,9]  studied  the  effects 
of  a  discharge  on  ultrasonic  sound.   They  used  barium  titanate 
transducers  with  resonant  frequencies  of  50  kHz,  100  kHz  and 
200  kHz  for  transraitte;r s  and  receivers.   The  signal  received  with 
no  discharge,  as  a  function  of  transducer  input  voltage,  was  the 
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reference.   By  dividing  this  into  the  results  obtained  with  the 
discharge  on, a  "spatial  growth  rate"  was  obtained.   A  definite 
change  in  the  growth  rate  was  observed.   It  was  found  to  be  de- 
pendent on  the  current  and  pressure  of  the  discharge  and  on  the 
polarity  with  respect  to  the  propagation  direction.   The  polarity 
effect  was  so  strong  that  under  conditions  such  that  while  propagation 
in  the  direction  of  electron  drift  motion  enhanced  the  signal  by 
50%;  no  signal  v/as  received  when  propagated  contra  to  the  electron 
drift.   The  variation  of  the  growth  rate  with  pressure  and  current 
was  not  monotonic,  as  the  signal  exhibited  a  maximum  enhancement 
at  intermediate  values  of  the  range  investigated.   The  conclusion 
drawn  by  the  authors  was  that  strong  damping  occurred  for  contra- 
propanating  signals  binder  all  conditions.   For  parallel  prooagation. 
there  is  a  critical  pressure  p   above  which  attenuation  occurs  and 
below  which  amplification  occurs.   They  attributed  the  directional 
effect  to  coupling  between  the  sound  and  the  drift  motion  of  the 
electrons,  a  coupling  v;hich  is  strongest  when  the  sound  velocity 
is  nearly  equal  to  the  electron  drift  velocity. 

Their  experiment  provides  experimental  support  to  the  con- 
cept of  energy  transfer  from  the  electric  field  to  the  neutrals 
via  the  electrons.   It  should  be  pointed  out  that  the  acoustic 

frequency  is  of  the  order  of  the  collision  frequency  at  the  low  end 

-2 
of  their  pressure  range,  10    to  .5  torr .   Indeed  it  was  at  this 

end  that  the  maximum  effect  was  observed. 

Macon  [26]  used  a  Helium  discharge  to  investigate  the 

effect  of  the  discharge  on  sound.   He  employed  resonant  cavity 
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techniques  with  frequencies  below  2000  Hz.   The  sound  source  was 
calibrated  to  provide  constant  amplitude  deflections  over  a  range 
of  frequencies.   A  rigidly  mounted  microphone  received  the  signal 
at  the  other  end  of  the  tube.   By  comparing  the  ratio  of  pressures 
at  adjacent  pressure  maxima  and  minima,  and  assuming  the  net  loss 
was  dominated  by  boundary  effects,  the  attenuation  constant  was 
determined.   The  value  for  the  "discharge  off"  condition  was  com- 
pared with  the  "discharge  on"  condition.   Measurements  for  the 
discharge  off  condition  were  taken  one-second  after  turning  the 
discharge  off  in  order  to  have  comparable  temperature.   The  tem- 
peratures, as  determined  from  sound  velocity  measurements,  dropped 
from  450  K  with  the  discharge  on  to  36O  K  with  the  discharge  off, 
with  one-'^econd.   The  results  indicated  a  decrease  in  a  L  Leuua  t  ion 
rate  of  30%;  21  value  too  large  to  be  accounted  for  by  the  pressure 
and  temperature  differences.   Although  the  experimental  results  for 
attenuation  and  amplification  rates  were  both  in  poor  agreement 
with  theoretical  predictions,  Macon  concluded  that  amplification  by 
the  electrons  was  definitely  occurring. 

These  experiments  indicate  that  there  are  amplification  and 
attenuation  effects.   Wojaczek's  experiment  does  not  hint  of  any 
directional  effect.   Macon  did  not  look  for  any.   The  directional 
effect  observed  by  Yatsui  is  explainable  by  models  for  other  fields, 
namely  electronics  and  solid  state.   The  Ingar d-Schulz  model  to  be 
discussed  in  the  next  section  does  not  take  directional  effects 
into  consideration. 
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II.   THEORY 

Physical  theories  go  hand  in  hand  with  experiments.   Theories 
are  developed  to  account  for  experimental  results.   Experiments 
are  conducted  to  evaluate  theories  by  testing  their  predictions. 
In  the  previous  chapter,  the  experimental  evidence  regarding  pro- 
duction and  modification  of  sound  by  a  discharge  plasma  was  re- 
viewed.  In  this  chapter,  the  theories  will  be  examined. 

A.   SOUND  PRODUCTION  BY    DISCH/\RGE  PLASMAS 

After  several  investigators  reported  observation  of  sound 
produced  by  electrically  modulated  discharges,  Ingard  [6]  performed 
a  theoretical  analysis  of  the  problem.   He  derived  the  wave  equation 
for  the  pressure  in  the  neutral  gas  in  which  the  interactionb  with 
the  electrons  and  ions  were  expressed  in  terms  of  sources  in  the 
wave  equation.   Neglecting  the  fluctuations  in  neutral  particle 
density  resulting  from,  the  imbalance  of  ionization  and  recom- 
bination rates,  Ingard  obtained  the  folloiving  wave  equation: 

c   ot  c 

In  this  equation,  p  is  the  acoustic  pressure,  defined  as  the  time- 
dependent  pressure  variation  about  the  ambient  pressure;  c,  the 
speed  of  sound;  and  y  is  the  ratio  of  specific  heats.   "H"  is  the 
energy  source  term,  representing  the  rate  of  energy  transfer  to  the 
neutral  gas  from  the  charged  particles  per  unit  volume.   The  rate 
of  transfer  of  momentum  to  the  neutral  gas  per  unit  volume  is  re- 
presented  by  the  momentum  source  terra  F. 
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Ingard   considered   the    effect    of   momentum    transfer    represented 
by   V'F   as    negligible   for    two   reasons.       First    the    electron    and    ion 
contributions    tend   to   cancel    each    other.      Second,    only    the    spatial 
variation    of   F   leads    to   sound,    and   for    most    cases,    this    is    very 
small.      For    the    energy    source   term,    he   wrote 

4m      /ra    V        I 

H    =   S         ^^     \  N   N    <ct)  (2.2) 

m^   \^     2      y      e   n 

where  v   is  the  average  electron  velocity  and  (ct)  is  the  average 
e 

collision  cross  section  for  elastic  scattering.   Ingard  considered 
the  ion  contribution  negligible  due  to  its  much  lower  mean  energy. 
Ingard  pointed  out  that  since  the  electrons  were  in  an  equilibrium 
state,  the  energy  lost  to  the  neutrals,  represented  by  H,  was  the 
same  as  the  energy  gained  from  the  field.   Therefore 

H  =  N  eE(M.  E)  (2  .3) 

e   ^  e  '  ' 

where  a   is  the  electron  mobility  and  E  the  electric  field, 
e 

Schulz  [lO]  used  equation  (2.1)  'as  a  basis  to  develop  an 
expression  for  the  sound  field  produced  by  current-modulating  a 
discharge  operating  between  two  planar  electrodes.   He  considered 
the  effects  of  the  two  source  terms  separately  in  calculating  the 
pressure  response.   For  both,  a  one  dimensional  problem  was 
assumed  and  it  was  assumed  that  the  time  dependence  was  given  by 
exp(-icot).   Under  these  assumptions,  the  wave  equation  inside  the 
source  region  becomes 


2=^ 


If  H  is  assumed  to  be  spatially  constant  within  the  region, 

then  the  general  form  of  the  solution  inside  of  the  region  is 

^      «  -ikz   „  ikz    .  y-l  ,,  ,_  ^. 

P.   =  Ae     +  Be     +  i  ^-—-   H.  (2.5) 

m  CO 

With    the   additional    assumption    of   symmetry   about    the   origin, 
chosen    to   lie   between    the    two   electrodes    which   are   separated  by 
a    distance  a,    Schulz's    expression    for    the  pressure    field    outside 
the   source   region    for    the    energy   source   is    obtained. 

P^=l2^.„,i„^||)e>^(ih|f).  (2.6, 

Schulz  assumed  tha't  the  momentum  transfer  acoustic  source 
arose  from  an  unbalanced  fluctuation  between  the  forces  exerted 
by  the  ions  and  the  electrons  on  the  neutral  gas.   He  stated  that 

this  force  per  unit  volume  was  given  by 

— »  — ♦ 

F  =  N.W.   --^^^^  -  N   "^en^^   .  (2.7) 

1  m  CO.  -100     e  -r-  ' 

m         CO   -xOO 
en 

With   the   assumptions    CO.   >^C0  and   CO     »C0.     ,    this    simplifies    to 
^  m  en         in  ^ 

F  ^    iN.  qE (CO/CO.     )     .  (2.8) 

--♦ 
If  the  electric  field  \l    is  assumed  to  be  constant  between  the 

two  electrodes  and  zero  outside,  then  V'F  can  be  expressed  as 


V-F  -  F  [5(z+|)  -  5(2-|)] 


where   5(z)    is    the    Dirac    delta    function.      This    leads    to   a   pressure 
field   of: 


cqN. 


(f)exp(.   icof^   .  (2.9) 


Schulz  concluded  that  the  momentum  source  term  would  be  negligible. 

He  compared  the  two  sound  fields  and  found  that  the  ratio  was  of 

CO, 


^in    /^e  \'^ 

— —    I       in    favor    of    th< 


the   order    of  (  ^T"    /       in    favor    of    the    energy    transfer    process. 

Schultz's    results   predict    that    the   energy    transfer    mechanism 
dominates.      This    result    is    in   agreement   with    Ingard's    earlier    [6] 
conclusion.      Schulz 's    expression    for    the   sound    field  predicts 

c  "I  ■p*  \r  r\TT 

that  the  frequency  dependence  is  like  where  x  =  V  —  . 

Heretofore  this  frequency  dependence  has  not  been  experimentally 
investigated.   However  there  is  experimental  evidence  which  refutes 
Schulz 's  results  and  contradicts  Ingard's  assumptions  concerning 
the  energy  source  term. 

As  mentioned  earlier^  the  experimental  evidence  suggests  that 
the  principal  source  term  lies  near  the  cathode  and  involves  the 
field.   Schulz 's  analysis  implied  that  the  positive  column  was 

the  source  region.   The frequency  dependence  also  shows 

that  the  sound  produced  by  high  frequencies  will  be  small.   The 
effective  frequencies  of  the  sound  pulses  mentioned  earlier  have 
been  in  the  range  of  20  kHz  to  1  MHz,  all  of  which  are  high 
acoustic  frequencies.   Perhaps  the  comparison  with  pulsed  dis- 
charges is  not  a  valid  one  because  of  the  difference  between  the 
experiments  and  the  conditions  assumed  for  the  derivation.   The 
derivation  presumes  small  perturbations  yet  the  current  changes 
in  a  pulsed  breakdown  discharge  can  not  be  small  perturbations. 

Since  sound  production  appears  to  occur  in  the  vicinity  of  the 
cathode  and  appears  to  be  associated  with  high  fields,  a  possible 
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source   terra    compatible   with    experimental    results    would   be   a    field- 
dependent    energy    source   term.      A    field-dependent    source    is    one    in 
which   the    energy    for    the   acoustic    wave   comes    from   the    field.      The 
dependence   on    this    term    enters    through    the   energy    equation    in 
contrast    to   Schulz's    treatment,    in   which    the    field   entered   through 
the   momentum    equation. 

When    the   processes    occurring    in    the   region    of    the   cathode  are 
considered,    a    mechanism   for    the   production    of   sound   suggests 
itself    to   the   author.      The    contribution    of    ions    to   the    energy 
transfer    is    usually    neglected.       Due    to    their    higher    mobility    the 
electrons    transfer    energy    from    the   field   to   the   neutral   particles 
much   more  readily.       It    should   be   recognized   that    due    to    their 
cominorj   macsos  ,    the    icns    and    neutrals    will   bo   ?.t    es*-ept  ^'^  ^  ^  ^'    I'Vio 
same    temperature.      The   unity    mass    ratio   means    that    during   a 
collision,    an    ion    will    transfer    a    large   fraction    of    its    energy. 
In    the    cathode   region,    this    factor    becomes    important.      There    the 
ions    are    the   majority    carriers    of   the    discharge   current.      Their 
density    is    higher    and   they    accelerate    toward   the    cathode   under    the 
high   field   of   the    cathode    fall.      The    increased   energy   will    be 
transferred   to   the   neutrals    through    collisions.      There   will    also 
be  a    net    momentum    transfer    since   the    ions    acquire   momentum   from 
the    field  and   transfer    it    also. 

A  rough   approximation   for    the   strength    of   such   a    source   term 

can   easily    be   developed.      The   average   energy    gained   from    the    field 

between    collisions    is    eE^ .      where   E   is    the    electric    field  strength 

in 

and  i.       is  the  mean  free  path  of  the  ions.   Since  approximately 

half  of  the  ions  energj'  would  be  transferred  in  a  collision,  the 
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energy    transfer    rate  per    unit    volume    is  • 

H.    "^  -^  eEi  .      CO.      N.  .  '  (2.10) 

1         2  in       m       1 

This    approximation    is    low   because    it    uses    the   average   energy    that 
would   be    transferred    if   the    ions    were   stopped  after    a    collision. 
It    does    provide  an    order    of  magnitude   estimate.       If  -^  .       is 
approximated  by    v/CO.     ,    then 

H.^^eE^.N..  (2.11) 

1         2  1    1  ^  ' 

The  ratio  of  this  term,  where  E  =  E  ,  the  field  in  the  cathode 

c 

fall  to  the  electron  source  term,  E  =  E  ,  the  field  in  the  positive 
column  is 


H .    '2eE   V.N. 
-2.^   £ i-i. 


H     N  eE  (M-  E  ) 
e    e  p  ^  e  p' 


(2.12) 


Now   eN  |i.  E   is  the  current  density  in  the  positive  column,  and 
e  e  p 

e  V.N.  is  approximately  the  current  density  in  the  cathode  region, 
11 

Since  the  two  current  densities  should  be  the  same, 

H.  E 
e  p 
The  approximations  used  in  arriving  at  this  result  are  equivalent 
to  saying  that  in  the  positive  column,  all  of  the  power  delivered 
goes  into  heating  the  neutral  gas  and  in  the  cathode  region,  part 
goes  into  heating  the  neutral  gas.  The  ions  accelerate  under  the 
action  of  the  field  in  the  cathode  fall  and  gain  kinetic  energy. 
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The  difference  between  the  energy  taken  from  the  field  and  the 
kinetic  energy  gained  represents  the  energy  transferred  to  the 
neutrals  through  collisions.   Since  the  current  density  is  con- 
stant but  the  field  varies  throughout  the  length  of  the  discharge, 
the  energy  transfer  rate  is  a  function  of  the  field  and  the 
fraction  of  energy  lost  in  collisions. 

B.   SOUND  MODIFICATION  BY  THE  DISCHARGE 

The  theoretical  problem  of  acoustic  waves  in  a  weakly  ionized 
gas  has  been  studied  by  several  authors  other  than  Ingard  l27-33J« 
Without  exception,  the  work  of  these  authors  needs  modification 
before  it  can  be  used  to  predict  the  results  of  an  experimental 
study  in  a  laboratory.   This  is  because  they  neglect  boundary  loss 
effects;  yet  at  the  low  pressures  of  most  discharge  experiments 
with  typical  tube  diameters  of  less  than,  or  of  the  order  of,  5  cm, 
boundary  losses  are  dominant.   The  work  of  Ingard  and  Schulz  L7J 
is  the  most  applicable  to  an  experiment.   In  addition,  it  is  a 
general  treatment  and  is  reasonably  susceptible  to  numerical 
evaluation . 

The  mechanism  of  amplification  is  embodied  in  the  concept  of 
energy  transfer  from  the  electrons  to  the  neutrals  through  col- 
lisions.  In  the  absence  of  an  acoustic  wave,  the  energy  transferred 
is  spatially  uniform  and  merely  serves  to  increase  the  temperature 
of  the  neutral  gas.   In  the  presence  of  an  acoustic  wave,  the  energy 
transfer  becomes  time  dependent  and  spatially  non-uniform.   This 
produces  an  acoustic  v/ave  which  can  be  in  phase  with  the  initial 
wave,  leading  to  amplification. 

30 


Ingard's  initial  analysis  of  the  ptoblem  [6]  followed  the 
experimental  observations  of  sound  produced  by  electrically 
modulated  discharges  L2 , 12 , 13  ,24] .   He  treated  the  problems  as 
belonging  to  two  categories,  sound  produced  by  current  fluctuations, 
or  an  acoustic  ivave  instability.   For  the  first  category,  he  derived 
equation  (2.1).   To  treat  the  second  category,  he  placed  a  terra 
representing  the  electron  density  perturbation  into  the  energy 
source  term  of  equation  (2.1).   His  results  indicated  that  the 
acoustic  signal  from  a  current  fluctuation  and  also  the  growth 
of  an  acoustic  fluctuation  would  be  small  at  room  temperatures  but 
enhanced  at  cryogenic  temperatures. 

Although  published  in  May  I966,  the  development  of  reference  6 
was  done  in  early  1965.   In  xhe  fall  of  19"6 ,  Macon  l26j  compleLeu 
his  investigation.   A  thesis  student  of  Ingard's,  he  attempted  to 
conduct  his  experiment  at  cryogenic  temperatures  but  had  to  use 
room  temperatures  because  of  noise  due  to  moving  striations.   He 
concluded  that  amplification  was  taking  place  at  temperatures 
greatly  above  room  temperature  (450  K).   His  values  of  growth  and 
attenuation  rates  iver e  not  in  agreement  with  theory  but  he  attri- 
buted this  to  his  experimental  arrangement.   The  more  comprehensive 
treatment  of  the  problem  by  Ingard  and  Schulz  [7]  was  published  in 
1967. 

The  Ingar d- Schulz  analysis  is  based  on  a  three  fluid  model  of 
the  ionized  gas  with  both  energy  and  momentum  transferred  between 
the  fluid  components.   The  effects  of  viscosity  and  heat  conduction 
within  each  component  and  v/ith  the  walls  are  also  taken  into 
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consideration.      The    conditions    assumed    in   making   approximations 

may    be   considered   typical    for    a    glow   discharge:    gas   pressures    of 

10  12         -3 

1-10    torr,    electron    densities    of    10      -    10         cm      ,    and   acoustic 

frequencies    10  ^   CO  ^    10^. 

The   equations    forming    the  basis    for    the   analysis    are    introduced 
as    linearized   space-time   Fourier    transforms    of    the    continuity 
equation    (2.l4),    the  momentum    transfer    equation    (2.15),    the    energy 
balance   equation    ( 2 . l6 )    and   the   equation    of   state    (2.17)    for    each 
fluid   component    (these   equations   are   presented   here   for    complete- 
ness   and  all    items    will    not    be   defined  at    this    time). 


-    icon  .    +    ikN  .v  .    =   0 
D  3    3 


(2.14) 


-    icov  .+    i-: —  D.    =   CO..  (v,  -v.)+   cu.^(v^-v.)+    (q./m.)e-/    .(k,a:)v.    (2.15) 
3        P  j   "  J  -^'^      ^      -^  J'^      ^       J  JJ  u  J 

(j,k,>C    -    n,i,e  and  permutations) 

■iCO   c     .Q  .    =  a.(n    /N    )+    B.(0    /T    )-    ik(KT./m.)v. 
v3    3  J      e      e'      ^j'    e      e^  3      3      3 


-/,.(k,»)c^.O. 


(j=n>i) 
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-icO  c      0      =    -    ik(KT   /m    )v      -  /^    (k.CO)    c      0 
ve   e  ^       e      e      e  2e  ve   e 
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(BjO    +B.p.)Q    /T       +   OL_p     (n    /N    )+   a;.p.(n./N.) 
^'n'n^i'^i'e      e  nn^n      n  I'l^x      i 

3     \^P  JT      3    3      l^T^      0      3  3      r J-  3      3 


(2.18) 


( j=n,i,e) 

The   frequency   CO  is    understood    to   be   real    and   positive   with    the 
waveform   taken   as      exp  (  ikx-icot )  .      This    means    that   a    sign    difference 
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between    the   real   and    imaginary    components    of    the   wave   propagation 
constant    k    corresponds    to  amplification.      Amplification    occurs    when 
the  amplification   mechanism   dominates    and   net    growth    of   the   wave 
occurs. 

In    the   momentum   transfer    equations    (2.15),    the   term   /    .(k,CO)v, 
represents    the    viscous    stress,  in    the    j-th    fluid   component.      A 
similar    term    /     .(k,CO)0.    accounts    for    heat    conduction    in    the    j-th 
component,    where   9.    is    the    temperature  perturbation.      Poisson's 
equation    is    used   to   relate    the   electric    field  perturbation    to   the 
density   perturbation. 

Ingar d   explicity   assumed   that    inelastic    collisions    and   ion- 
neutral    energy    transfer    could  be    ignored.      Other    assumptions   are 
implicit    in    the   analysiG    because    of   noglected   terms:    the   d.c. 
electric    field    is    not   present    in    the   energy    balance    equation,    and 
arabipolar    diffusion    to    the   walls    is    neglected. 

By   algebraic   methods,    the    equations    were   reduced   to   a    set    of 
three   homogeneous    equations    for    the    component    velocity   perturbations, 
written   here    in   matrix    form. 

ABC 


D        E         F 


G         H         I 


=    0 


(2.19) 


Approximate   expressions    for    these   co-efficients   were   given   by 
Ingard  and   Schulz,    in    the   main    body    of   their    paper,    as    being 
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consistent  with  the  acoustic  mode  of  motion  in  a  gas  with  a  low 
degree  of  ionization.   These  equations  are: 
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+  7- ;  B  =-    -    iCO    .  /CO;     C 

CO         y  ni 


ne       1 
-    1  + 


D  ^    -    iu).    /CO;     E 
m 


nn 


CO.  /       ,        CO.    +C0.    \ 

1  ./CO  m      ie\      ^ 

-f'*  M— *— ir-j'  ^ 

CO  \     XI  / 


CO      cor  ' 


CO. 

~co 


CO. 
le 

CO 


;     (2.20) 
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Using    the    co-efficients    given    in    equation    (2.20),     Ingard  and 
Schuiz   made    additional    approximations    and   expressed   the   dispersion 
relation    in    the    form: 

/      \o    r      /n 
2^  /  CO    \^    I       '    '" 


\^n)  ['-{-^^-^jt:? 


where 


-,J_        ^n      6^        1  1 

^  ^^n    '     '^     1.5^    ^^1.6^ 


n   =  (CO  .-*-  CO    )    5  =  J-Is!!!^  . 
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n      n      n      n 
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T'      P 


(2.21) 


+    (C0r]/2pr^)    +    (OJK   /2/)c   r^)       . 


It  was  stressed  that  this  expression  was  valid  only  if  the 
second  and  third  terms  in  the  brackets  were  small  with  respect 
to  unity  . 
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In  the  above  expressions,  the  adiabatic  speed  of  sound  for 
the  j-th  component,  c.,  is  defined  as: 

c  .  =  (y  .KT  ./m  .)   , 
J      J   J  J 

where:   y.    =    c     ./c     .,    the   ratio   of   specific    heats;    K,    the   Boltzman 
J  pj       vj' 

constant-         T.    and  m  .    are    the   temperatures    and  mass    of    the    j-th 
D  D 

component.       The   term  CO..(jV-^)    is    a    constant    of  proportionality 

v;hich    is    related   to   the   effective    collision    frequency    of   a 

particle    of    the   j-th    type  moving    through   the    fluid   component   i 

with   Nff    particles   per    unit    volume. 

The  Kirchoff   wave  attenuation    caused  by    viscosity   and   heat 

conduction,    which   was    introduced   through    the    terms    /,      and    f^ 
'  In  ' 2n 

in   a    neutral    gas    is    taken    into  account    by   use   of    the    characteristic 

frequency   CO      :     defined  by 
nn 


03 


CO 


nn 


COT 


~    (/         +   ^^    /       )     . 
CO     ^^In  y       ■'2n^ 


(2.22) 


The    characteristic   frequency    for    ions,    60..,    is    defined    in   a 

similar    manner.      The    ion    and    electron   plasma    frequencies    are 

denoted   by   CO.    and  CO    .      The    term   1/coT  expresses    the   effect    of    the 
1  e  ^ 

heating    of   the   neutrals    by    the   electrons. 

In    the   analysis    of    their    results,    Ingard  and  Schulz   point    out 

that   wave   amplification   will    occur    if   the    terra   in    the    imaginary 

2 
part    of   k    ,    Eq.    (2.21),    involving   f    dominates.      They    give  a    simple 

argument    and   develop    a    criterion    for    amplif icatioii    involving    the 

electron-neutral    collision    frequency  CO 

^  en 


(2.23) 


.,,,.,      r  CO 
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They  also  point  out  that  contrary  to  the  result  obtained  by 
Ingard  [6]  earlier,  this  more  general  treatment  indicates  that  a 
decrease  in  neutral  gas  temperature  does  not  always  increase  the 
amplification.   They  show  that  the  optimum  temperature  can  be  close 
to  room  temperature,  based  on  the  approximate  expression  Eq.  (2.20), 

Solving  equation  (2.19)  for  the  ion  and  electron   perturbation 

velocities  in  terms  of  the  neutral  velocity  and  using  approximations, 

it  is  shown  that  the  assumption  often  used,  that  v  ,  v.  and  v   are 

■^  e         1  n 

all    the    same    in    the  acoustic    mode,    is    not    generally    valid.       It    is 
true    only    for    frequencies    below   some   characteristic    frequency.      At 
higher    frequencies,    the    ion    and   electron    velocities    decrease    in 
amplitude  and  are   brought    out    of  phase   with    the   neutrals. 

Maccnc    observation    of  ^n   amnl i f i cat  ion    effect   at    increased 
temperatures    does    not    contradict    the    Ingar d-Schulz    development; 
the   observations    of   Yatsui   et.al.    do.      Since   Ingar d-Schulz    did 
not    consider    directional    effects,    and    did   not    set    up    the    equations 
so  as    to    include  any    drift    motion,    their    results    could   not    have 
predicted   Yatsui's    results.      The   question   arises    whether    the 
directional    effects    are  present   at    lower    frequencies    and  higher 
pressures.       In    the   Yatsui    experiment,    the    frequencies    were   high 
and   the  pressures    low,    resulting    in    conditions    v/her  e    the   wave- 
lengths   were    on    the    order    of   the   mean   free  paths. 
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III.   NUMERICAL  EVALUATION  OF  INGARD-SCHULZ  EQUATIONS 

The  despersion  relations  developed  by  Ingard  and  Schulz  predict 
a  change  in  attenuation,  possibly  even  amplification  for  an  acoustic 
wave  in  a  weakly  ionized  gas.   To  improve  the  chances  of  observing 
what  might  be  a  small  effect,  it  was  felt  that  the  optimum  con- 
ditions must  be  determined  and  the  experiment  conducted  under 
similar  conditions.   "Conditions"  here  refers  to  the  values  of 
•controlable  discharge  parameters. 

The  purpose  of  the  experiment  was  to  measure  a  change  in  the 
attenuation  constant.   When  the  terms  in  Ingard  and  Schultz's 
expressions  were  analyzed,  the  obvious  variables  were  the  micro- 
scopic nAr;^rneters  :  neutral  narticle  densitv  N  .  electron  density  N  , 
r         L  1  '   n '  e 

neutral  temperature  T  ,  and  electron  temoerature  T   and  the  angular 

^  n  '■  e 

frequency  of  the  acoustic  wave  CO.   The  problem  of  determining  the 
optimum  conditions  was  thus  the  mathematical  problem  of  minimizing 
a  general  function  of  five  variable?  with  the  variables  constrained 
to  a  certain  region,  the  region  of  experimental  feasibility.   This 
is  best  done  with  a  computer . 

Since  as  presented,   the  Ingar d-Schulz  equations  were  not 
suited  to  such  an  analysis,  the  process  was  one  of  analysis  and 
synthesis.   Analysis  was  needed  to  express  the  terms  as  relations 
involving  only  the  above  variables  and  physical  constants  of  the 
gas  such  as  cross  section  or  molecular  diameter.   These  terms  were 
then  evaluated  numerically  and  combined  with  each  other  to  compute 
the  predicted  value  of  the  attenuation  constant  for  the  input 
values  of  the  variables. 
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Using  a  computer,  the  equations  were  evaluated  for  a  wide 
range  of  values  of  the  input  variables.   This  provided  a  more 
thorough  evaluation  of  the  Ingar d-Schulz  equations  than  they  had 
made.   The  results  showed  that  their  approximate  expression, 
equation  (2.21)  predicted  net  amplification  over  a  large  region  of 
variable  space  although  the  more  complete  expressions  never  did, 
even  when  Ingar d's  criterion  for  the  validity  of  the  approximation 
was  satisfied. 

The  process  of  analysis  is  described  in  Section  IIIA  following. 
The  procedure  vised  in  evaluating  the  terms  and  determing  the 
optimum  values  for  the  variables  is  described  in  section  II IB.   The 
results  of  the  evaluation  are  presented  and  discussed  in  section 
lilC.  . 

It  should  be  pointed  out  that  this  section  discusses  the 
evaluation  of  the  Inga c d-Schulz  equations  before  the  experiment 
was  conducted.   Afterwards,  the  program  was  used  to  compare  the 
predicted  values  of  the  attenuation  constant  with  the  observed 
values  by  substituting  known  or  estimated  values  for  the  variables. 

A.   ANALYSIS 

1 .   Obtaining  the  Attenuation  Factor 

The    system   of   homogenous    equations    (2.19)    has    a    non-trivial 
solution    only    if   the    co-efficient    matrix    is    singular,    i.e.    its 
determinant    is    zero  [341.      We   can    therefore  write: 

lIco-efficients|l    =   A(El-HF) +B(GF-DI ) +C(  DH-GE)    =    0  (3.1) 


38 


We    define   the   complex    quantity    n    =    ck/co.       It    can    be   considered 

2 
an   acoustic    index    of  refraction.      Since    in   Eq. (2.20)    n      appears 

only    in    the   co-efficient   A,    we  may    solve   for    it    by    defining   A.,    as: 

/     .  CO    .+0)     \ 

2         ,       .    /    CO  ni      ne\ 

A^=   A-^n       =    1    h-x    f  _    .   ^^—j 

and  substituting  this  into  equation  (3.1).   After  a  rearrangement 

2 

of  terms,  n   is  found  to  be: 

2       ,  CiDH-BG)  ^  B(GF-DI) 

1  EI-HF       EI-HF  ^^  ' 

2 
Now   n      can   be    written    in    the    form 

2 
n      =  \+iii  (X  ,a    r  eal) 

and   the    desired  result    is    t.hp>   wave  propagations    constant   k    in    the 

from: 

k   =  0!  +  ip,  {a.,^    real)  (3-3) 

so   that    the  problem   reduces    to   solving    for    the   square   root    of   a 
complex    number.       In   particular,    since  OL    is    defined   to   be  positive, 

the   root    of    interest    is    the    one   for    which    f>    has    the    same    sign   as   \i> 

2 
This    means    that    if   n      is   written    in  polar    form: 

2  I     2|     iO 
n      =    I  n    !  e       , 

the  solution  is  the  root  at  Q/2,  i.e.  the  principal  root.   The 

attenuation  factor,  p,  the  ultimately  desired  result,  is  given  by 

p  =  ^  Im(n^)  (3.4) 

2 

where  n   is  the  principal  root  of  n  . 
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2  .       Acoustic   Losses 

An    acoustic   wave    is    damped   through    two  physically    different 
mechanisms,    viscosity   and    thermal    conductivity.      The    damping    can 
also   be   distinguished  by    the    location   where   the   losses    occur,    i.e. 
at    the  boundaries    or    in.    the  bulk    of    the   gas.      The    latter    charac- 
terization   has    the   advantage    that    it    divides    the   losses    into   two 
terms    of    different    frequency    dependence.      When    the    losses    are    high, 
the   two   effects   are    not    simply   additive   but    instead  are    coupled 
together    in    the    differential    equations    for    an   acoustic   wave.      As 
this    would    only    be   the    case   for    viscous    fluids    or    for    tubes    of 
small    diameter^  neither    of   which   are    treated   here,    we   shall    con- 
sider   the    effects    to   be   additive.      Since    the   mean   free  path  and 
therefore    the   boundary    layer    is    inversely   proportional    to   densixy, 
at    sufficiently    low   pressures    the   assumption    of   additive    effects 
may    have    to   be  re-examined. 

A  plane   wave    traveling    in   a    duct    attenuates    expontent ia lly 
with    distance   with   attenuation    co-efficient    "a".      With   small    losses, 
a    is    the    sum   of   the   bulk   attenuation    co-efficient    a      and   the   wall 
attenuation    co-efficient   a    .       For    a    circular    pipe,    these   co- 
efficients    are   given    by   Morse   and   Ingard   [35,p989]    as: 
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where  K,  the  thermal  conductivity,  is  given  by 


and 


^  pc        .        .. 

— ' —       IS    the   co- 


efficient    of   viscosity. 


The  mean    free  path  Z    is 


i    = 


nJTttd  N 


(3.6) 


In  these  expressions,  p    is  the  mass  density,  c  the  adiabatic  speed 

of  sound  and  D  the  atonic  diameter.   When  these  expressions  are 

used  for  a  three  component  fluid,  c  will  be  written  c  ;  however 

there  is  no  difference  between  the  two,  c   ^  c. 

n 

Sii^ce  thp  pnergy  of  a  wave  is  given  by  the  square  of  its 

amplitude,  the  energy  of  a  wave  attenuates  with  co-efficient  2a. 

With  a  velocity  of  c  ,  the  characteristic  time  for  energy  dis- 

n 

sipation    by    the   neutral    component,    T    ,     is    related   to   the   atten- 
uation   co-efficient   a    by: 


=   ^ac    . 


(3.7) 


Therefore,    for    the   small    damping   approximation    of  a    plane   wave 

in   a    circular    pipe,    the    characteristic    time    for    energy    dissipation 

3 
is    given   by    : 
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The   expression    does    not    exactly    agree   with    the    expression    in 
Ref .    7   but    it    is    consistent    with   the    derivation    there  and  agrees 
with    the   same    expression    in    Morse   and   Ingard  [35].    The   difference 
may   be   due   to  an    error    in    typesetting. 
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Substituting    for  JC   and  Tj  ,    using  y    =    5/3    (a    good  approximation   for 
a    noble   gas),    the   above   expression    simplifies    to: 


\ 


2.4   00     -   .- 

n  ' 


J  2.4  co^c^ 


(3.9) 


This    term   enters    the   dispersion   relation   as   CO/W        defined   by 

nn 

equation    (2.22), 
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3 .       Zero    Ionization   Limit 

In    order    to  be   consistent,    the   attenuation   predicted  by 
Ingard's    equations,    for    the    limiting    case   of  zero    ionization, 
should   be   comparable   with   previous    results    for    acoustic   atten- 
uation.      For    the    st?(ndard    of   comparTson    we   shall    use    the 
attenuation    co-efficient    "a"    which   was    used  as    an    input    in    these 
equations . 


a    =   a      +   a 
W          B 

1                        ?          - 

...            (J 

c 

n 

j2.4  CO  — 

'                   c 

n 

(3.10) 


For  zero  ionization,  equations  (2.21)  and  (3-2)  reduce 

2 
to  the  same  expression  for  n  , 


^  ^  ^  ^r  • 


(3.11) 


Now,    both    equations    (2.21)    and    (3.2)    were   based   on 
approximations.      When    the   expressions    for    the    co-efficients, 
Eq.    (2.20),    were  presented,     Ingard  and   Schulz    had  already   made 
several    approximations.      Going   back   to    the   full    expressions    for 
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the    co-efficients,    which   are    contained    in    the   appendix    of    Ingard 
and   Schulz's   paper,    an    examination    shows    that   when    the    ionization 
goes    to   zero,    only    one    independent    co-efficient   A^   remains,    given 
by: 


A      =    1    - 
o 


c^^k 

.    CO 


1- 


r-1    i  ^'2n/^ 


y     i  +  i  ;■     /CO 

2n 


+    1 


^n 


OJ 


Even    though    there    is    only    one    co-efficient,    it   must    equal 

zero    for    the    equations    to   have  a    non-trivial    solution    for    V    . 

2         ^'^n 
.Setting   A      equal    to  zero,    and   solving    for    n      = ,    we    obtain 


00 


_  l  +  i    /^     /CO 

2  In 


1      r-1^2n/'^        ' 


(3.12) 


where  /'.   and  /^   are  the  viscous  and  thermal  conduction  com- 
in       2n 

ponents  of  the  acoustic  losses  given  earlier  in  section  IIA. 

They  are  hidden  in  the  expression  for  the  attenuation  co-efficient 

The  attenuation  time  T   can  be  expressed  in  terms  of  these  com- 

n  ^ 


ponents  as 


T  ^  In     y   ^  2n 

n 


Using  the  same  expressions  and  approximations  as  before 
they  can  be  written  as: 


xfco^  I  -.^\jj2.4  COic 


/ 


2n 


5     2  i        5 
f  CO  -  +  -^  \ 
3    c   3 


J  2.4  ^-oic 


since   there   were   fewer    approximations    involved   in 
Eq.     (3.12),    the    value   of    [3    obtained    is    a    stronge).'    comparison    test 
then    the    use    of    equation    (3.11)-      Using   both   equations    (3.11)    and 
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(3- 12),  p  was  calculated  for  several  different  values  of  N  ,  T 
and  CO.   The  results  were  compared  with  the  value  of  "a"  calculated 
for  the  same  input  values  using  equation  (3.10).   In  all  cases,  the 
results  agreed  within  less  than  1%.       It  is  therefore  concluded  that 
in  the  limit  of  zero  ionization,  the  Ingar d-Schulz  dispersion 
relations  predict  the  correct  attenuation. 
^ •   Neutral  Heating  by  Electrons 

The  effect  of  energy  transfer  to  the  neutrals  by  the 
electrons  is  expressed  through  the  term  1/WT.   Ingard  and  Schulz 
define  T  by  the  relation: 

7  =  (y-1)  ^  .  (3.13) 

c 


TVi 


he  rate  of  energy  tran'=;fer  per  unit  mass.  H,  is  defined  by 


8  V.  [""^2  /kt   >>  3/2 
H  =/  ~W-^W  — r  N    cr  .  (3.14) 


3f,  M  ,  ■■      - 

3  /        m     I  ni      /  e      en 

n     /     \     n/     \        e  j 

where  CT       ,    the  average    cross    section    for    energy    transfer,    is 
en 

defined  as : 

^      X'  1  dcr  2  4 

/  3/,  r>\       en    d(cos0)f(v    )v   dv    . 

(3.15) 

The   integral    over    cos    9    is    the    total    cross    section 
(velocity    dependent)    for    momentum    transfer,    q    ,    as    defined  by 
Present    [36].       It    is    given    by: 


4       .       . 
This    IS    the   correct    expression   for    O       ,    vice   the    expression 
r    1  on 

given    m    L7J    as    confirmed  by   private   communication   with   M.    Schulz 

(due   to   a    typographical    error). 
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da 
'^D^^'e^  "/   2n(l-cos  9)   ^-^  d(cos  0)        (3.I6) 

-7  —  1 


A  Maxwellian  distribution  of  electron  velocities  is  a 
good  approximation  to  the  actual  distribution  in  a  gas  discharge 
only  at  extremely  low  values  of  E/p  [39J.   It  is  recognized  that 
E/p  is  not  low  in  the  usual  case  as  evidenced  by  figures  125  and 
126  of  von  Engle  l4o].   However  the  Maxwellian  distribution  is 
computationally  convenient  to  use  and  is  valid  as  a  first  approxi- 
mation to  the  real  distribution  function.   It  is  felt  that  the 
itn'^'r  o vsTTicnt  in  accurac"  resultino  f  r  on  the  use  of  '\    b'^tter 
approximation  such  as  the  Druyvesteyn  distribution  does  not  justify 
the  increased  computational  difficulties.   The  value  obtained  for 

0"    under  this  approximation  is  smaller  than  the  actual  value 
en 

because  the  approximation  indicates  a  greater  number  of  high 
energy  electrons  than  actually  exist  and  the  cross  section  de- 
creases at  higher  electron  energies.   Substituting  the  expression 
for  the  Maxwellian  distribution  into  Eq.  (3.1?)  gives: 


(J         =   v^q  (v  )e   2KT    , 

en   \  m   /    o  j    e  D^  e        e   dv 

The  cross  section  for  momentum  transfer  is  available  in 
terms  of  electron  energy  expressed  in  electron  volts.   By  defining 


45 


E      -  KT   /e,    the   energy    equivalent    election    temperature,    and 

effecting   a    change   of   variable,    0"         can   be   expressed    in    the 

en 


simple  form 


00 

J  o 


00 

^en=ti       qD(''ET)>'V''dR  (3.18) 


2 
m   V 

e   e  .  . 

where   R   =        ■  is    the    variable    of    integration    by   which    the 

e 

average   cross    section    for    energy    transfer    can   be   evaluated   from 

the   momentum   transfer    cross    section    data. 

5.      Collisional    Losses 

If   we   refer    back   to    equation    (2.2),    we   see    that    the   terms 

involving    the   constants    of  proportionality   CO        arise   from   the 

collision    integral    for    momentum   transfer.       The   collision    integral 

can    be   written: 

In   order    to   give  a    numerical    value   to  CO      ,    it   must    be   expressible 
in    terms    of   the   variables    and    gas   properties.      To   do    that,    its 
definition   must    be   understood.       Since   the   momentum   transferred 
from   species    1    to   2    through   collisions    is    equal   and    opposite    to 
that    transferred   from   2    to    1    through   collisions,    then   we  can 
equate: 

This    shows    that    the    constants   CO.,     are   related   to  CO,    .    as: 

Jk  kj 

N        m 

CO..     =  -~     —  CO,    .    .  (3.19) 

ik         N  .      ra .       ki  ^  ' 

3         3 
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Notice  that  due  to  the  mass  ratio  factor,  CO.,  and  60,  .  cannot  both 

jk  kj 

be   collision   frequencies.       If    in   Eq.     (3.19)    the   j-th   component    is 
taken    to   be   the   heavier    component   and    the   k-th    component    the    ion 
or    electron    component,    the   resulting    expressions    will   be   con- 
sistent  with    the   relations    given   by    Ingevrd  and   Schulz.       In 
addition,    when    the   heavier    component    is    the   neutral    gas,    this    is 
consistent   with    the  picture    of  CO        being    the   collision   frequency 
of   a   particle    of    the   j-th    type  moving    through    the   k-th   fluid 
component    since    the  neutrals    are  much   more   numerous.       We    can    there- 
fore   develop   analytical    expressions    for    CO      ,    CO.      and  CO    ..      The 
f  -^  ^  en         in  ex 

other    terms    will    bo  related   by    Eq.     (3.19). 
a.       Ion-Neutral    Collision   Frequency 

If    it    is    assumed    that    the    ion-neutral    collision    cross 
section    is    essentially    the   same   as    for    a    neutral-neutral    collision, 
then    it    follows    that    the    ion-neutral    collision    frequency    is    the 
same   as    the   neutral-neutral    collision   frequency.      The   neutral- 
neutral    collision   frequency    is    very    closely   approximated   by    the 
elementary    kinetic    theory    result    for    a    hard   sphere    ideal    gas.      The. 
ion-neutral    collision    frequency    can    therefore   be    expressed  by: 

T    c 


^in        i        \yu    i       ' 


(3.20) 


b.       Electron-Neutral    Collision    Frequency 

In   order    to    obtain   an    expression    for    the   electron- 
neutral    collision    f r equency,  we   must    start    with    the    integral   form   of 
the   expression   given    by   Present    [36]    for    the   collision    rate   for    two 
intermixed   species. 
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EN 


J  G   J         n 


c3Cr(|  V    -V    I  ) 


e      n 
"dTT 


!v    -V    |f     (v    )/     (V    ) 
e      n      n      n      e      e 


For    electron-neutral    collisions,       v    -v    \^   v      and    the    integrals 

'     '     e      n'  e 

can   be   separated.      The    frequency    of   collision    of   an    electron   with   a 

neutral   particle,    (j^      ,    is    0^^  /N    .      Substituting   also   the    expression 
^  en  EN      e  ^  ^ 

for    the   Maxwellian    distribution    for    / ( v ) ,    the   collision    frequency 

can   be  written:  „ 

3/2     /    CO 


CO 


en 


mv 
2KT 


V      Cr(  V    )e 
e      ^    e 


dv 


Finally    this    can    be    expressed    in   a    form    similar    to   the    expression 


for   C         as 
en 


LO 


en 


2KT 


n     \|     TTm 

'  e      -^o 


/" 


c.      The   Electron-Ion   Collision    Frequency 

The    electron-ion    collision    frequency    for    momentum 

transfer    CO    .    is    given   by   Tanenbaum    [4l]    as: 
ei 


CO    .     =    3.62    X    lO"-"-^    N.    T 
ei  1      e 


-3/2 


^   -A.     (N.     in    cm"^)       (3-22) 


where  -A.   =    12n  r^  N      is    the   Coulomb     cutoff    parameter 

D     e 

and  r     ,    the    Debye   length    is    given    by: 


D 
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Expressions  for  the  plasma  frequencies  for  electrons 


60  ,  and  for  ions  CO.  , 
e  1 


are   well    known   and  present    no   computational    difficulty. 
6 .       Ion-Acoustic    Losses 

Losses    caused  by    viscosity    and    thermal    conductivity    in 
the    ion    fluid   component    are   represented  by   the    term  CO/co.  .    which 


•is    defined  by 


y  .-i 


CO.  .  CO      V    rl         y  .  2i/ 


The    quantity   /    .v.    represents    viscous    stress    in    the    ion    fluid.      The 
other    term    f ^ • ,    when   multiplied   by    © . ,    the    ion    temperature  per- 
turbation,   accounts    for    heat    conduction    in    the    ion    fluid.      Adapiting 

a   heuristic    approach,    assume   that   CO/co.  .    has    the   same    form   as 

11 

CO/CO        given    by    equation    (3.9).       Under    such   an   assumption,    it    can 
be   expressed  as: 
2 


^"    =  J-  =    JTT^  ,0^  _i  ^  J.  I  ^17:1;  ^ 


CO.  .  T 

11  1 


Si 
c 


3r 


.  c  . 
1    1 


(3.23) 


where  the  subscripts  indicate  that  these  are  ion  terms. 

It  is  not  unreasonable  to  assume  that,  since  their  masses 
are  essentially  the  same  and  their  temperatures  nearly  the  same,  the 
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ion  sound  speed  c.  is  the  same  as  the  neutral  sound  speed  c  . 
However  the  mean  free  path  for  ions  H  . ,    must  be  unlike  that  for 
neutrals,  since  one  must  consider  the  effect  of  Coulomb  inter- 
actions.  Wu  [42]  expresses  the  mean  free  path  for  Coulomb 
interactions,  which  have  been  corrected  for  Debye  screening,  as: 

L.    =     ' — ~- ,  (3.24) 

N      r      Sn(x/t    ) 
e      m      ^    P      m 

2 

where      r       =  -e — -       is    the  mean    distance   of   closest   approach   for 

m         4ttg       KT  ^  ^ 

o 

two    like   charges    q.      Wu   refers    to   T   as    the   mean   plasma    temperature 
Since  we  are    considering    ion    losses,    the   temperature    to   be  used    in 
the    expression    for    r       is    the    ion    temperature.       However    the    Debye 
length    is    determined   by    the    electrons    due   to    their    higher    tem- 
perature and   greater    mobility. 

Because   the   mean   free  path    of    ions    in   a   weakly    ionized 
gas    is    limited  not   by    Coulomb    interations    but    through    collisions 
with   neutrals,    this    expression   alone   would   give   erroneous    results. 
The  mean    free  path,    \,    can   be    expressed    in    terms    of    the   mean 
velocity    v   and   collision   frequency   v    as 

V 
In   a   mixture   of  gaseous    components,    the    total    collision 
frequency    of   one    component    is    the   sum    of    the  partial    collision 
frequencies.      This    leads    to   the   result: 

1  111 

+     r + 


X    ^^  "  X .      "  X  .  .       X  .      • 

eff  le  11  m 
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Iqnorinq  lA  .   and  making  the  approximation  that  X  .       =  X   ,  the 

^       ^      le  ^         r-x-  ^j^      j^j^ 

mean  free  path  for  ions  can  be  written  as: 


B.   COMPUTATION  AND  EVALUATION 

1.   Defining  the  Limits  of  the  Variables 

Since  the  goal  of  the  analysis  is  a  prediction  of  results 
for  a  given  experimental  arrangement,  it  would  have  been  desirable 
to  formulate  the  problem  in  terms  of  filling  pressure  p  ,  filling 
temperature  T  j  tube  radius  r,  frequency  v  and  current!  .   Unfor- 
tunately the  conversion  of  p  ,  T  ,  and  I   to  the  variables  T  ,  T 

''^o'       o  o  en 

and  N      is    not   a    trivial    problem.       For    that    reason,    the    independent 

variables    chosen   were:    neutral    density   N    ,    neutral    temperature   T    , 

n  n 

electron  density  N  ,  electron  temperature  T   and  angular  frequency 
CO.   The  limits  of  these  variables  were  then  chosen  to  correspond 
to  the  values  expected  in  a  feasible  experiment. 

A  practical  range  of  values  for  the  tube  radius  is  from 
.5  cm  to  10  cm.   This  value  was  fixed  at  2  cm  except  when  investi- 
gating its  effects. 

The  lower  limit  of  the  angular  frequency  was  fixed  at 
100  radians/second  because  it  was  experimentally  possible  and 
numerically  convenient.   The  assumption  of  plane  waves  was  implicit 
in  the  formulation  of  the  problem  by  Ingard  and  Schulz.   From  an 
experimental  point  of  view,  in  order  to  avoid  the  effect  of  higher 
modes  propagating  in  a  tube,  this  sets  an  upper  limit  for  (x)  of  the 
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order  of  2nc/r  [  35] •   However  most  pulse  experiments  involve 
higher  frequencies  and  Yatsui,  et.al.  [8,9]  worked  with  50 >  100 
and  200  kilocycle  continuous  and  burst  sineisoidal  signals.   The 
limit  of  validity  of  the  continuum  description  is  when  the  wave- 
length is  of  the  order  of  the  mean  free  path.   This  corresponds  to 
about  10   radians/sec.  at  room  temperature  and  pressures  around 

1  torr  and  was  selected  as  the  upper  limit  for  CO. 

l4   -3 
The  neutral  density  limits  were  set  at  3*5  x^  10   cm 

and  1 . 5  X  10   cdi   ,  corresponding  approximately  to  filling  pres- 
sures of  .01  torr  and  50  torr  at  room  temperature  respectively. 
The  electron-density  limits  relative  to  the  neutral 
density  were  selected  arbitrarily,  with  a  logical  basis,  as 

10~  N   to  lO"  N   v'ith  ^.bsoT-'te  limi+c  of  10  nm"   to  l.S  Y  10   rm"  . 
n  n 

For  currents  of  milliamperes  at  low  pressures  around  1  torr,  the 
degree  of  ionization  is  about  10    L40,43].   A  lower  limit  two 
orders  of  magnitude  below  this  is  reasonable  and  a  higher  limit  of 
three  orders  of  magnitude  greater  would  correspond  to  ampere  cur- 
rents.  These  relative  limits  were  tempered  by  the  current  densities 
which  they  implied  for  large  pressure  changes. 

The  neutral  gas  temperature  is  not  truly  independent 
but  actually  depends  on  the  discharge  current,  pressure,  tube 
radius,  and  external  cooling  effects.   For  the  purpose  of  this 
analysis,  it  is  sufficient  to  assume  the  temperature  can  range  from 
the  boiling  point  of  the  gas  to  the  upper  limit  of  temperature 
which  can  reasonably  be  obtained  in  an  oven,  approximately  700  Kelvin 
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The  electron  temperature  T   is  definitely  not  an  in- 
dependent variable,  being  determined  by  the  discharge  conditions 
specified  through  the  other  variables  [4o].   The  only  analytical 
expression  available  implies  that  the  electron  temperature  is  a 
function  only  of  the  gas  and  the  similarity  parameter  pR.   The 
expression  is  of  limited  validity,  falling  at  very  low  values  of 
pR  and  high  p.   Rather  than  restrict  Tg  with  an  artificial  con- 
straint, it  was  felt  that  treating  it  as  an  independent  variable 
would  yield  more  valid  results.   The  liiaits  were  chosen  to  cor- 
respond to  those  normally  encountered  e:q:5er  imen  tally  [40,44], 
10^  to  10^  Kelvin. 

It  should  be  pointed  out  that  permitting  these 
vari-'^ble?  to  mo\;e  fr(^e^y   within  their  limits  places  less  con- 
straints on  their  possible  values  than  actually  exist.   The  re- 
sult may  predict  physically  unrealizable  results. 

2 .   Computation  of  Attenuation  and  Determination  of  Minimum 

For  the  computation,  the  formulas  were  first  placed  in  the 
order  in  which  their  values  were  needed.   They  were-  then  translated 
for  computer  processing  and  a  program  written  using  the  FORTRAN  IV 
programming  language.   The  values  for  the  fundamental  constants 
are  from  [45],  data  for  Argon  are  from  Cook  [46]  except  the  atomic 
diameter  and  scattering  cross  section  for  momentum  transfer  which 
are  from  McDaniel  [39]. 

The  scattering  cross  section  data  were  taken  from  figure 
4-3-l(c).of  McDaniel,  using  a  transparent  overlay  grid  and  dividers 
to  establish  a  table  of  values  of  q   versus  electron  energy  in 
electron  volts.   A  subroutine  was  written  to  use  this  table  to 
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perforin   the    numerical    integration    to    determine  O         using    Eq.     (3-l8) 

and   to    determine  CO        using    Eq .     (3-21)  .       In    computing   CO       ,    the 

en      =»   -1   V     J  ^     ^   en 

momentum  transfer  cross  section  data  were  used  to  approximate  the 
collision  cross  section.   The  momentum  transfer  cross  section  is 
from  20%  to  50%  smaller  than  the  collision  cross  section  for  all 
values  of  electron  energy.   This  substitution  causes  the  electron 
collision  frequency  to  be  low  by  ^  25%.   If  this  was  found  to  be 
critical,  a  correction  could  be  applied.   This  procedure  was  con- 
sidered to  be  justified  by  the  increased  convenience  in  setting  up 
the  problem  and  in  reducing  the  number  of   complete  integrations 
required  by  s.    factor  of  two. 

A  set  of  fixed  values:  CO  =  lO"^,  N   =  10   ,  N   =10   , 

n  e 

o  4  .     . 

T   =  300  Iv  and  T   ~  10  K.  were  used  to  '^vo.  lua  +  e  r^  1  "■  o^n^r  oc;c;  t  on  <:;  with 
n  e 

an  electronic  calculator.   This  served  as  a  cross  check  on  the 
computer  results.   To  check  the  validity  of  the  expressions  used  to 
calculate  attenuation,  a  comparison  was  made  of  the  attenuation 
factor  for  bulk  losses  with  a  value  from  another  source.   Using 
Eqn .  (3.9)  and  (3-7),  we  can  write 

^BULK        1       / 1\  ,rr  ^ 

— 2 ^~2  {  T~  ]         -  r^  ~2  ' 

CO  2cC0       \      n  /  BULK  C 


therefore 


f  c 


Evaluated  at    standard   temperature  and  pressure,    this    yields 

m      /sec 


K=    1.89  X   10-"      -1-        -1 


5i 


2  -11 

which  agrees  with  the  value  of  p '//  ,  1 • §7  x  10    ;  for  Argon 

in  Kinsler  and  Frey  [471. 

3 .   Discussion  of  Results  of  Calculations 

At  the  point  used  for  checkout  of  the  computer  program, 
the  agreement  between  the  attenuation  constant  calculated  using 
Eq.  (3.1)  and  the  Ingar d-Schulz  approximation,  Eq.  (2.21),  was 
very  good  with  differences  on  the  order  of  a  few  parts  per  million. 
However  the  agreement  was  not  so  good  in  other  areas.   The  use  of 
Eqn.  (2.21)  yielded  predictions  of  amplification  over  a  wide  range 
of  conditions.   In  contrast,  results  using  equation  (3.1)  show  no 
amplification  within  the  entire  range  ox  the  variables,  but  merely 
an  increase  or  decrease  in  attenuation,  with  respect  to  that  with- 
out ionization.   The  amplification  predicted  by  Eq .  (2.21)  is 
therefore  considered  erroneous,  arising  from  factors  being  neglected 
in  the  approximation. 

The  minimum  attenuation  was  found  to  occur  at  the  limits  of 
the  allowed  ranges  of  the  variables.   Tliose  variables  at  their  high 
limits  were  neutral  density,  electron  density  and  electron  tempera- 
ture.  The  frequency  and  neutral  temperature  were  at  their  low 
limits . 

Figure  1  shows  the  attenuation  as  a  function  of  angular 
frequency,  displayed  on  a  semi-log  plot.   The  curve  parameter  is 
neutral  density.   Neutral  temperature,  electron  temperature  and 
degree  of  ionization  have  been  held  constant  at  their  optimum  values, 
A  more  revealing  display  is  shown  in  Figure  2  which  is  a  plot  of  the 
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Figure    1.      Predicted  Attenuation    Factor    in    Discharge 

Plasma    in    Dependence   on   Ci)  for    Various    Values 

of   Neutral    Density.  Parameters    are:    T      =    8Sk, 

T      =    105k,    and   N    /N      =    10-5.  " 

e  en 

Curve    1:    N      -■  10^5    cm"     . 
n 

Curve    2:N      =6,2xlO-^cm. 
n 

16        -3 
Curve    3:N      =3.9x10         cra-^. 
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Figure  2,  Relative  Change  in  Attenuation  due  to  Ionization 
in  Dependence  on  CO  for  Various  VAlues  of  Neutral 
Dens  ity  . 
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fractional  change  in  attenuation  due  to  ionization.   The  curve 
parameters  are  the  same  as  for  Figure  (1).   Positive  values  cor- 
respond to  an  increase  in  attenuation;  negative  values,  a  decrease. 

The  ion  perturbation  velocity  was  found  to  remain  the  same 
magnitude  as  the  neutrals  and  in  phase.   The  electron  perturbation 
velocity  was  the  same  magnitude  as  the  neutrals,  and  in  phase,  at 
low  frequencies,  and  decreased  to  zero  and  went  out  of  phase  at 
high  frequencies.   The  frequency  at  which  the  transition  began 
shifted  to  higher  frequencies  as  the  density  increased.   This  is 
displayed  in  Figure  3. 

The  results  of  the  computation  indicate  that  net  amplifi- 
cation is  a  false  prediction,  arising  from  approximations;  however, 
attenuation  is  aflected  to  a  measurable  degree.   The  laaximum  de- 
crease occurs  with  high  pressure,  low  temperature,  low  frequency 
and  a  high  electron  temperature,  highly  ionized  discharge. 
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Figure    3-       Magnitude    of    Ratio    of    Electron    to    Neutral 
Perturbation    Velocity    in    Dependence    on   CO 
for    Various    Values    of   Neutral    Density. 
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IV.   EXPERIMENTAL  ARRANGEMENT  AND  PROCEDURE 

The  original  concept  was  to  use  an  auxiliary  discharge,  with 
a  gated  sine  wave  of  current  to  produce  gated  sine  waves  of  sound, 
the  sound  amplitude  was  to  be  measured  before  and  after  passing 
through  the  main  discharge  area.   From  a  comparison  of  the  signal 
amplitudes  at  the  two  points  for  various  frequencies,  the  attenu- 
ation factor  and  its  variation  with  current  magnitude  and  polarity 
would  be  determined.   The  analysis  of  a  preliminary  experiment 
showed  that  other  factors  had  to  be  considered  and  that  a  Fourier 
analysis  of  the  signals  was  needed.   Compromising  between  the  need 
for  a  signal  with  low  frequency  components  and  one  with  all  signi- 
ficant features  appearing  before  the  reflection  of  the  signal  from 
the  tube  ends,  a  one  millisecond  current  pulse  was  selected  to 
provide  the  sound  pulse.   Although  a  two  am.pere  power  supply  was 
available  for  the  main  discharge,  the  high  temperature  of  the  main 
discharge  column  limited  the  experiment  to  half  an  ampere. 

A.   SOUND  SOURCE 

The  sound  source  was  a  small  auxiliary  discharge  far  removed 
from  the  main  discharge  and  the  two  diodes  bounding  the  main  dis- 
charge.  The  auxiliary  discharge  was  always  on  when  in  operation. 
The  quiescent  state  was  a  normal  glow  with  a  current  of  less  than 
10  ma.   The  peak  current  during  a  pulse  was  around  150  to  200  ma. 
which  corresponded  to  an  abnormal  glow.   The  current  in  the  auxiliary 
discharge  was  controlled  through  a  transistor  load.   During  pre- 
liminary work,  it  was  ascertained  that  the  sound  produced  v,-as 
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closely  related  to  the  current  in  the  auxiliary  discharge  but  not 
to  the  voltage  or  to  the  power  dissipated.   Therefore  only  the 
current  was  controlled  and  monitored. 

Although  a  pulse  was  used  for  the  final  data,  the  current  con- 
trol circuitry  was  designed  for  sinusoidal  currents,  either  con- 
tinuous or  gated.   To  avoid  extinguishing  the  discharge,  a  minimum 
current  flow  was  required  at  all  times.   This  required  a  bias  signal 
equal  to  the  sinusoidal  signal.   The  continuous  mode  was  not  opera- 
tional duo  to  rapid  deterioration  (in  times  less  than  or  of  the 
order  of  one-second)  of  the  discharge  current  when  run  at  the  high 
currents  required  for  a  good  acoustic  signal.   The  cause  was  not 
thoroughly  investigated,  but  was  believed  due  to  excessive  power 
dissipation. 

The  current  control  circuitry  is  shewn  (enclosed  in  dashed 
lines)  in  Figure  4.   It  consists  essentially  of  two  circuits  feeding 
the  base  side  of  two  transistors  operating  in  parallel.   One  cir- 
cuit consists  of  a  small  d.c.  battery  with  a  variable  resistor  for 
controlling  the  quiescent  state  of  the  discharge.   The  other  cir- 
cuit sums  two  input  signals  to  provide  the  a.c.  signal. 

Because  of  their  high  voltage  capability  (300  V),  and  their 
high  dynamic  resistance,  Motorola  MJ3202  transistors  were  used. 
Parallel  operation  was  chosen  to  give  a  greater  current  range  to 
the  discharge  while  remaining  within  the  power  rating  and  linear 
range  of  the  transistors.   Emitter  biasing  resistors  were  used  to 
balance  the  load  between  the  two  transistors.   Although  some 
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Figure  4.   Sound  Source  Circuit. 
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negative  feedback  occurred  due  to  the  al)sence  of  bypass  capacitors, 
it  did  not  have  a  noticeable  effect. 

A  General  Radio  GR1217-C  Unit  Pulse  Generator  provided  the  d.c. 
pedestal  for  the  gated  sine  wave.   The  degree  of  modulation  by  the 
gated  sine  wave  and  the  amplitude  of  the  a.c.  discharge  currents 
were  adjustable  through  variable  resistors.   For  ease  of  control, 
multiturn  potentiometers  were  used  for  c:ontrol  of  bias,  quiescent 
current  and  a.c.  current  amplitude. 

The  General  Radio  equipment  used  required  the  low  side  of  the 
transistor  circuit  to  operate  at  ground  potential.   To  minimize 
electrostatic  interaction  with  the  diode-  microphones,  the  anode  of 
the  auxiliary  discharge  was  grounded.   This  meant  that  both  ter- 
minals of  the  d.c.  power  siipply  were  floating  above  ground    Sii^ce 
the  floating  potential  varies  greatly  during  the  current  modulation, 
capacitive  coupling  of  the  terminals  to  ground  was  found  to  distort 
the  discharge  current  waveform.   The  power  supply  used  (Universal 
Electronics  Model  G3050M)  was  chosen  because  of  its  low  coupling. 
However  its  maximum  voltage  was  30^  volts.   This  did  not  permit 
using  the  full  current  capability  of  the  circuit  at  low  pressures. 
Consequently  operation  was  limited  to  pressures  above  5  torr  for 
good  acoustic  signals,  although  use  of  another  power  supply  would 
have  extended  the  operating  range.   Figure  5(a)  shows  a  typical 
discharge  current  pulse  waveform  with  a  gated  sine  wave. 

The  electrodes  for  the  auxiliary  discharge  were  two  tantalum 
cylinders'  ( 1^^  cm  long  by  2  cm  dia.)  with  faces  separated  by  3-2  cm. 
They  were  mounted  axially  in  the  tube,  supported  by  .070  inch 
diameter  KOVAR  leads.   The  leads  were  beaded  ;vith  glass  for  electri- 
cal insulation.  63 


(a) 


(b) 


Figure  5.      Sound  source  current  waveform  and  resultant 
sound.      Ambient  pressure   is    8.7   torr .      Gas 
is   Argon.      Frequency   of   sine  wave  being   gated 
is    3000   Hz.      Time   scale   is    .5  ms/cm. 


(a)    Sound  source  current 
100  ma /cm. 


Vertical  scale   is 


(b)    Sound  produced  by    discharge   current. 

Vertical   scale   is   uncalibrated.      The  second 
pulse   of  sound   is    the  reflection   from   the 
near    end  of  the   tube. 
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B.  SOUND  DETECTION 

Mobility-limited  thermionic  diodes  were  used  to  detect  the 
sound  waves.   They  were  mounted  axially,  one  on  each  side  of  the 
discharge  region.   Electrostatic  shielding  was  provided  by  a 
cylinder  of  heavy,  2.3  mil.,  copper  foil  which  was  wrapped  around 
the  tube  at  the  location  of  the  diode.   The  shield  extended  several 
inches  on  both  sides  of  the  diode  and  was  connected  to  ground. 
Since  the  cables  and  connectors  were  shielded,  the  only  part  of  the 
circuit  not  shielded  was  a  small  portion  of  feedthroughs   where 
they  penetrated  the  glass .   Also  the  diode  was  not  shielded  on  the 
ends.   For  a  more  complete  discussion  of  the  diode  and  its  sus- 
ceptibility to  noise,  see  Appendix  A.   Figure  5(b)  shows  a  typical 
diode  response  signal. 

For  the  experiment,  filament  power  for  both  diodes  was  provided 
by  the  same  power  supply,  a  Power  Designs  36IOOR  regulated  power 
supply  operating  in  the  constant  voltage  mode.   Figure  6  shows  the 
circuit  employed  for  power  and  signal  detection. 

C.  TUBE    DESIGN    AND  CONSTRUCTION 

There  are  two  features  of  the  tube  used  in  the  experiment  which 
are  worthy  of  note.   One  is  the  method  of  mounting  the  elements. 
The  other  is  the  length  of  the  tube  and  the  placement  of  the  elements 
therein. 

The  tube  elements  were  mounted  on  feedthroughs   which  did  not 
distort  the  tube  walls  significantly.   The  feedthroughs   were 
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Figure  6.   Sound  Detection  Circuit. 
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sealed  with  a  special  glassblowing  technique  v;hich  is  discussed 
in  the  appendix.   The  result  was  a  tube  which  could  be  treated  as 
an  acoustic  waveguide  for  frequencies  much  higher  than  those  used. 

To  reduce  the  effects  of  reflections  and  electrostatic  inter- 
actions, the  elements  were  physically  separated  so  that  the  acoustic 
time  of  flight  separated  the  unwanted,  from  the  wanted,  signals. 
This  required  a  half  meter  buffer  region  at  the  tube  ends  and  a  one 
meter  region  between  the  sound  source  and  the  nearest  diode. 

In  addition,  it  was  discovered  during  preliminary  experiments, 
that  there  were  positively  and  negatively  charged  particles  ex- 
tending into  the  supposedly  neutral  region  far  beyond  where  they 
were  expected.   A  series  of  experiments  established  that  they  were 
real  and  tb^t  both  typp-s  were  present.   The  buffer  region  reduced 
the  effect  but  screening  electrodes  were  used  for  further  reduction 
Since  with  cold  filaments,  any  diode  signal  was  attributable  to 
charged  particles  rather  than  sound,  the  screnning  effectiveness 
was  studied  with  cold  filaments,  in  which  case,  the  diodes  acted 
like  double  Langmuir  probes.   The  cold  filament  effect  was  strong 
in  Neon  but  barely  noticeable  in  Argon.   Although  interesting  in 
itself  because  these  signals  were  detected  at  distances  from  the 
discharge  of  30,  40  and  50  cm  in  Neon  at  pressures  -^  10  torr ,  the 
phenomenon  was  studied  only  enough  to  eliminate  its  interference. 
A  schematic  drawing  of  the  tube  appears  in  Figure  7  sind  the  lo- 
cations of  the  tube  elements  appear  in  Table  1. 

Copper-constantan  thermocouples,  of  5  niil  wire,  were  used  to 
measure  the  temperature  on  the  tube  axis  at  three  points.   For 
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Figure  7.   Schematic  Drawing  of  Discharge  Tube. 

Identification  of  Symbols  and  Locations 
are  Shown  in  Table  I. 


TABLE  I 

LOCATION 

OF  TUBE  ELEMENTS 

ITEM 

SYMBOL 

LOCATION* (cm) 

End  1  t 

-149.7 

Anode  1 

Al 

1  m  /, 

■   J-  ^  ^  .  -T 

Anode  2 

A2 

-98.7 

Thermocouple  1 

TCI 

-50.3 

Diode  1 

Dl 

0 

Screen  1 

SI 

26.3 

Cathode/Anode  1 

CI 

37.3 

Thermocouple  2 

TC2 

102 

Cathode/Anode  2 

C2 

122,2 

Screen  2 

S2 

133.4 

Thermocouple  3 

TC3 

139.3 

Diode  2 

D2 

152.5 

End  2 

203  * 

*  Location  given  is  location  of  center  of  element 

T  Tube  ends  taper  from  4  cm  to  1  cm  diameter  over  5  cm. 
Distance  refers  to  beginning  of  taper. 

+  For  calculating  tube  volume,  the  effective  length  is  3^3  cm 
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support  and  insulation,  the  leads  were  first  sealed  in  fine  twin 
capillary  glass  tubing  with  ULTEK  VACSEAL  epoxy .   The  wire  to 
extend  into  the  tube  was  left  exposed.   During  construction  of 
the  main  tube,  a  small  piece  of  tubing,  large  enough  to  hold  the 
capillary  tubing,  was  attached.   The  capillary  tube  was  sealed 
inside  of  this  tube  with  the  epoxy.   The  seal  was  good  and  the 
thermocouples  worked  well.   Such  a  method  is  not  recommended  for 
the  main  discharge  area.   This  is  because  the  temperatures  en- 
countered in  the  vicinity  of  the  discharge  may  exceed  the  maximum 
recommended  for  the  epoxy,  15O  C.   In  this  experiment,  such  heating 
caused  contamination  of  the  discharge  and  limited  operation  of  the 
main  discharge  to  currents  below  .5  ampere. 

The  pumT>ing  tram  used  for  evacuating  the  tube  consisted  of  a 
VEECO  Model  EP2A1  diffusion  pump,  rated  at  85  ^/sec,  pumping  to  a 
DUOSEAL  forepump.   A  NRC  model  720  ionization  gauge  control  with  a 
type  536-SK  ionization  gauge  was  used  to  monitor  the  pressure  above 
the  trap,  which  was  kept  filled  with  liquid  nitrogen.   A  one  inch 
line  was  used  between  the  trap  and  the  high  vacuum  isolation  valve 
(a  GRANVILLE  PHILLIPS  type  C  value)  used  to  connect  the  discharge 

tube  to  the  pumping  system.   The  ultimate  vacuum  obtained  was 

-8 
around  5  x  10    torr  as  measured  by  the  gauge. 

The  gas  filling  system  was  connected  to  the  tube  via  a  (CEC 

type  25)  high  vacuum  valve.   Pressure  was  measured  with  a  manometer 

filled  with  Octoil  S.   A  trap  located  between  the  filling  system 

and  the  tube  was  partially  filled  with  liquid  nitrogen  whenever  the 

connection  to  the  tube  was  open. 
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D.   MAIN  DISCHARGE 

In  order  to  determine  the  effect  of  polarity,  the  electrodes 

at  both  ends  of -frie  main  discharge  needed  to  serve  as  cathode  or 

anode.   Therefore  they  were  of  identical  design,  consisting  of  a 

filament  and  a  cylindrical  anode.   The  filament  was  used  as  the 

hot  cathode  for  the  main  discharge  or  as  the  hot  cathode  for  a 

"Pupps  discharge"  at  the  anode  of  the  main  discharge.   It  was  a 

\   inch  diameter,  7  turn  helix,  directly  heated  tungsten  dispenser 

cathode  designed  for  laser  work  and  produced  by  SPECTROMAT.   It  is 

rated  at  an  emission  current  of  5A/cra   at  1150  C  .   The  anode 

B 

cylindrical  electrode  which  was  mounted  concentric  with  the  fila- 
ment, was  located  so  that  the  first  two  turns  of  the  filament, 
on  the  discharge  side,  were  exposed.   This  was  done  so  that  the 
Pupps  discharge  would  tend  to  go  toward  the  leading  edge  of  the 
anode,  and  thus  eliminate  anode  oscillations.   A  small  low 
ripple  (20V,  1.5A)  D.C.  power  supply  was  used  for  the  Pupp's  dis- 
charge power  source. 

Filament  power  for  the  main  discharge  cathode  was  supplied  by 
a  KEPCO  KS-8  lOOM  power  supply.   Filament  power  for  the  Pupps  dis- 
charge cathode  v;as  supplied  by  a  POWER  DESIGNS  PD36IOOR  unit.   The 
load  for  the  main  discharge  was  two  high  power  variable  resistors 
connected  in  parallel.   The  main  discharge  power  supply  was  a 
KEPCO  77OB,  600  V  d.c.  supply. 

A  reasonably  simple,  yet  effective  system  was  used  to  elec- 
trically isolate  the  main  discharge  system  from  the  remainder  of 
the  circuitry.   The  60  cycle  input  power  for  all  units  was  isolated 
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using  transformers.   All  equipment  for  the  main  discharge  was 
mounted  on  the  same  rack,  which  contained  no  other  equipment.   The 
chassis  ground  for  all  of  this  equipment  was  connected  to  a  common 
point  which  was  then  connected  to  ground  via  a  one  megohm  resistor 
in  parallel  with  a  10  /if  capacitor.   The  a.c.  and  d.c.  potentials 
across  these  pioints  were  monitored  with  high  impedance  devices. 
This  precaution  was  taken  to  prevent  the  buildup  of  a  high  elec- 
trostatic potential  between  the  two  ground  points.   The  a.c.  current 
in  the  discharge  was  monitored  by  using  a  current  probe  with  a 
sensitivity  of  1  mv/ma.   To  prevent  a  large  potential  difference 
between  the  two  ground  systems,  arising  from  interaction  of  the 
fields,  a  potentiometer  was  used  to  vary  the  potential  of  the  main 
electrodes  with  respect  to  the  secondary  ground.   This  is  iJiiOvvu 
in  Figure  8  as  part  of  the  main  discharge  circuit. 

The  discharge  was  ignited  with  a  Tesla  coil  in  conjunction 
with  a  copper  wire  wrapped  around  the  discharge  tube  between  the 
anode  and  the  cathode.   This  was  necessary  because  of  the  high 
breakdown  voltages  required  at  the  pressures  used. 

E.   ELECTRONIC  MEASUREMENTS 

The  three  waveforms  of  interest  -  the  sound  source  current, 
and  the  two  diode  signals  -  together  with  a  trigger  signal  were 
recorded  on  a  four  channel  FM  tape  recorder  (Precision  Instrument 
Co.,  PI-6200)  for  later  analysis.   Recording  was  done  at  37.5  ips . 
Record  and  playback  specifications  are  contained  in  Appendix  C- 
The  frequency  response  did  not  appear  to  distort  the  two  diode 
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signals  but  it  did  affect  the  trigger  and  discharge  current  signals 
since  these  had  significant  high  frequency  components.   Figure  9(a) 
shows  a  typical  diode  signal.   Figure  9^b)  shows  a  typical  playback. 
of  a  diode  signal.   It  should  be  noted  that  these  figures  are  not 
derived  from  the  same  sample. 

The  trigger  signal  was  provided  by  a  General  Radio  GR  1217-C 
Unit  Pulse  Generator  operating  at  a  pulse  repetition  rate  of 
10  Hz.   The  output  signal  v;as  a  one  millisecond  rectangular  pulse 
which  was  used  to  trigger  the  scope  and  was  also  recorded.   The 
"delayed  sync"  output  was  used  to  trigger  another  unit  pulse 
generator  from  which  a  one  millisecond  pulse  was  fed  to  the  cur- 
rent control  circuit.   The  discharge  current  was  observed  by 
measuring  the  voltage  across  a  j   ohm  precibion  le^isloi. 

The  output  signals  from  the  diodes  were  amplified  with  two 
differential  amplifiers,  TEKTRO?\IX  RM122 ,  with  an  input  impedance 
of  10  megohms.   The  amplifiers  could  provide  amplification  of 
either  lOOX  or  lOOOx.   Both  output  signals  were  recorded  on  tape. 
The  diode  output  contained  very  sharp  noise  spikes  which,  after 
amplification,  adversely  affected  the  recorder  operation.   To 
eliminate  this,  a  slicing  circuit,  which  clipped  all  signals  whose 
value  exceeded  +  1.5  volts,  was  employed  betiveen  the  amplifier 
output  and  the  tape  recorder.   This  circuit  permitted  an  improvement 
in  the  signal  to  noise  ratio  since  the  noise  could  not  exceed 
1.5  volts  and  the  signal  could  be  amplified  to  a  level  just  bel 
that. 


.  ow 
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(a) 


(b) 


Figure  9.   Comparison  of  Typical  Diode 
Signal  with  Playback  of 
Diode  Signal . 

(a)  Direct  display  of  diode  signal. 

(b)  Playback  of  diode  signal.   Due 
to  triggering  problems,  the 
playback  is  not  a  playback  of 
(a) .   The  playback  does  show 
how  little  the  recording  process 
distorted  the  signal. 


74 


F.   EXPERIMENTAL  PRCX:EDURE 

After  preliminary  experimental  work  had  established  the  back- 
ground, the  progress  of  the  experiment  can  be  divided  into  five 
stages:  (1)  construction  of  the  tube,  (2)  preparation  of  the  tube, 
(3)  calibration  of  the  diodes  (4)  data  accumulation,  and  the  last 
stage,  (5)  analysis  of  results.   The  analysis  of  results  will  be 
discussed  in  the  next  chapter  . 

1 .  Construction 

Due  to  its  great  length,  the  tube  was  constructed  in  three 
sections.   After  construction  and  before  annealing,  each  section 
was  leak  tested  with  a  helium  mass-spectrometer.   The  sections 
were  annealed  in  a  vacuum  and  then  sealed  under  vacuum  until  time 
for  assembly.   The  support  system  built  for  the  tube  consisted  of 
a  long  aluminum  beam  v;ith  six  rigid  tube  supports.   This  permitted 
the  sections  to  be  rotated  while  being  sealed  together  so  that  the 
joint  was  reasonably  smooth  with  no  sags  or  large  irregularities. 
After  the  sections  were  joined  together,  the  tips  of  the  thermo- 
couple support  tubes  were  broken,  the  thermocouples  inserted  and 
sealed  with  epoxy.   The  epoxy  was  then  cured  for  two  hours  with 
heat  lamps.   After  curing,  the  tube  was  placed  in  its  final 
position  and  sealed  to  the  system  via  a  section  of  glass  tubing 
containing  a  metal  bellows. 

2 .  Preparation 

Preparation  of  the  tube  involved  bakeout   and  outgassing 
of  the  tube  and  its  elements  so  that  a  high  vacuum  could  be  ob- 
tanied  followed  by  filament  activation.   Bakeout   was  perfcri.ed 
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with  heating  tapes  at  200  C  for  over  24  hovirs  with  the  system 
under  evacuation.   The  cold  trap  was  kept  filled  with  liquid 
nitrogen  following  bakeout .   Outgassing  was  performed  by  induction 
heating,  after  bakeout.   Since  the  leads  were  KOVAR ,  vjith  a  much 
lower  melting  point  than  the  tantalum  cylinders,  outgassing  of  the 
tantalum  at  elevated  temperatures  was  done  in  short  bursts.   In 
the  case  of  the  diodes  and  main  discharge  electrodes,  the  filaments 
were  activated  first  and  then  the  anodes  were  outgassed.   Filament 
activation  consisted  of  slowly  raising  the  temperature  by  varying 
the  filament  voltage  until  the  temperature,  as  measured  with  an 
optical  pyrometer,  was  at  a  certain  value.   For  the  diode  filaments, 
this  temperature  was  1750  C  ,  (subscript  B  denotes  brightness 
temperature)  the  jl  ecouiiueuded  temperature  for  thoriatcd  filaments- 
Following  the  manufactur ers instruct  ions ,  the  main  cathode  fila- 
ments were  only  heated  to  1190  C  .   This  temperature  was  held  for 
five  minutes . 

Following  outgassing  and  filament  activation,  the  filaments 
were  aged  for  several  hours  at  their  normal  operating  temperatures. 
Then  the  main  discharge  anodes  were  cleaned  by  electron  bombard- 
ment.  This  was  accomplished  by  applying  300  volts  to  the  anode,  at 
which  time  the  anodes  turned  bright  red  under  the  absorption  of 
45  watts  of  power . 

3.   Calibration  of  the  Diodes 

After  the  tube  was  prepared,  with  all  elements  outgassed, 
the  diodes  had  to  be  calibrated.  Since  the  mobility-limited  cur- 
rent is  less  than  the  space-charge  limited  current,  and  the 
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transition  is  a  function  of  gas  pressure,  it  is  best  first  to 
attain  space  charge  limited  operation.   This  was  easily  determined 
by  preparation  of  a  log-log  plot  of  the  plate  current  versus 
plate  voltage  to  a  value  of  plate  voltage  (30  volts),  well  in 
excess  of  that  to  be  used  (13-2  volts).   This  was  found  to  be  a 
straight  line,  indicating  that  the  filament  voltage  selected 
provided  a  good  operating  point. 

Calibration  consisted  of  varying  the  pressure  in  the  tube 
and  measuring  the  plate  current  as  the  plate  voltage  was  varied. 
The  filament  voltage  was  held  constant.   As  the  pressure  increases, 
there  is  increased  filament  cooling.   It  is  not  practical  to  use 
an  optical  pyrometer  to  maintain  constant  filament  temperature  at 
all  pressures.   Since  the  power  applied  increases  as  tne  temyjera- 
ture  decreases  for  constant  voltage,  the  constant  voltage  mode  is 
semi -compensatory . 

The  filaments  were  thoriated  tungsten  and  at  first  operated 
extremely  well.   After  the  main  discharge  was  first  run,  which 
overheated  the  epoxy  seal  at  the  thermocouple,  the  diode  current 
became  noisy  and  started  to  decrease.   It  became  necessary  to 
increase  the  current  to  a  high  value  and  operate  the  diode  fila- 
ments as  pure  tungsten.   This  required  a  r ecalibrat ion . 
4 .   Data  Accumulation 

It  was  decided  to  take  data  at  three  pressures.  The  low 
pressure  was  selected  to  be  around  8  torr  because  a  larger  power 
supply  was  needed  to  obtain  a  good  current  pulse  at  pressures  much 
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lower  than  that.   The  high  pressure  was  l6  torr  because  of  the 
difficulty  in  getting  the  main  discharge  started  at  higher 
pressures.   The  third  pressure  chosen  was  12  torr. 

In  order  to  be  able  to  detect  the  existence  of  any  trends 
(in  the  data),  due  to  time  the  order  of  experiments  was  chosen 
with  a  random  number  table.   First  the  order  of  the  pressures  was 
chosen  and  then  the  order  of  the  polarity  of  the  main  discharge. 
The  order  of  runs  is  shown  in  Table  II.   Due  to  the  difficulty 
of  filling  to  the  same  pressure,  the  runs  were  not  completely  in- 
dependently ordered. 

TABLE  II 


""""^--..^^^^Polar  ity 
Pr  essur  e^^^v^^^ 

A 
L 
L 

+  1 

0 

-1 

8 

1 

2 

1 

3 

12 

3 

3 

2 

1 

16 

2 

1 

3 

2 

A  set  procedure  was  followed  when  taking  the  final  data. 
First  all  equipment  was  turned  on  and  allowed  to  warm  up.   Then 
the  tape  recorder  was  calibrated.   The  calibration  of  the  oscil- 
loscope was  checked.   For  each  pressure,  the  steps  were  the  same 
except  for  minor  variations  due  to  the  different  order  of  the 
polarities  within  each  pressure.   The  trap  in  the  filling  line 
v/as  partially  filled  with  liquid  nitrogen.   After  the  tube  was 
filled  with  Argon  at  the  desired  pressure,  the  valve  to  the  filling 
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manifold  v/as  secured.   The  diode  and  main  discharge  filaments 
were  given  plenty  of  time  to  warm  up.   Wlien  the  main  discharge 
was  on,  the  Pupps  discharge  was  adjusted  for  minimum  AC  current 
variation  in  the  main  discharge.   This  was  observed  on  the 
oscilloscope  using  the  current  probe.   The  d.c.  diode  signal  was 
measured  after  adjustment  so  that  the  amplified  a.c.  signal  did 
not  exceed  the  clipping  level  of  the  slicing  circuit.   When 
everything  was  ready,  the  tape  recorder  was  started.   While  it 
was  running,  the  temperatures  were  recorded.   Between  runs,  the 
tape  recorder  calibration  was  checked  and  the  output  zero  adjusted, 
but  the  input  gain  was  not  touched.   After  the  experiment  v;as  over, 
a  5000  Hz  tone  burst  sine  wave,  stepped  down  to  2  millivolts  peak 
to  peak  wai>  fed  to  the  amplifierc  ir:  parallel.   This  wa?  r<:>rorded 
for  later  use  as  a  calibration  aid.   After  each  set  of  runs  at  a 
single  pressure,  the  tube  was  pumped  dov"n  to  a  hard  vacuum  and 
refilled  for  the  next  set  of  runs. 
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V.   THEORETICAL  ANALYSIS 

This  section  treats  the  analysis  of  data  theoretically.   It 
develops   the  expressions  needed  to  utilize  the  data  to  obtain 
information  from  it  concerning  the  processes  under  investigation. 
In  doing  so,  it  takes  into  consideration  the  progagation  of  errors. 

A.   ENSEMBLE  AVERAGING 

Nine  experimental  runs  were  made,  each  with  different  conditions. 
An  event,  or  experiment,  consists  of  a  current  pulse,  and  the  sub- 
sequent sound  production.   An  experimental  run  consists  of  many 
events  or  replications  of  the  e>qDeriment.   The  data  recorded  for 
each  event  were  the  triggering  signal  and  the  waveforms  at  three 
points.   Ailei  convei-tjion  Lo  uiyilal  xuiiu,  eat_h  wavexOiiu  was  re- 
presented by  a  set  of  samples.   Since  a  run  consists  of  a  set  of 
events,  the  data  consists  of  a  set  of  sets  of  samples  which  shall 
be  referred  to  as  an  ensemble  of  sample  sets. 

Consider  just  one  of  the  signals.   The  samples  in  any  one  of  the 
sets  were  not  taken  randomly,  but  have  fixed  relationships  to  each 
other  which  reflects  their  time  dependence.   For  the  set  of  samples, 
the  zero  of  time  is  chosen  to  correspond  to  the  triggering  event. 
Thus,  the  n-th  sample  represents  the  value  of  the  signal  at  time 
t  =  nT  where  T  is  the  period  of  the  sampling  frequency. 

This  sample  represents  not  only  the  time  signal  at  that  time, 
but  also  something  else.   It  represents  the  value  of  the  cross 
correlation  between  the  signal  at  that  time  and  a  Dirac  delta 
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function  signal  v/hich  coincided  with  the  triggering  event.   This 
is  a  subtle  point  which  was  used  for  advantage  in  eliminating 
noise. 

The  cross  correlation  function  /(x)  of  two  real  functions 
g(x)  and  h(x)  is  given  by: 


CO 


/(x)  =       g(u-x)h(u)du 


—CD 


Let  g  be  the  Dirac  delta  function  6(t)  and  h  be  the  signal  s(t). 
Then 


CD 


/(t)    =  (  6(u-t)s(u)du 


—  oo 


.on 


and   therefore, 

/(t)    =    s(t). 

By  ensemble  averaging  the  sets  of  samples,  we  obtain  Lhe 
ensemble  average  of  / ( t )  ,  which  is  the  cross  correlation  functit 
of  the  signal  with  the  tr igger ing  event.   For  a  large  enough 
ensemble,  this  should  approach  the  time  average  of  the  cross  cor- 
relation function.   In  the  time  average,  contributions  from  non- 
correlatable  portions  of  the  signal,  such  as  those  due  to  noise, 
random  or  non-random,  will  vanish. 

Because  of  their  intimate  relation,  /(t)  and  s(t)  may  be  used 
interchangeably  except  that  it  must  be  recognized  that  in  each 
sample  set,  s(t)  has  a  different  reference  time  while  /(t)  has  the 
same  reference  (for  this  reason,  s(t)  will  be  used  for  /(t)  except 
when  the  special  features  of  /(t)  are  needed). 
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When  /(nT)  is  ensemble  averaged,  the  mean  value  and  its 
standard  deviation  are  obtained.   Since  there  is  no  reason  to 
assume  that  the  distribution  of  noise  or  fluctuations  in  value 
of  /(nT)  should  depend  on  n,  the  standard  deviation  is  considered 
a  constant.   Then,  if  /.(nT)  represents  the  value  for  the  i-th 
set  of  samples,  the  process  of  ensemble  averaging  can  be  repre- 
sented by 


1  '^ 
/(nT-)  =  ^^_^   /.(nT) 


(5.1) 


and 


a  2 

n 

1        /^ 
N-l    (^i  =  l 

02 

=  —  L        cr 

M     ,^rrl           n 

(5.2) 


where  M  is  the  number  of  points  in  the  sample.   It  is  possible 
to  consider  this  ensemble  average  /(nT)  and  the  standard  deviation 
O    as  the  expected  value  of  the  signal,  s ( nT )  and  its  standard 
deviation,  cr . 

B.   FOURIER  TRANSFORMATIONS 

When  a  harmonic  waveform  is  fed  into  a  linear  time-invariant 
device,  the  output  is  also  harmonic  but  in  general  with  a  dif- 
ferent amplitude  and  phase  [48].   This  is  true  whether  the  device 
be  a  filter,  or  transducer  or  some  other  system  with  an  input  and 
an  output.   Such  a  device  can  be  completely  specified  by  a  fre- 
quency dependent  complex  quantity  Z(v)  given  by  the  ratio  of  the 
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output  to  the  input.   If  the  input  v  (t)  is  given  by 

v^(t)  =   Ae 

and  the  output  by 

i2nvT+i0 


then 


v^{t)    =    Be- 


Z(V,  -   f  e^^ 


If  an  arbitrary  input  v  (t)  is  applied  to  the  device  and  it 
is  desired  to  calculate  the  output  v  (t),  v  (t)  is  analyzed  into 
its  spectrum;  each  spectral  component  multiplied  by  the  corre- 
sponding transfer  factor  to  obtain  the  spectrum  of  the  output,  and 
then  v  (t)  is  synthesized  from  its  spectrum.   Thus 

v^Cv)  -  Z(V)V^(V) 
and 

v^{t)    -      Z(V)V^(V)G      dv. 

(Here  and  elsewhere,  where  a  quantity  is  viewed  in  both  the 
frequency  and  time  domains,  the  upper  case  symbol  will  be  used 
for  the  frequency  domain  and  its  lower  case  for  the  time  domain.) 
If  the  process  of  convolving  g(x')  and  h(x),  is  represented  by 
g*h,  defined  by 


fix)    --   g^h  =(       g(u-x)h(u) 


du , 


and    the  process    of   Fourier    transformation    is    represented  by    ^    , 
the    determination    of   v    (t)    can   be  presented  by    the    diagram   below, 
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v^(t)    o        V^(V) 


z(t)     O    Z(V) 

i  i 

V2(t)    D    V  (V) 


(multiplication) 


In  this  experiment,  there  are  three  time  dependent  waveforms, 
i(t),  s  (t)  and  s  (t),  corresponding  to  the  sound-producing  dis- 
charge current,  the  sound  at  point  1  and  the  sound  at  point  2 
respectively.   The  object  of  the  experiment  was  to  obtain  the 
transfer  function  of  the  conversion  of  current  to  sound,  Z  (v), 
and  the  transfer  function  of  the  spatial  region  between  points  1 
and  2,  Z, „(v).   The  entire  process  can  be  represented  by  the 
diagram  belov;,  where  the  time  representation  of  the  transfer 
functions  and  the  convolution  process  has  been  omitted  for  clarity 

i(t)   r>    I(V) 


Z  (V) 
s^  ' 


SQ(t)  =>        Sq(v) 


i 

s^(t)   :>   S^(v) 


Zi2(v) 


S2(t)   :,   .S^;(V) 


8Zf 


The  application  of  Fourier  transform  theory  to  this  problem 
assumes  linearity  of  the  processes.   That  is,  if  s   (t)  is  the 
output  at  point  2  to  an  input  of  i  (t)  and  if  s   (t)  is  the 
response  to  i  (t),  then  the  response  to  i(t)  =  i  (t)+i  (t)  is 
So('t)~  s   (t)+s   (t),  regardless  of  the  choice  of  i   and  i  . 
In  the  transform  domain,  the  process  from  current  to  sound  at 
point  2  is  represented  by  the  set  of  equations: 

S^(V)  =  I(V)  Zg(v) 

s^(v)  -  Sq(v)  Zq^(v) 
s^Cv)  =  s^(v)  z^2(^) 

since  i(t),  s  (t)  and  s  (t)  are  knovvn ,  their  respective 
transforms  are  known.   Therefore  these  con^sLitute  three  equcticnc 
with  four  unknowns,  S  (v),  Z  ^(v),  Z^, (v)  and  Z  (v).   For  the  case 
where  the  main  discharge  is  on,  only  (z  „(v)l  is  desired.   This 
can  be  solved  for  easily. 


To  determine  Z  ,  when  the  discharge  is  off,  acoustic  theory 
provides  a  model  from  which  a  fourth  equation  can  be  obtained.   For 
sound  propagating  in  a  tube,  the  transfer  function  between  two 

points  is  given  by  ik(v)d^2 

Z^2(V)  =  e 

where  d    is  the  distance  between  1  and  2.   With  this  equation  and 
the  additional  assumption  that,  in  the  absence  of  a  discharge  in 
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ik(v)x^ 


the  intervening  space,  k(v)  is  the  same  throughout  the  tube,  the 
system  reduces  to  two  equations  in  two  unknowns, 

S^(V)  =  I(V)  Z^{'^)e 

ik(v)x 
S2(V)  =  I(V)  Zg(v)e 

The  propagation   constant   k(v)    is    complex, 
k(v)    =  a{v)+    ip(v) . 

Under    the    exyier imental    conditions,    boundary    losses    dominate.       In 

this    case,    P(v)    =    p     \l~^  •       I^  phase    information    is    neglected, 

meaning    that    the    delay   between    i(t),    S    (t),    and   S„(t)    is    neglected 

then    p      can   be   obtained   from   S    (v),    S„(\')    by 

/,  ..^ 

\   1 S^ ( V  )  I    / 


°      ^'2-^1  d(v-^)       U^i(^^ 


Is,(v)| 
's^""^'    ~~   |i(v)| 


\z(v) 


C.   ACOUSTIC  TRANSFER  FUNCTION 

After  ensemble  averaging,  the  two  diode  signals  are  represented 
by  sample  sets  and  standard  deviations 

s^(nT),  O^ 

n  =  1,M 
and  s^{nT),    O^ 

Since   the   standard    deviations    are   considered   time    independent,    each 
represents    a    constant    signal    in    time,    the    transform    of   ivhich,    T.  , 
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is   a    zero   frequency    component    of   the    same   magnitude   as   CT .      There- 
fore,   in    the   transform   domain,    the   two   diode   spectra    are 

s^(v.) 

and  22^^i^ 

2  2 

where  S      and  E^    affect    only    S(v      =    0). 

The   standard   deviation    or    error    in   S,(v.)    and   S^(v.)    is    that    due 

to   the    noise,    plus    that    due   to    other    sources.       (The    error    in    the 

value   of  a    quantity    such   as    S,     is    represented  by    A   preceeding    the 

symbol    for    the  primary    quantity,    e.g.    A5)    .  ) 

Let   Y(v.)    represent    the    natural    loaorithm    of   Z,„(v.). 

Then  |s    (v.)| 

Y(V.)    .^~£--i_-  (5.3) 


e 


and 

For    the    discrete   case,    the   earlier    expression    for    p       is    not    very 
useful.       Another    e>qDression   which    can   be   written    containing    p   '    is 

Y(v.)    =  -   P^(x'2-x^)v."   .  (5.5) 

This  expression  can  be  solved  for  B   at  each  value  of  v . .   However 

o    .  1 

a  better  approach  is  to  use  all  the  information  available  and 

minimize  the  error  in  6  .   If  all  values  of  Y(v.)  were  known  to 

^o  ^     i' 

the   same    degree    of  accuracy,    then    a    simple    linear    least    squares 
fit   would  provide   the   optimum   estimate   for    p    .       However,    from 
the   ensemble   averaoe   and   another    consideration    to   be   discussed 
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later,  the  error  in  Y(v.)  is  known.   A  weighted  average  will 
therefore  improve  the  estimate  for  p  .   The  values  of  Y(v . )  are 
best  weighted  with  a  weighting  function 

w  =  — ^ (5.6) 

AY(V.)^ 
where  A  is  a  normalizing  constant  given  by 

A  =  S    AY(V.)'^  • 
i  =  l      1 


D.       SOURCE  TRANSFER    FUNCTION 

When    there    is    no    discharge   between   points    1   and    2,    it    is 
assumed   that    the   transfer    function    is    the    same   throughout    the    tube, 
This   means    that    the    value   of   p      may   be   used   to   extrapolate    the 
souiid  spectruiu   observed  at   pomtii    1   and   2    to   the    oriLyiu.      There    xt 
may   be   compared  with    the   current    spectrum    of    the   sound   source 
discharge    to   obtain    the    source   transfer    function. 

The   spectral    components    of   S      and   S      projected   to    the    origin 


are   given    by: 

(3=1,2) 

If  B   is  used,  the  implication  is  that  the  estimate  of  B   is 
o  "^  o 

independent  of  the  S..   But  since  B   was  derived  from  the  estimates 

3  '  o 

for  S^  and  S..  ,  this  is  not  true.  Recall  that 


^     Po^  (^2-^l) 


S   =  e  -     -   ^   .  (5.8) 

^2 
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therefore  1 


This    can   be   substituted    into   Eq.     (5.7)    to   give 

^1 


(5.9) 


^l\  ""2-^1 


Sio(^i)    =   ^l^S^ 


and  ^_ 

/Si\x2-xi 

S20(^i)    ^   ^2(i^j 

Since   these   are   equivalent,    the   estimate   for    S,.    is 

^2  -    ^'l 


^2"^'l  ^2""^! 

Sq    =    S^      ^      -"    S2        ^  .  (5.10) 


The  standard  deviation  for  S^  is  given  by 


{{- 


^2  \^   /^^iV    /  ^1  \^  /^^2^^^  ^ 


-o-o<(.p^   -   M^?^   -^'^  ■''•"' 


The  transfer  function  of  the  source  is  obtained  with  a  simple 
division, 

s„(v  ) 

1 

Whereby  \2    Z      \  2 

E.   EFFECT  OF  THE  DISCHARGE 

In    order    to   make   a    valid   comparison    between    experimental    results 
and   theoretical   predictions,    the  predictions    should   be   based    on 
values    which    represent    the    experiment.       In    this    case,    to   corapare 
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the  Ingar d-Schulz  predictions  with  experiment,  the  effects  of 
longitudinal  and  radial  thermal  gradients  must  be  taken  into 
consideration.   Since  the  electron  temperature  and  electron  density 
are  not  known,  they  must  be  estimated.   Also  it  should  be  noted 
that  the  discharge  does  not  fill  the  region  between  the  two  diodes 
but,  using  the  values  from  Table  1,  occupies  only  55 •^'>o   of  the 
space.   When  there  is  no  discharge,  this  does  not  affect  results 
but  when  the  discharge  is  on,  the  region  is  inhoraogenous  and  the 
attenuation  constant  obtained  is  the  average  value  for  the  entire 
space  between  the  diodes. 

Several  assumptions  must  be  made  in  order  to  predict  the 
experimental  results.   The  first  is  that  there  is  no  transition 

predicted  value  for  the  effective  attenuation  constant  should  be 
the  average  of  the  value  for  no  discharge  and  that  in  the  discharge, 
where  the  values  are  weighted  by  the  lengths  of  the  two  regions, 
E>ue  to  the  nature  of  the  interface  between  the  regions,  reflections 
are  considered  negligible.   Temperature  measurements  taken  between 
the  discharge  and  diode  2  support  this  assumption.   They  showed  only 
a  one-degree  higher  temperature  than  measurements  far  removed  from 
the  discharge.   The  two  diodes  are  separated  by  a  distance  d„  and 
the  main  discharge  elements  by  d  .   Thus,  if  p    and  p    denote 
the  attenuation  constant  in  the  discharge  and  outside  of  the  dis- 
charge respectively,  the  effective  attenuations  constant  p'  is  given 

d         d 
f^O    d.  "^Ol    d„  P02 


wher  e  d   =  d   -  d  . 


0         0 
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other  assumptions  concern  the  effects  of  the  discharge  through 
thermal  gradients.   During  their  investigation  of  the  neutral 
temperature  in  medium  pressure  discharges,  Hayess,  et.al.  L49J 
observed  a  radial  temperature  gradient  under  all  discharge  con- 
ditions.  This  will  affect  the  propagating  sound.   Hayess  L 50 J 
considered  the  problem  in  conjunction  v;ith  his  experiment.   He 
concluded  that  for  frequencies  below  the  first  radial  resonant 
frequency,  only  the  fundamental  mode,  a  plane  wave  would  pro- 
pagate.  He  further  concluded  that  the  temperature  difference 
would  not  affect  the  attenuation  greatly  since  only  its  square 
root  appears  in  the  attenuation  factor.   The  cut  off  frequency 
for  this  experiment  is  about  9  kHz.   The  discharge  region  can 
be  characterized  by  two  temper  a  1  ujl  ei^  for  two  different  purpoccc. 
The  longitudinal  gradient  in  the  tube  affects  the  pressure  in 
the  tube.   For  this  purpose,  the  dischaige  can  be  characterized 
by  its  mean  temperature  T  .   Since  the  temperature  outside  of 
the  discharge,  T  ,  did  not  change,  the  tube  pressure  with  the 
discharge  on  is  given  by 


1 
p.  =   p   =   p 


o      d-  /     T  \   ^o  /    T   \ 

where  d   is  the  length  of  the  tube.   The  uncertainty  in  the 
pressure  is  given  by 


"M 


'^0)^'°    ■  , 
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Since      ■!_■  ■  =  p'  — ,    this  can  be  expressed  as 

°^M       0  T 
M 


M 


2 


Under    the   assumption    that   boundary    effects    are   the    dominant 
loss   mechanism,    the   temperature  at    the   u-all,    T    ,    is    important    for 
the    discharge    region.      This    is   because   these    losses    occur    within 
two    or    three   mean    free  paths    of   the   boundary. 

To    calculate    the  predicted  attenuation    constant   with    the 

Ingar d-Schulz    equations,    the    experimental   parameters    must    be 

converted  to    the    variables    used    in    the    computer    program:    (x),    N    , 

N    J    T    ,    and   T    .       Tlje    wall    tenipei-atuj- e    can    be    used    foi.    T     .       Tlje 
e         e  n  '  n 

neutral    density    can   be   calculated   from   p'    with   T. , ,    using    the 


express  ion 


1  8 
N      =    9.6642   X   10        p/T 
n  ^ 


-3  .  .  . 

where  p    is    m    torr    and   N      m    cm       .      The   relative   uncertainty    m 
'  n 

the    density    is    given   by 

AN  V    /^  \2   .  /_\2 

^f)  if)  if)  ■ 

since  they  were  not  measured,  estimates  must  be  used  for  the 
electron  temperature  and  the  electron  density.   With  the  similarity 
parameter  Rp  as  an  entry,  T   may  be  obtained  from  von  Engle  [4o] 
with  an  estimated  uncertainty  of  10%.   Estimating  the  electron 
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density  requires  a  more  complicated  process.   Venzke  L51J  has 
plots  which  give  E/p  as  a  function  of  i/R  with  Rp  as  a  curve 
parameter.   The  value  of  E/p  thus  obtained  can  be  used  to  obtain 
the  electron  drift  velocity  from  page  5''^5   of  MdDaniel  [39].   With 
the  assumption  of  a  uniform  current  density  across  the  tube,  the 
electron  density  is  given  by 

I 


^e    2nr  v  n   * 
o  D^e 

Of   course    it    is    known    that    the   current    density    is    not    uniform 
across    the   tvAbe   but    neither    are   other    factors    such   as    neutral 
density.      The   entire  acoustic   problem    is    treated    in   a    one-dimensional 
sense.       Because   of   the   approximations    involved,    the   electron    den- 

The   relation    of   p    ,    the   attenuation    constant,    to    p,    the 
attenuation    co-efficient,    can    be   viewed    in    two   different    ways. 
If   p    =    p     ^To"  ,    the    two  ways    are    equivalent.      One    way    is    to  write 

Po   =    P(W  =    1). 
The   other    way    is    to  write 

Pq-  P(w)/  \(aJ. 

The    latter    form    is   preferred   since  a    departure   from    the   assumption 
of  boundary    loss    domination   will    be  readily    apparent    from   a   plot 

of   p      as    a    function    of  OJ.       It    should   be   noted   that    Figure   2,    which 

Pl"^UI 
is    a   plot    of     — ,    V70uld   not    be   changed  by    switching    p      for    p. 

^UI  ° 

Since    p      should   have    no   dependence   on   to,    the   value    of   OJ  should   not 

be    important.      To   be   certain,    the    value    corresponding    to   the    highest 
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and   lowest    experimental    values    of  60  should  be  used  and   the 

dependence   cliecked. 

Since   alD.    of   the    input    quantities    are    estimates,    with   an 

uncertainty   which   can   be   treated  as    a    standard   deviation,    the 

estimate    of    B      obtained   will    be   uncertain.      The  uncertainty    in    6 
^o  '  ^  o 

is    given    by 


^n    \b^J 


The  partial    derivatives    can    be    obtained   with    the   computer    program 
by    the   expression 

be  e    (x.=X.+AX.)    -    B    (x.=X.) 

bx.  AX. 

If  Ap   is  particularly  sensitive  to  any  of  the  variables,  the 
process  of  evaluation  should  be  reviewed  to  see  if  the  uncertainty 
can  be  reduced. 
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VI .   TREATMENT  OF  D\TA 

The  choice  of  data  handling  procedure  was  strongly  influenced 
by  problems  associated  with  the  experiment.   The  strong  reflection 
from  the  tube  ends  was  one  of  the  most  important  factors  in  all 
considerations.   The  need  for  gated  signals  rather  than  continuous 
waves,  in  order  to  study  the  effect  of  the  discharge,  meant  that 
the  data  could  not  be  taken  with  meters.   The  multi-component 
spectrum  associated  v;ith  the  bias  for  the  gated  sine  waves  meant 
that  the  signal  had  to  be  Fourier  analyzed. 

Due  to  the  strong  reflections,  the  signals  could  not  be 
analyzed  in  the  continuous  mode  with  a  Correlation  Function 
Computer  such  as  the  PAR  lOOA.   When  the  gated  mode  was  used,  the 
time  base  was  limited  to  the  period  from  the  triggering  signal  to 
the  start  of  the  reflected  signal.   Since  this  time  base  was 
sampled  only  100  times,  the  resolution  was  very  poor.   Preliminary 
analysis  of  the  signals  indicated  that  the  source  transfer  function 
decreased  at  low  frequencies.   Since  this  could  help  to  check 
theories  of  the  production  mechanism,  the  low  frequency  components 
were  of  particular  interest  and  resolution  was  needed. 

A  Boxcar  integrator  was  also  used  to  recover  the  signal  on  a 
scale  compjatible  with  transcription  into  numbers.   However  like 
the  Correlator,  this  utilized  only  one  signal  waveform  at  a  time 
and.it  wasted  most  of  that.   In  addition,  the  data  had  to  be  con- 
verted into  numbers,  and  card  punched  for  computer  processing. 
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Even    using    digital    voltmeters    to   obtain    the    data,    obtaining   more 
than    one   hundred   twelve  points    from  a   waveform  was    very    difficult. 

For    these   reasons,    the  process    of   recording    on    tape   and 
direct    conversion    to   digital    form  was    finally    selected.       It    offered 
the   greatest    resolution,    processing   possibilities    and   error    re- 
covery   capability. 

A.       mXA    REDACTION   PROCEDURE 

1 .      Ana].og    to    Digital    Conversion 

After    completing    the    data    recording,    the   data    tapes    were 
reviewed  with   an    osci] loscope.      One  minute   segments    were   cut    from 
each   of   the    original    tapes    and   spliced   together.       They   v;ere 
separated   by    long    sections    of   white   leader    tape   for    easy    identifi- 
cation   of    Lhe    segments. 

A  hybrid   computer    composed   of  a    CCMCOR    5000   analog    computer 
interfaced   vjith   an    SDS    9300    digital    computer    was    used   for    the    con- 
version   of   the   analog    data    to    digital    form.       The   analog    side    of    the 
hybrid   computer    was    set   up    to   receive   the   analog    signal   and  provide^ 
the  appropriate   signals    to   the    digital   portion   for    digitizing    the 
information.      The    gain    of    the    data    channels    was    50.      A   channel 
corresponds    to   the   waveform   recorded   on    it.      Logic   boards      were 
used   to   start    sampling      at   a    frequency      of    10   kHz    when    the   trig- 
ger   signal    reached      a    variable    trigger    level.         Two   switches    were 
used   to   start    and   stop    the    conversion   process.      The    number    of 
samples    taken   following    each    trigger    signal   v;as    controlled  by    the 
digital   program.      After    converting    1200    samples    of    each   of    the 
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three  data  channels  to  digital  form,  the  information  was  written 
on  a  magnetic  tape.   The  program  used  for  the  conversion  is  a 
standard  program  available  to  all  users  of  the  hybrid  computer. 
The  parameters  for  the  program  are  the  number  of  channels  to  be 
sampled  and  tlie  number  of  samples  per  channel  to  be  made  for  each 
record.   The  logic  involved  in  the  sampling  process  is  different 
for  most  conversions  and  this  required  a  slight  amount  of  patching 
on  the  logic  patch  boards. 

2  .   Seven  to  Nine  Track_  Conver  s  ion_s 

The  information  written  on  the  7  track  magnetic  tape  by 
the  hybrid  computer  is  in  the  form  of  2!^    bit  binary  words.   To  be 
processed  by  the  IBM  360,  these  must  be  converted  to  32  bit  binary 
words.   The  only  programming  involved  in  IIjIs  is  a  program  tc  read 
the  7  track  tape,  call  a  conversion  routine  and  write  the  informati 
on  another  tape,  usually  9  track  for  increased  efficiency,  for  later 
processing.   For  this  job,  this  step  was  also  used  to  separate  the 
samples  by  channel. 

3 .   Ensemble  Averaging  with  Data  Rejection 

The  diode  signals  contained  noise.   The  noise  was  random  and 
consisted  chiefly  of  positive-going  spikes.   A  low  pass  filter  could 
not  eliminate  them  without  distorting  a  significant  portion  of  the 
acoustic  waveform.   In  addition  there  were  a  few  bad  trigger  signals 
recorded.   Thus  it  ivas  known  from  the  beginning  that  some  of  the 
data  records  were  bad  and  should  not  be  included  in  the  processing 
but  it  was  not  known  which  they  were.   Since  the  entire  data  set 
consisted  of  ten  files,  each  of  which  contained  from  80  to  iCO 
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records  of  3^00  words,  screening  became  a  major  task.   The 
following  procedure  v/as  finally  found   to  be  most  efficient  and 
effective. 

Each  file  corresponded  to  an  experimental  run  and  all 
records  in  that  file  represented  the  same  waveforms  plus  noise. 
The  records  in  a  file  were  all  averaged  together  and  the  standard 
deviation  for  each  channel  was  obtained.   Following  this,  all 
records  were  compared  with  the  average.   Those  with  less  than  75% 
of  their  points  within  2.8  standard  deviations  of  the  average  were 
dropped  from  further  processing.   Using  only  the  acceptable  records, 
a  new  ensemble  average  and  standard  deviation  was  obtained  for 
each  channel.   These  were  used  in  a  final  ensemble  averaging  pro- 
cess in  which  only  those  Doinxs  within  2  stanuaru  uevlatioiiS  of 
the  average  were  used.   This  provided  an  ensemble  average  of  the 
waveform  and  the  standard  deviation  of  the  points  involved. 
^ •   Delay,  Overlay  and  Zero 

At  this  point,  the  data  consisted  of  1200  samples  for  each 
of  the  three  waveforms,  representing  their  average  values,  and  the 
standard  deviation  of  the  waveforms.   For  the  sampling  process 
1200  samples  were  taken  so  that  all  of  the  desired  waveform  of  the 
signal  at  point  2  was  sampled.   This  meant  an  undesired  signal  was 
included  in  the  samples  of  the  signal  at  1,  namely  the  reflected 
signal.   In  order  to  make  a  valid  comparison  of  spectral  components, 
the  time  samples  of  the  waveform  which  is  transformed  must  contain 
the  same  number  of  samples.   Reducing  the  number  of  samples  of  s   to 
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those  necessary  to  avoid  the  inflections  in  s   means  the  starting 
point  for  the  transformation  must  be  delayed.   Since  phase  infor- 
mation v.'as  not  going  to  be  considered  in  the  analysis,  and  time 
delay  merely  changes  the  phase  of  the  components,  delaying  the 
signals  would  not  affect  the  results.   Therefore  the  signals  s 
and  s   were  delayed  by  choosing  a  starting  point  for  each  such 
that  600  samples  included  the  waveform  up  to  the  point  where  the 
reflection  started. 

The  diodes  interacted  electrostatically  with  discharge. 
This  caused  a  small  coherent  square  wave  in  the  diode  output 
which  coincided  with  the  time  of  the  current  pulse.   Its'  inclusion 
in  the  waveform  to  be  transformed  would  give  erroneous  spectral 
compouents  for  the  acoustic  signal.   Delaying  v;ac  sufficient  to 
eliminate  this  unwanted  signal  from  s  .   In  order  to  eliminate  it 
from  s  ,  without  introducing  unwanted  spectral  components,  the 
values  of  the  samples  in  the  region  containing  this  unwanted 
signal  were  replaced  with  values  taken  at  random  from  the  adjacent 
region.   Since  this  adjacent  region  contained  only  samples  of  the 
waveform  before  arrival  of  the  acoustic  signal,  they  only  repre- 
sented background  noise. 

Since  the  initial  portions  of  the  signals  S,  and  S 
represent  no  signal,  i.e.  the  period  before  arrival  of  the  acoustic 
signal  from  the  sound  source,  its  average  value  should  be  zero. 
The  signals  were  offset  to  zero  by  first  determining  the  average 
value  of  the  no  signal  region  and  then  subtracting  it  from  all 
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values.   For  the  current  waveform,  this  region  followed  the 
current  pulse.   The  three  operations  of  delaying  the  signals, 
overlaying  the  electrostatic  signal  and  establishing  the  zero 
level  were  performed  in  one  program. 

B.   OPERATIONS  ON  THE  TRANSFORM 

Two  sequential  programs  were  used  to  transform  the  vjaveforms 
and  obtain  acoustic  informations  from  them.   The  first  program 
took  the  input  waveforms,  corrected  them  for  gain  and  diode 
sensitivity  and  determined  their  spectra.   The  second  used  the 
spectra  to  determine  the  transfer  functions  of  the  discharge 
region  and  tlie  source. 

Input  data  for  the  transform  program  were  the  three  ensemble- 
averaged  waveforms,  i(nT),  s  (nT),  s  (nT)  and  their  standard 
deviations,  Ai,  AS,,  AS„ ,  and  the  diode  sensitivities. 

Two  uncalibrated  amplifiers  were  used  in  amplifying  and 
recording  the  diode  signals.   These  were  calibrated  as  part  of  the 
processing.   While  recording  data,  both  amplifiers  were  fed  a 
common  signal  consisting  of  2  millivolt  peak  to  peak,  5000  Hz, 
tone  burst.   This  calibration  signal  was  processed  in  the  same 
manner  as  the  other  signals.   In  the  first  part  of  the  transform 
program,  the  spectrum  of  this  signal,  in  the  diode  channels,  was 
Fourier  analyzed.   Then  the  spectral  component  closest  to  5000  Hz 
(50^0)  was  used  to  calibrate  the  gain  for  each  diode  channels  so 
that  this  component  was  equal  to  one. 
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The  diode  sensitivities,  supplied  as  inputs,  were  calculated 
earlier  from  experimental  data  and  expressed  in  terms  of  milli- 
volts/torr.   The  sensitivity  of  the  current  detector  was 
5  X  10    volts/milliamp ,   The  sensitivities  and  corrected  gain 
factors  were  used  to  convert  the  values  into  the  correct  magni- 
tudes of  the  physical  quantities  involved,  milliamps  for  the 
current  and  torr  for  the  acoustic  signals. 

The  transforms  of  the  corrected  waveforms  were  then  calcu- 
lated using  a  computer  library  subroutine  "RHARM".   This  program 
requires  that  the  number  of  points  be  an  integer  power  of  two. 
There  were  6()0  points  available  and  the  nearest  power  of  two  was 
512,  leaving  88  excess  points.   This  means  that  there  is  some 
uncertaiijty  as  to  Vi'uere  to  start  the  arialysis.   For  this  reason, 
the  waveforms  were  analyzed  five  times,  each  time  advancing  the 
starting  position  by  10  places.   The  average  value  and  variance 
of  each  component  was  calculated  from  the  five  spectra. 

The  output  of  the  transform  program  was  the  average  spectrum 
and  the  variance  of  its  components.   Only  the  first  'twenty  five 
of  these  were  passed  to  the  next  program.   The  magnitudes  of  the 
remaining  components  were  less  than  5%   of  the  fundamental  and  were 
characterized  by  a  great  deal  of  scatter. 

The  next  program  used  this  data,  and  calculated  the  acoustic  and 
source  transfer  functions.   The  equations  of  Chapter  V  were  used 
for  the  calculations.   A  computer  library  subroutine  "LSQP]^2"  was 
used  to  obtain  the  least  squares  determination  of  p  .   For  this 
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calculation,  the  data  was  v;eighted  inversely  proportionally  to 
the  Square  of  its  uncertainty,  with  weights  normalized  such  that 
their  sum  was  equal  to  the  number  of  points. 

For  the  data  corresponding  to  no  discharge,  Eq.  (5.10)  was 
used  to  obtain  an  estimate  for  the  sound  spectra  at  the  source. 
This  was  divided  by  the  sound  source  current  spectra  to  obtain 
the  transfer  function  of  the  sound  source. 
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VII.       PRESENTATION    AND    DISCUSSION    OF    RESULTS 

A.  AVERAGE  WAVEFORMS 

The  averaging  and  screening  method  worked  very  well  although 
the  program  could  be  rewritten  to  perform  the  job  more  efficiently. 
The  three  waveforms  for  the  discharge-off  case  at  7-^5   torr  are 
shown  as  Figures  10,  11  and  12  just  before  being  transformed.   For 
contrast,  the  first  2/3  of  the  average  \;aveform  for  s   for  the 
same  run  is  shown  as  Figure  13  with  the  same  scale.   The  second 
pulse  is  the  reflected  signal  from  the  end  close  to  the  sound 
source.   This  shows  the  effect  of  the  delaying,  overlay  and  zeroing 
operations.   Unless  otherwise  specified,  the  remainder  of  the 
figures  displayed  in  this  chapter  were  derived  from  these  wave- 
forms . 

B.  ACOUSTIC  TRANSFER  FUNCTION 

The  points  of  the  acoustic  transfer  function  showed  a  rea- 
sonable amount  of  scatter  about  the  predicted  form.   The  greatest 
degree  of  scatter  was  found  to  be  associated  with  those  components 
whose  contributions  to  the  signal  were  small.   Figure  l4  shows 
the  first  twenty  five  spectral  components  of  the  sound  at  point  2. 
Figure  15  shows  the  transfer  function  in  a  linear  form.   Figure  l6 
is  the  same  information  expressed  as  attenuation  of  sound  in 
decibels  per  meter  as  a  function  of  frequency.   The  curve  repre- 
sents the  least  squares  fit  to  the  data  under  the  assumption  that 
boundary  losses  dominate.   The  goodness  of  fit  found  in  all  cases 
supports  the  validity  of  this  assumption. 
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FigutG  10.  Sound  Source  Discharge  Current.  P  =  7-65  torr 
Main  Discharge  off.  Time  Scale:  2,048  msec/in 
Vertical    scale    is    arbitrary. 
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Figure    11.       Sound  Signal    at   Point    1.    (delayed,    overlayed 
and   zeroed).      P    =    7.65   torr.       Main   discharge 
off.       Time   scale:    2.0^8   msec/in.       Vertical 
scale    is    arbitrary. 
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Figure    12.       Sound    Signal    at    Point    2.     (dolayed,    overlayed 
and   zeroed).      P    =    7.65   torr .      Main    discharge 
off.      Time    scale:    2.048   msec/in.    Vortical 
scale    is    arbitrary. 


Figure    13.    Sound   Signal   at    Point    1.    (Uncorrected). 

P    =    7.65    torr.      Main    discharge    off.      Time 
scale    =    2.048   msec/in.    Vertical    scale    is 
the   same   as    P'igure    11. 
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Sj.ectral  Decomposition  of  Sound  Signal  at  Poirt  2 
This  is  the  spectral  decomposition  of    " 'ivr  ■ 
form  shown  in  Figure  1.-.   Dashed  line  is  a 
simple  connection  of  points. 

106 


Figure  15.  S  (v)/S  (v)  in  Dopendence  on  V ^  for  a  Pressure 
of  7.65  torr .   Main  discharge  off.   Error  bars 
represent  one  standard  deviation.   Curve 
represents  the  weighted  least  squares  fit  to 
the  data. 
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Figure  l6.   Attenuation  in  Dcp)cndence  on  Frequency  for  a 
pressure  of  7-65  torr .   Main  discharge  off. 
Curve  represents  the  least  squares  fit  to  the 
data.   This  figure  corresponds  to  Figure  15. 
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C.       SOURCE   TRANSFER    FUNCTION 

The   source   transfer    function   was    found    to  have    the    low   fre- 
quency   drop    off    indicated  by    the   preliminary    observations.      At 
:ffequencies    above   a   pressure-dependent    transition    frequency,    the 
frequency    characteristic    did   follow   the  predicted   form.      The 
magnitude   was    also   found   to  agree   with    the   predicted   values.      The 
transfer    functions    are    displayed  as    Figures    17,    l8   and    19   for    the 
filling  pressures    7.65   torr ,    12.7   torr    and    l6    torr    respectively. 
The   structure    is    believed   to   be    due    to   scatter    in    the    data    since 
the  points    with    the   greatest    departure   from   the    smooth    curve 
correspond   to    those   spectral    components    with   the    largest    relative 
error  . 


Two   aspects    of    the   results    can   be   used   to   test    the    theory 
of   sound  production    from   a    discharge.      One    is    the    frequency    de- 
pendence  and    the    other    is    the   magnitude. 

a.      Frequency    Dependence 

The    low-frequency    drop    off    is    contrary    to    the    fre- 
quency   dependence  predicted   by    Schulz's    results.      This    effect    is 
believed   to   be   due    to   neutral   particle    diffusion    under    the    in- 
fluence  of   thermal    gradients.       It    is    a    low-frequency,    low-pressure 
effect    which    should   appear    only   when    thermal    energy    sources   provide 
the   acoustic    energy.      This    will    be    discussed   later. 

According    to   Schulz's   prediction,    the    transfer    function 

should   vary   as 

.       ,  0^  . 
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Figure  1?.   Source  Transfer  Function  in  Dependence  on 

Frequency  for  p  =  7-65  torr ,   Main  discharge 
off.   Vertical  scale  in  units  of  (newton/m-) 
per  ampere.   Solid  curve  represents  predicted 
form  fitted  at  one  point.. 
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Figure   l8. 


Source  Transfer    Function    in    Dependence    on 
Frequency    for    p    =    12./    torr .      Main    discharge 
off.      Vertical    scale    in    units    of    ( newtons/ra" ) 
per    ampere.       Solid   curve    represents    predicted 
form    fitted   at    one   point. 
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Figure  19.   Source  Transfer  Function  in  Dependence  on 

Frequency  for  p  =  l6  torr .   Main  discharge  off, 
Vertical  scale  in  units  of  (newton/m^)  per 
ampere.   Solid  curve  represents  predicted  form 
fitted  at  one  point. 
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This  prediction  was  tested  by  fitting  a  curve  given  by 

through  the  transfer  function  data  points.   The  coefficient  C 
was  evaluated  by  causing  the  curve  to  pass  through  the  point  at 
V  =  2450  Hz . 

Schulz's  derivation  was  based  on  planar  electrodes 
separated  a  distance  "a".   The  experiment  used  cylindrical  elec- 
trodes of  length  1^  cm  separated  by  a  distance  of  3  cm.   Therefore 
the  effective  length  of  the  discharge  could  lie  anywhere  between 
3  and  6  cm.   The  value  of  "a"  was  adjusted  within  these  limits  to 
obtain  a  curve  which  fitted  the  data  reasonably  well.   For  7.65 
torr .  the  value  of  a  was  3 • 5  cm;  for  16  torr ,  the  value  was 
4.6  cm.   Since  the  pressure  affects  the  discharge,  there  was  a 
possibility  that  the  value  of  "a"  for  12.7  torr  could  be  inter- 
polated from  the  other  two  values.   The  value  obtained  by  inter- 
polation was  3.93.   This  gave  a  good  fit  to  the  data.   The  smooth 
curves  shown  in  the  figures  are  the  fits  of  equation  (7.1)  to  the 
data  . 

b.   Magnitude 

Using   Schulz's    equation,    the   source   transfer    function 
is   predicted  as  sin(-r— ) 

The   energy    transfer    rate   from   the    electrons    H    is    the 
same  as   used   earlier    in    conjunction   with  wave   amplification.       It 
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may  be  expressed  in  terms  of  the  characteristic  time  for  electron 

heating  as  2 

^  c 

H  =  - 


T(y-l)  • 

2  , 
The  calculated  value  of  c  /T  v;as  obtained  with  the 

computer  program  used  for  evaluating  the  Ingar d-Schulz  equations. 

This  was  accomplished  for  the  transfer  function  at  v  =  2540  Hz 

and  p  =  16  torr .   The  discharge  current  spectral  component  at 

this  frequency  was  6.4  ma  which  was  estimated  to  correspond  to 

10-3 
an  electron  density  of  2  x  10    cm   .   In  determining  this  density, 

the  cross  sectional  area  of  the  electrodes  was  used  since  the  dis- 
charge did  not  extend,  in  the  radial  direction,  significantly  be- 
yond them.  The  procedure  used  was  the  same  as  outlined  in  Sec.  V.E. 

The  predicted  value  of  Z^  obtained  was 

S  ,    2 

S  ampere  ^  "^ 

The   corresponding    value,    experimentally    obtained,    was 

Z.   -   33.2   "^/"'' 


'S  ampere 

The  agreement  of  the  two  values  is  good.   Considering 
the  number  of  approximations  made,  order  of  magnitude  agreement 
would  have  been  considered  satisfactory.   No  error  analysis  was 
conducted  on  this  computation  for  that  reason. 

Except  for  the  departure  at  low  frequencies,  the  re- 
sults are  consistent  with  the  theoretical  predictions  based  on  the 
assumption  of  energy  transfer  by  the  electrons  in  the  positive 
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column.      There   are   aspects    of    the   problem  which   were  not    examined 
by   this    experiment    but   which    should   be   studied   to   gain   complete 
understanding    and   to   understand   other    e>rper  imental    results.      A 
more    complete    investigation   would    include  phase   as    well   as    magni- 
tude   information    in    the   transfer    function.       In   addition,    the 
source    of   the    sound   could   be    localized  and  momentum   transfer 
aspects    investigated.      Although    this    experiment    indicates    that 
electron   heating    of   the   neutrals    in    the  positive    column  provides 
a   mechanism   for    sound  production    in    glow    discharges,    it    has    not 
discounted  the   roll    of    ions    in    sound  production    under    other 
cir  cums tances . 

2 .      Low   Frequency    Departure 

During   preiiiuinaxy    vvuik    in    coijj'unct  ion    with    this    e:cpGr  i- 
ment ,    a    large    number    of  acoustic   waveforms    were    observed.       One 
characteristic    of   the   sound   signal   produced   by   a    pulse   v;as    noticed, 
which    is    reflected    in    the    drop-off   of    the   transfer    function   at 
low   frequencies.       Figure   20    shows   part    of    the   waveform   received  at 
point    1    during    one    of   the   runs.       EXie   to   attenuation    of   the   high 
frequency    components,    the   features    are   not   as    sharp   as    some   which 
were   observed  at    closer    positions    during    earlier    work.      The   main 
features    are   visible   however. 

The   first   and  most    significant    feature    is    that    the   vertical 
distance   from    1    to   2    is    approximately    equal    to    the   distance   from 
3    to   4,    with    the    latter    always    slightly   greater.       The    next    feature 
is    that    the   slope  and    curvature   of    the    curve   from   2    to   3    is   approxi- 
mately   the   negative   of   that    for    the    curve   from   4    to    5.       The   two 
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Figure   20.       A  Typical    Sound    Pulse   Produced   by   a 
Rectangular    Current   Pulse. 
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curves  suggest  an  exponential  process  with  a  long  time  constant. 
A  last  point  is  that  point  5>  where  the  slope  is  approximately 
zero,  is  lower  than  point  1,  the  level  before  the  signal.   The 
polarity  of  the  figure  is  such  that  point  1  corresponds  to  the 
ambient  pressure  with  upward  deflections  indicating  an  increase 
in  pressure  in  the  discharge  region. 

The  following  description  of  events  occurring  in  the  dis- 
charge region  during  the  creation  of  the  sound  pulse  is  postulated. 
When  the  discharge  current  increases,  the  energy  transferred  to 
the  neutral  gas  increases,  thereby  increasing  the  temperature  of 
the  neutral  particles.   This  temperature  increase  is  so  rapid  that 
the  netural  particle  density  remains  essentially  constant.   As  a 
consequence,  the  prccsurc  changes  almost  direc+ly  with  thp.  tem- 
perature.  The  neutral  particles  move  away  from  the  discharge, 
radially  as  well  as  axially,  under  the  temperature  and  pressure 
gradients.   As  the  concentration  decreases,  the  pressure  decreases 
toward  the  ambient  value.   The  pressure  change,  Ap, ,  which  would 
take  place  because  of  the  diffusion  is 

Ap   =  KANT 

where  AN  is  the  change  in  density. 

When  the  heating  processes  are  removed,  the  neutral  gas 
cools  rapidly  to  the  ambient  temperature.   Since  the  cooling 
process  is  much  more  rapid  than  the  particle  diffusion,  the  pres- 
sure at  ambient  temperature  would  be  the  equilibrium  pressure 
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minus  the  pressure  difference  caused  by  the  density  change.   This 
pressure  difference,  Ap„,  is  given  by 

Ap^  ^  KANT^. 
Since  T  "^  T  ,  Ap,  "^   Ap„.   The  neutral  particles  would  then 
diffuse  back  in,  to  return  the  region  to  the  ambient  density. 

The  diffusion  process  is  slow,  with  a  long  time  constant. 
Under  the  action  of  a  high  frequency  signal,  there  would  not  be 
time  for  a  significant  amount  to  occur  between  cycles.   Consequently, 
the  acoustic  pressure  would  not  reflect  any  density  changes.   At 
the  other  extreme,  for  very  low  frequency  temperature  variations, 
the  pressure  changes  would  be  negligible  but  the  density  changes 
would  be  large.   Therefore  it  can  be  predicted  that  the  effect 
will  become  appreciable  at  loiv  frequencies. 

Notice  that  in  the  postulated  process,  no  assumptions  were 
made  concerning  the  heating  process  except  that  it  is  spatially 
uniform  in  one  region  and  that  particles  are  free  to  move  in  and 
out  of  the  region.   Thus  the  assumption  of  one  dimensionality  for 
such  a  process  occurring  in  a  tube  could  be  a  valid  assumption. 
Therefore  Schulz's  assumptions  will  be  considered  valid  and  the 
failure  of  the  prediction  at  low  frequency  attributed  to  the 
wave  equation.   The  wave  equation  fails  because  it  does  not  take 
diffusion  into  consideration. 

If  diffusion  is  the  important  factor  in  the  discrepancy  at 
low  frequencies,  then  the  effects  will  become  important  for  fre- 
quencies whose  period  is  of  the  order  of  or  longer  than  a  character- 
istic diffusion  time  T  .   This  time  is  the  time  for  particle 
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diffusion  for  a  mean  displacement  of  some  characteristic  dimension 
of  the  discharge  region.   Therefore,  the  characteristic  frequency 
is  given  by  McDaniel  [39]  as 

=  i  =  -2. 

^C  -  T  -  ^2 

where  C  is  the  characteristic  dimension.   The  effective  diffusion 

coefficient  D  is  given  by  Present  [36]  as  D  =  -^  where  v  is  the 

mean  velocity  and  H    is  the  mean  free  path. 

Since  the  mean  free  path  varies  inversely  with  pressure, 

this  predicts  that  the  frequency  at  which  the  effect  occurs  will 

vary  inversely  with  pressure.   An  examination  of  Figures  15,  I6 

and  17  shows  that  this  is  the  case:  at  ].6  torr  ,  v   ^  500  Hz;  at 

12.7  torr,  V  ="  750  and  at  7.^5  toxr  ,  v  ~    ICOC  \\z ,    the  product 
c  c 

pv   being  roughly  constant  as  it  would  for  the  predicted  dependence 
on  density.   However,  substitution  of  the  discharge  length  "a"  for 
the  characteristic  dimension  gives  a  frequency  which  is  three  or dei s 
of  magnitude  too  low.   Thus,  although  diffusion  effects  may  be 
responsible  for  the  low  frequency  effects,  the  simple  argument 
used  to  determine  the  critical  frequency  did  not  predict  the 
correct  result. 

This  diffusion  effect  is  a  new  problem  for  acoustics.   It 
is  a  low  frequency,  low  pressure  effect  and  will  probably  only 
occur  in  conjuction  with  a  thermal  acoustic  source.   The  excitation 
of  sound  by  a  heat  or  energy  source  is  self -rectify ing  since  the 
energy  transfer  into  the  system  cannot  be  balanced  by  the  removal 
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of  heat  or  energy  from  the  system  by  the  same  process.   This 
means  that  the  energy  source  term  H  in  equation  (2.1)  can  not  be 
represented  by  the  form 

H(r,t)  =  H(r)e^'^ 

since  such  a   representation  implies  that  energy  can  be  removed  by 
the  source  with  the  same  efficiency  that  it  transfers  energy  into 
the  source. 

To  predict  the  correct  low  frequency  characteristics  of 
sound  produced  by  a  discharge  plasma,  the  derivation  of  the  wave 
equation  must  take  diffusion  and  heating  effects  into  consideration 
The  results  of  the  derivation  should  be  a  differential  equation 
capable  of  predicting  the  neutral  density  and  temperature  when 
solved  for  the  zero  frequency  case.   Such  a  derivation  is  con- 
sidered beyond  the  scope  of  this  investigation. 

D.   EFFECT  OF  THE  DISCHARGE 

1 .   Predicted  Attenuation 

When  the  measured  values  of  axial  temperature  were  com- 
pared with  the  predictions  of  Venzke  et.al.  L5l]j  the  agreement 
was  within  1%.      This  verification  of  the  predictions  of  that  paper 
was  considered  sufficient  justification  to  use  its  predictions  for 
T  ,  the  wall  temperature  and  T.,  the  mean  temperature  with  the 
same  relative  error.   These  values  were  substituted  into  the  ex- 
pressions, of  Section  V.E.  to  obtain  the  neutral  density  in  the  two 
regions  for  each  run.   The  steps  outlined  in  that  section  were 

followed  to  obtain  estimates  for  N   and  T  .   Vvith  the  density, 

e       e 
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wall  temperature,  electron  temperature  and  electron  density  as 
inputs,  the  attenuation  constant  was  calculated  using  the  Ingard- 
Schulz  equations  for  CO  =  31,400  and  CO  =  600. 

The  results  obtained  were  that  the  zero  ionization  value 
of  the  attenuation  constant  was  essentially  the  same  as  the  value 
predicted  for  a  discharge  under  the  experimental  conditions.   The 
value  of  p   calculated  in  the  manner  mentioned  earlier  was  larger 
at  the  higher  frequency,  which  indicated  that  bulk  losses  were  not 
completely  negligible.   The  largest  difference  between  the  values 
of  p   for  the  two  frequencies  occurred  at  7.65  torr  where  the 
difference  was  slightly  less  than  3%  for    the  un-ionzed  value.   The 
greatest  ionization  effect  predicted  was  a  .5%   reduction  in 
attenuation  for  the  two  discharge  runs  made  with  a  filling  proc;- 
sure  of  7.65  torr.   In  other  words,  the  Ingar d-Schulz  equations 
predicted  that  no  plasma  effects  would  be  discernable  under  the 
experimental  conditions. 

For  the  plasma  associated  with  a  filling  pressure  of 
7.65  torr,  the  electron  temperature  T   was  estimated  to  be 

1.4  X  10  K  and  the  electron  density  was  5  x  10    cm  ^  .   With  these 

17    -3 
values  and  the  values  for  the  neutrals,  N   =  2.4  x  10    cm    and 

n 

T   =  323K,  the  values  of  some  of  the  plasma  parameters  were  cal- 

9 
culated.   The  characteristic  collision  frequencies  are: -CO   =  2.9  xlO  , 

^  '       en 

7  -1 

and  60.       =   fjO        =    5.87   X    10      sec       .      The   electron   and    ion  plasma 
xn  nn 

frequencies    (CO   ,    60.  )  »    are  4.00   x    10        and    1 . 5  x    10      cm      .      The  mean 
-1  ^    e         1 
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-4 

free  path  of  a  neutral  particle  is  7.0  >:  10    cm.   While  these 
values  are  different  for  the  other  discharges,  the  values  are 
representative  of  the  orders  of  magnitude  of  the  various  para- 
meters . 

2 .   Comparison  of  Predicted  and  Experimental  Values 

The  aissumption  of  boundary  loss  domination  was  made 
earlier.   Since  the  failure  of  this  assumption  would  occur  at 
higher  frequencies,  where  the  data  were  weighted  the  least  due 
to  their  larger  r elat ive  err  or ,  the  fit  to  the  data  should  reflect 
the  low  frequency  characteristics.   This  was  tested  at  the  value 
for  p   predicted  for  CO  =  600.   The  predicted  values  are  shown  in 
Table  III  under  the  column  marked  p'.   The  experimentally  deter- 

inincd  V-lu'^  of  ^      a^d  its  c  +  pnriArH  r!<^\Ha  + -ion  p^Tf^    shown  in  thp> 
•  o 

columns  marked  B   and  AB  . 
•^  o       ^  o 

The  ratio  of  the  difference,  between  the  predicted  and 

experimental  value  of  p  ,  to  the  standard  deviation  of  p   is 

°     p  -p' 
^  o  ^  o 

shown  m  the  column  headed  ~ .   The  magnitude  of  these 

AB 
'  o 

numbers  are  all  less  than  unity,  indicating  that  the  differences 
between  the  predicted  and  the  experimental  value  are  all  within 
one  standard  deviation  of  the  result.   The  greatest  difference 
between  the  predicted  and  experimental  values  is  13%. 

However  the  signs  of  the  difference  are  correlated.   The 
second  column,  headed  "polarity",  indicates  the  polarity  of  the 
discharge  with  respect  to  sound  propagation.   A  plus  sign  means 
sound  propagation  in  the  direction  of  electron  drift  motion.   A 
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zero    is    used   to    indicate   the   discharge-off   state   and  a   minus    sign 
indicates    sound  propagation    counter    to    electron    drift.       The    sign 
of   the    differences    is   positive   for    the    no-discharge   cases.       For 
the    discharge   cases,    it    is    negative.       The   negative    sign    indicates 
that    the  predicted  attenuation    is    greater    than    the    observed 
attenuation . 

3 . "   Discussion  of  the  Results  of  the  Comparison 

.  The  results  show  that  for  the  ionized  cases,  the  predicted 
value  of  the  attenuation  constant  is  greater  than  the  observed 
value.   Although  the  two  values  lie  within  one  standard  deviation 
of  the  experimental  value,  the  consistency  of  the  result  indicates 
that  it  is  due  to  either  a  systematic  error  or  to  a  plasma  effect. 
The  temperature  dependence  of  the  attenu? +  3  op  con^t-=>nt  at  rnnstant 
pressure  is  such  that  if  the  v/all  temperature  used  for  the  com- 
putation were  too  low,  the  predicted  values  of  p   would  be  too 
high.   The  three  parameters  which  influence  the  predicted  values 

to  any  degree  are  CO,  N   and  T  .   The  error  must  be  due  to  either 

n  n 

N      or    T      since   although    the    value    of   CO   is    the   same    for    the    ionized 
n  n 

and  un- ionized   cases,    the   results    of   the   comparison    show   that    the 
sense   of   the   errors    involved   is    different. 

In   calculating    the  predicted   values,    the  propagation    of 
error    was    followed.      Although    the   mean    temperature    in    the    dis- 
charge  was    assumed   to   be  uncertain   by    10%,    it    did    not   affect    the 
error    of   the   corrected  pressure   significantly.      That    error    was 
affected  primarily   by    the   uncertainty    in    the    value    of   the    filling 
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pressure,  which  was  estimated  to  be  5%.   If  the  error  in  filling 
pressure  was  systematic,  it  would  have  affected  all  cases  in  the 
same  manner.   Therefore  neither  an  error  in  the  mean  temperature 
nor  an  error  in  the  filling  pressure  is  believed  to  be  a  source 
of  systematic  error. 

The  only  remaining  source  of  systematic  error  would  be  the 
value  of  wall  temperature  used.   This  affects  the  results  through 
the  density,  calculated  using  T   and  p',  and  directly.   The  error 
in  the  wall  temperature  was  conservatively  over-estimated  to  be 
10  ,  which  means  a  25%  error  with  respect  to  the  increase  above 
ambient  temperature  but  only  a  3%  error  with  respect  to  the  ab- 
isolute  temperature.   Representing  the  error  in  p'  due  to  the 
error  in   T   by  up  ,  it  is  given  by 


^T   I  bN    bT  OT/ 


'^^w 


Since 

N   =  -^ 
1         KT 

for  constant  pressure, 

5n  _  -N 
bT  T  ' 

which  gives 

^T   \^  T  ON     bT  J       CO 

Using  this  expression,  6p   was  calculated  for  AT   =  10.   The 
results  are  displayed  in  Table  III.   It  is  apparent  that  these 
values  are  not  large  enough  to  account  for  the  discrepancy  betv;een 
results  except  for  runs  7  s^nd  9.   It  should  be  noted  that  if  the 
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value  of  the  filling  pressure  used  for  predicting  p      for  runs  7, 
8  and  9,  were  high,  then  the  correct  value  for  the  prediction  would 
be  higher  for  all  three  cases.   This  would  decrease  the  dis- 
crepancy for  the  no  discharge  case  and  cause  the  results  for  the 
ionized  cases  to  be  more  like  the  others.   The  conclusions  to 
be  drawn  from  this  discussion  are  that  while  a  systematic  error  in 
the  value  of  the  wall  temperature  used  for  predicting  the  results 
is  possible,  it  is  not  likely  to  be  the  cause  of  the  difference 
between  the  experimental  and  predicted  results. 

The  two  runs  where  this  error  should  be  the  smallest  are 
runs  2  and  3}  since  the  temperature  increase  is  less  at  lower 
pressures.   The  difference  for  these  two  runs  is  not  significantly 
lower.   If  the  resuixs  for  filliug  pxessurtes  7-65  and  l6  tcrr  r,lcnc 
are  compared,  then  it  appears  that  the  polarity  docs  have  an  effect. 
However,  there  is  insufficient  replication  to  support  such  a  con- 
clusion. 

The  conclusion  supported  by  the  results  is  that  for  these 
experimental  conditions,  amplification  of  the  sound"  may  have 
occurred  and  a  polarity  effect  may  be  present.   The  observed  effects 
are  within  the  experimental  error. 

E.   EVALUATION  OF  THE  YATSUI  EXPERIMENT 

Yatsui,  et.al.  [8,9]  conducted  their  experiment  in  Argon  at 
low  pressures  and  high  frequencies.   For  a  comparison  of  their 
experimental  results  with  the  Ingar d-Schulz  prediction,  the 
Ingar d-Schulz  equations  have  now  been  evaluated  v;ith  paramoter 
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values  corresponding  to  two  points  in  the  Yatsui  experiment.   For 

V  =  10   Hz  and  1  =  5  rnA,  they  reported  maximum  amplification 

-2 
occurred  at  1.5  x  10    torr  and  damping  occurred  for  pressures 

greater  than  .31  torr.   For  the  condition  of  maximum  amplification, 

the  ratio  of  the  received  signals  with  and  without  a  discharge  was 

3.3.   The  predicted  value  for  this  set  of  experimental  parameters 

will  be  discussed  first. 

Since  the  signal  traveled  only  a  fraction  e  of  its   total 

travel  in  the  discharge,  the  ratio  of  the  two  signals  is 


y-  -(e{3p+(l-e)p)d    -£(Pp-p)d 


S-p  _  e 


Vg  ^-Pd 


-  e 


where  d  is  the  total  distance  between  source  and  receiver.   There- 
fore, the  difference,  p  -p,  between  the  two  attenuation  coefficients 
can  be  determined  from  their  data  by 

-  1    ^^ 
Pp  ^    ed    ,   S 

The  denominator,  ed,  is  the  distance  traveled  in  the  discharge, 
38  cm.   Therefore  the  experimentally  determined  value  is 

p  -p  =  ^^   3.3  =  -  3.12  x  10"^  cm"-"-. 

Using  the  procedures  outlined  in  Section  V.E.,  the  electron 

4 
temperature  was  estimated  to  be  5-5   x  10  K  and  the  electron 

10    -3 
density  was  estimated  to  be  5.1  x  10    cm   .   With  these  values, 

and  CO  =  6.28  X  lO'',  N   =  5  x  10    cm""^,  and  T   =  300K,  the  quantity 

n  n  ^ 
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p  -p  was  computed  using  the  computer  program.   This  calculation 

-4 
gave  p  -p  =  -  5.3  x  10   ,  which  is  lower  than  the  experimental 

value  by  two  orders  of  magnitude.   Using  liberal  estimates  of  the 

errors  in  the  input  variables,  the  uncertainty  of  this  result  is 

-3 
A( P  ~P )  ~  1 • 52  X  10   .   This  is  not  enough  to  give  order  of 

magnitude  agreement. 

The  Ingar d-Schulz  predictions  were  also  evaluated  for  the 

discharge  conditions  corresponding  to  the  critical  pressure, 

17    -3  10 

.31  torr,  using:  N   =  1.1  x  10  '  cm   ,  N   =  5.5   x  10   , 

n        .  e 

k 
T   =  1.6  X  10  K  and  T   =  3OOK .   The  value  predicted  was  found  to 
en  ^ 

^®  P  -f  =  -  9.54  X  10"^°  cm"^  . 

P 

Since  this  difference  would  hardly  be  detectable,  it  is  in 
agreement  with  the  Yatsui  results.   Further  agreement  with  the 
results  is  the  prediction  that, for  pressures ' above  this,  damping 
occurs . 

Thus,  although  the  Ingar d-Schulz  equations  do  not  predict  a 
directional  effect,  its  predictions  are,  to  a  certain  degree,  in 
agreement  with  experiment. 


128 


VIII.   SUMMARY  AND  CONCLUSIONS 

This  study  has  examined  the  creation  and  modification  of 
sound  by  a  discharge  plasma  both  theoretically  and  experimentally. 
Following  a  review  of  earlier  experiments  on  the  creation  of  sound, 
a  mechanism  for  sound  production  by  ions  was  proposed  and  shown 
to  be  feasible.   Then  a  numerical  analysis  of  the  equations  of 
Ingard  and  Schulz  was  performed.   This  analysis  showed  that  net 
amplification  would  not  occur.   It  also  showed  that  the  approxi- 
mate form  of  the  dispersion  relation  was  not  useful  because  the 
criterion  for  its  applicability  was  invalid. 

Incidental  to  the  experimental  study,  the  theory  of  the 
mobility-limited  thermionic  diode  used  as  a  microphone  was  de- 
veloped beyond  lib  previous,  eleiaentary  form.   This  analysis  and 
the  subsequent  successful  usage  as  an  acoustically  transparent 
microphone  for  two  point  measurements  has  increased  its  usefulness 
as  an  experimental  tool. 

The  diodes  were  used  in  an  experimental  investigation  of  the 
modification  and  creation  of  sound  at  medium  pressures.   Fourier 
transforms  and  computer  data  processing  techniques  were  used  in 
the  analysis  of  the  data.   The  investigation  into  the  modification 
effect  of  the  discharge  was  in  a  different  region  of  current, 
pressure  and  frequency  from  earlier  experiments.   The  results  in- 
dicate that  an  amplification  mechanism  may  be  operative  and  that 
there  may  be  directional  effects.   Since  the  Ingar d-Schulz 
equations,  evaluated  for  the  experimental  conditions,  predicted 
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no  effect  would  be  detectable,  these  results  do  not  support  the 
equations.   The  Ingar d-Schulz  equations  were  also  evaluated  for 
the  conditions  of  the  Yatsui  experiment.   The  predicted  values 
were  not  found  to  be  in  agreement  with  the  experimental  results 
but  the  trend  with  pressure  was  in  agreement.   The  conclusion 
drawn  was  that  the  theory  needed  more  development. 

The  major  contribution  of  this  study  is  believed  to  be  the 
information  gained  from  what  was  originally  the  secondary  experi- 
ment, the  conversion  of  discharge  current  fluctuations  into  sound. 
This  is  the  first  known  investigation  in  which  the  current  magni- 
tude was  compared  to  the  sound  it  produced.   The  results  were  in 
reasonable  agreement  with  the  Ingar d  formulation  of  the  problem. 
A  low  frequency  effect  was  observed  which  was  interpreted  as  the 
first  observation  of  the  effects  of  particle  diffusion  in  acoustics 

The  general  agreement  of  Ingards  predictions  in  both  areas 
indicates  that  electron  heating  of  the  neutrals  is  the  dominant 
source  of  acoustic  energy.   However  several  experiments  have 
raised  questions  not  answered  with  the  present  theory.   One  of 
these  is  the  directional  effect  on  discharge  modification.   Two 
others,  associated  with  sound  production,  are  the  effects  of  high 
fields  and  the  localization  of  the  sound  source  near  the  cathode. 
These  stand  out  as  subjects  of  further  experimental  and  theoretical 
development.   The  interpretation  of  the  low  frequency  behavior  of 
the  sound  source  as  a  diffusion  effect  also  warrants  further  work, 
both  experimental  and  theoretical.   In  addition,  the  study  of  the 
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Ingar d-Schul2  equations  indicates  that  the  increased  attenuation 
to  be  expected  under  certain  conditions  could  be  a  fruitful  area 
for    investigation. 
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APPENDIX  A 

THE    MOBILITY-LIMITED  THERMIONIC    DIODE   AS   A   MICROPHONE 
Standard  microphones    have    limitations    which   hinder    their 
usage  when   studying   sound  at   very    low  pressures.      This    is   parti- 
cularly   true    in   conjunction   with   plasma-acoustic   work.      There  are 
several    reasons    for    this.      A   very    important    one    is    contamination 
from  gassy    components.      Another    consideration    is    that    the   micro- 
phones  were    designed  for    operation   at   atmospheric  pressure. 
Operation    at    reduced  pressures    may    require   circuit   modification 
since,    for    condenser    microphones,    breakdown    can    occur.       Chalupnik 
et.al.    r52J    reported   changes    in    sensitivity    and   frequency    response 
with  pressure    in    two    commonly    used   condenser    microphones    (Western 
Electric    640AA,    Bruel    and  Kjaer    type   4111).      The   lack    of   a    sensi- 
tive  sound    detector    without    the  problems    outlined  above   has 
slowed  plasma-acoustic    studies.       Recently    Dayton    et.al.    L 53l    re- 
ported  the    use    of  a   mobility-limited   diode  as    a   microphone   without 
the    drawbacks    of   conventional    microphones   at    low  pressures.      How- 
ever   heretofore   the   theory   has    not    been   developed  beyond    its   most 
elementary    form. 
1.       Theory 

If   the    cathode    is    sufficently    hot,    so   that   the    operation 
is    not    emission    limited,    the   current    flow    in   a    thermionic    diode, 
operating    in    a    high   vacuum,    is    limited  by    the   effects    of   space 
charge.      The    current    density    is    given    by   the   Langmuir-Childs    law 
which    says    that    the   density    is   proportional    to    the   cathode-plate 
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potential    to   the   three-halves   power.      The   constant    of  pro- 
portionality   is    dependent    on    the   geomentry    of   the   current-flow 
situation.       The   three   halves   power    dependence  assumes    an    equi- 
potential    filament;    a    condition   which    is    approximately    true  pro- 
vided  the   voltage    drop    across    the   filament    is    small    compared   to 
the   cathode-anode  potential.       If   this    condition    is    not   met,    then 
the  current    density    varies    with    the  plate   potential    to   a    higher 
power,    up    to   5/2    C54'J,    the  actual    value    of   which    is    determined  by 
the   definition    of   the    cathode   to   anode  potential, 
a.      Mobility-limited   Operation 

In   the  presence    of   a    background  gas    of   neutral   parti- 
cles,   the    electrons    will    suffer    collisions    and   their    average 
velocity   v.'ill    therefore   deer  ease  c      This    affects    the    space   charge 
field  and    the    current    density.       In    the    limiting    case,    where   the 
anode-cathode   separation    is    much    larger    than   a   mean    free  path,    the 
velocity    of   the   electrons    is    limited   by    their    mobility    in    the   gas. 
Cobine   L55J    derives    an    expression   for    the    current    density    in   a 

parallel-plate   diode   with   a   mobility-limited   current: 

-,  2 

j    =    9.92    X    lO"-^^      til—      (amp-cm'    )  (A.l) 

d"^ 

where  d  is  the  electrode  separation  in  cm  and  [i    is  the  mobility. 
The  derivation  assumed  the  electron  velocity  is  proportional  to 
the  field,  i.e.  the  mobility  is  field  independent.   This  result 
means  that  the  current  is  proportional  to  the  square  of  the  anode- 
cathode  potential.   A  derivation  for  a  cylindrical  situation  under 
the  same  assumptions  yields  the  same  dependence  on  plate  voltage. 
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Although  a  general  treatment  has  not  been  attempted,  the  similarity 
of  the  derivation  to  that  of  the  Languir -Childs  law  suggests  that 
the  result  is  geometry  independent.   With  the  additional  assumption 
that  the  mobility  u  is  inversely  proportional  to  pressure,  the 

current  in  a  mobility  limited  diode  is  given  by 

2 

I  a  ~  (A. 2) 

P 

For    weak   fields,    the  mobility    is    not    field-independent 
L40J    and   the   drift    velocity    is    given   by 

h 

v^a(^)     .  (A. 3) 

In  such  a  case,  it  is  easy  to  show  that  the  diode  current  for 
parallel  plate  and  cylindrical  diodes  is  given  by: 

I  a  -V   •  (A-'^) 

'2. 
P 

The   field-dependent,    the   field-independent    and   any 
intermediate    case  as    well    as    the   transition    from   space   charge   to 


m 


obility  limited  operation  are  included  in  the  general  form 


I  =  Cp  (A. 5) 

where  C,  k   and  k   are  dependent  on  the  situation  and  must  be 
empirically  determined. 
b.   Sensitivity 

To  determine  the  sensitivity  of  the  diode  as  a  micro- 
phone, consider  the  simple  circuit  shown  in  Figure  21.   If  the 
effect  of  the  sound  on  the  diode  is  observed  by  measuring  the 
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V  ^ 


R  ^   V, 


OUT 


Figure  21.       Diode  Circuit.      A   simple   circuit   which    can    be 
used   to   study    sound   with   a   mobility-limited 
thermionic    diode. 


NEON 


4      5    6    7  8  9  10 

PRESSURE    (TORR) 


20 


30      40 


Figure    22 


Diode    Current    in    De|)Gndence    on   Pre.so' 
Typical    Diode    in   Argon   and   Neon. 
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voltage  across  R,  then  the  sensitivity,  S,  of  the  diode  circuit 
can  be  defined  as : 

Using  equation  (A. 5),  it  is  easy  to  shov;  that  this  is  given  by: 


(A. 7) 


•^1 

"dc 

p 

1 

*s 

where  V    is  the  d.c.  component  of  the  output  voltage.   This 
expression  shows  that  the  sensitivity  increases  with  the  d.c. 
voltage  across  R,  in  an  approximately  linear  relation,  until  the 
output  voltage  becomes  comparable  with  the  supply  voltage  V.   The 
signal  voltage  V  produced  by  an  acoustic  pressure  wave  p  is  there- 
fore given  by  ^ 

V  =  Sp.  (A. 8) 

c.   Frequency  Response 

Equation  (A. 7)  is  not  correct  for  short  xvavelengths . 
Implicit  in  the  discussion  thus  far  has  been  the  assumption  that 
the  finite  size  of  the  diode  would  not  affect  results.   Actually 
the  diode  averages  the  pressure  over  its  length.   When  fringing 
effects  are  ignored,  the  diode  current  is  proportional  to  its 
length.   Therefore  a  diode  of  length  bi    contributes  a  current 
5i  given  by  . 

where  I  is  the  current  for  a  diode  of  length  i,    considered  uniform 
The  total  current  consists  of  the  d.c.  component  plus  the  a.c. 
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component  caused  by  a  pressure  disturbance.   The  current  fluct- 
uation is  given  by 

I  =  f  P  •  (A. 10) 

^  .     ,   ~       ikx-iCL)t    ,.  •  ^  . 

For  a  pressure  wave  given  by  p  =  p  e        ,  whxch  is  propagating 

in  the  axial  direction  of  the  diode,  the  total  current  I  is  given 

by  integrating  equation  (A. 9)  over  the  length  of  the  diode,  wherein 

the  real  part  of  the  pressure  must  be  used. 


o 


cos  kx  dx, 


therefore:  ~        S  sin   k-S  (A. 11) 

T  =  —  D    : —   . 

R  ^o     kii 

^^      ^                                                ^  ^1-  J-  ,   •   *v-    ^      sin(2njeA 
.The  frequency  response  of  the  dxsde  is  therefore   -   )  y. . 

d.   Temper dk tur e  Effect 

The  equations  developed  thus  far  indicate  that  the 
diodes  can  be  calibrated  in  terms  of  absolute  pressures  by  extra- 
polation from  the  zero  frequency  case,  i.e.  the  d.c.  characteristic 
of  the  diode.   This  requires  that  the  pressure  be  varied  and  the 
current  measured  as  a  function  of  plate  voltage  at  the  pressures 
used.   The  pressure  measurement  is  usually  made  at  some  point 
removed  from  the  diode.   Since  the  filament  is  hot,  the  effects 
of  the  thermal  gradient  between  the  diode  and  the  measuring  point 
must  be  considered. 

Many  properties  of  particles  in  gases  which  are  fre- 
quently expressed  as  functions  of  pressure  are  in  fact  primarily 
dependent  on  the  particle  density.   In  such  cases,  temperature 
effects  can  be  compensated  for  by  correcting  the  pressure  j'cr 
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temperature.      For    example,    mobility    is    given   as    a    function   of 
pressure    even    though    it    is    a    function    of    density. 

If   the    data    is    for    pressures    at    temperature  T     and   the   temperature 
of   the  point    of    concern    is    T    ,    and   the   pressure    is    p    ,    then    the 
effective  pressure,    p    ,    for    entry    into    the    data    is 

Now  consider    the   effect    of   the    temperature   difference 
on    the    diode   characteristics,    starting    from  Gauss's    law: 

V^    =    -   ^  (A. 12) 

o 

and  the  expression  for  current  density 

j  =  pv. 
Assume  electrons  moving  in  a  gas  under  a  weak  field  so  that 

where  C,  is  a  constant  peculiar  to  the  gas.   Consider  cylindrica] 
symmetry  and  negligible  end  effects.   Since  the  current  I  is  a 
constant,  and  related  to  j  by  geometry,  the  density  can  be  expressed 
in  terms  of  the  current  and  other  factors  then  substituted  into 
Eq.  (A. 12).   The  result  can  be  written  in  cylindrical  coordinates  • 


as 


h 


2_  ,_d_  /r  dv  )  ^  L.  RZ  1 

r  dr  (^   dry    2ntZ      c    /i  dV  *>  '2 
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This    can   be   re-expressed  as 


k'2i/  -1,  ,\  -1  ^ 


dV    -^  _d_      r    dV   K^  _J_     HJ: 
v^      dry      dr    \^      dr  j        2n^      C^  l^-ij) 

The  pressure  p   can   be   expressed    in    terms    of    the 
pressure,    p    ,    at    the   measuring   point   as 

P   =  Pi/(T(r))    =  p/(r)     . 

With    this    substitution,    Eq.     (A.l6)    can    be    integrated 


dry  2    2n£      C^    )   ^     ^       ^ 


Therefore 


dr  2   2tt£    r       ' 


F(r )    =    I         /  2(r  '  )r  '  dr 


from   which   the   result 


, =^  v3/4r°  ^ii^j'^' 


3Pi"  ^o 

can    easily    be   obtained.      This    can    be  written   as 

3/2 
V 

I   =  A  r   . 

Pi 

A  review  of  the  steps  taken  to  arrive  at  this  result  shows  that 
even  had  some  other  dependence  on  pressure  been  used,  the  functional 
dependence  does  not  change  because  of  a  temperature  effect. 
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A   theoretical    treatment    has    given   several   predictions 
as    to    the   operation    of   a   inobility-liraited   thermionic    diode   as    a 
microphone.      The   results    may    be   summarized  as    follows. 

1.  The   diode    d.c.    current   will    vary   with  pressure  and  plate 

voltage   as  v      v 

I    =   Cp         V 

•     where        %  ^    k,    ^    1    and  ^  ^    k^   ^    2 . 
1  2  2 

2.  For    a    simple   circuit    like   Figure     21,    the    sensitivity,    S, 
is    given   by 


s  = 

\      V 

P                    ^DC 

3.      The   sensitivity   as    a    rr.icr  cphonc   will    fall    off   v.'ith  '.wave- 
length  as 

sin    (2n^/X) 
(2nX/A) 

2.   Experiment 

An  experiment  was  performed  using  several  diodes  which 
investigated  the  first  two  predictions.   The  frequency  response 
was  not  tested. 

a.   Experimental  Arrangement  and  Procedure 

The  diodes  used  were  of  similar  construction.   They 
consisted  of  a  tantalum  cylinder,  1  cm  long  by  1  cm  diameter,  as 
plate  and  a  tungsten  or  thoriated  tungsten  filament  as  cathode. 
The  diodes  were  mounted  several  decimeters  from  the  tube  ends  and 
the  sound  source,  in  a  tube  4  cm  in  diameter.   A  current-modulated 
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glow  discharge  acted  as  sound  source.   The  a,c.  signal  was 
measured  using  a  circuit  like  Figure  6  with  a  differential  amp- 
lifier to  amplify  the  signal  which  was  then  observed  with  an 
oscilloscope  or  with  a  boxcar  integrator. 

After  outgassing  the  system  and  activating  the  fila- 
ments, space  charge  limited  operation  in  a  high  vacuum  was  verified, 
This  was  done  by  measuring  the  plate  current  as  a  function  of  plate 
voltage  to  a  value  of  30  volts,  which  was  twice  the  highest  ex- 
pected operating  voltage.   If  a  log  I,  log  V  plot  gave  a  straight 
line,  the  selected  operating  point  was  considered  valid.   If  not 
and  the  current  dropped  below  a  straight  line  at  30  volts,  then 
the  diode  was  considered  emission  limited  and  appropriate  steps 
taken.   If  the  operatirjg  point  was  valid,  then  that  filament 
voltage  was  maintained  throughout  the  calibration.   For  the  re- 
mainder of  the  calibration,  the  pressure  was  changed  and  the  I 
versus  V  measurement  made.   For  this  portion  of  the  calibration, 
if  the  log  log  plot  showed  that  the  current  at  higher  voltages  was 
above  the  linear  fit  to  the  lower  values,  than  ionization  was 
taking  place  and  subsequent  operation  below  these  values  of  plate 
voltage  was  required. 

After  calibration,  the  a.c.  response  to  various  pres- 
sure waveforms  was  noted.   The  sensitivity  was  measured  with  a 
boxcar  integrator  set  to  read  the  same  point  in  the  signal  as  the 
d.c.  signal  was  varied  by  changing  the  resistance  of  the  load 
resistor.   The  sound  source  amplitude  was  kept  constant  during  the 
measurements . 


lifl 


b.   Noise 

The  diode  was  found  to  be  susceptible  to  electrical 
noise  from  external  and  internal  sources.   The  external  noise  v;as 
reduced  to  an  insignificant  level  by  shielding  with  heavy  copper 
foil.   The  foil  was  wrapped  around  the  tube  at  its  surface,  ex- 
tending beyond  the  diode  approximately  5  cm  at  both  ends.   The 
shield  was  grounded  with  a  separate  head.   All  leads  to  the  diode 
were  also  shielded. 

In  addition  to  the  external  noise,  the  diode  v;as  sub- 
ject to  an  internal  noise  which  seemed  due  to  spurious  ionization. 
This  noise  v;as  characterized  by  fast  rising  positive  spikes,  re- 
presenting a  sudden  increase  in  current,  followed  by  a  more  gradual 
fall.   These  spikes  always  appeared  when  contamination  was  suspected, 
r.   D.C.  Characteristics 

Dayton  et.al.  LSbI  reported  results  for  a  diode  oper- 
ating in  N„  whereas  Argon  and  Neon  were  used  in  this  experiment. 
The  results  obtained  were  similar  to  those  of  Dayton.   A  typical 
plot  for  Argon  and  Neon  is  shown  in  Figure  22.   The  values  of  k 
and  k   obtained  are  displayed  in  Table  IV  along  with  information  on 
N„  taken  from  Dayton's  paper.   The  values  shown  are  the  averages 
computed  from  data  for  several  diodes  for  the  pressure  range  5  to 
30  torr . 

TABLE  IV 


Gas 

^1 

^2 

^ 

Ne 

.60  +   .01 

1.67  +  .04 

1.5 

Ar 

.63  +    .02 

1.03  +  .09 

1 

No 

.62 

1.98 

1.7 
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d.   Sensitivity 

A  plot  of  signal  voltage  versus  d.c.  voltage  for  con- 
stant acoustic  pressure  constitutes  a  measurement  of  the  sensitivity 
in  arbitrary  units.   Figure  23  shows  the  results  of  such  a  plot. 
The  solid  curve  was  calculated  using  equation  (A. 17)  with  k   and 
the  amplitude  determined  by  a  least  squares  fit  to  the  data  points. 
The  upper  half  of  the  points  were  neglected  in  the  fit  because  of 
a  definite  shift  in  the  slope  of  a  linear  plot  of  the  data  when  the 
plate  voltage  became  comparable  to  the  filament  voltage. 
3.   Discussion  of  Results 

The  value  of  k   is  very  similar  in  all  three  gases  but 
there  is  a  large  difference  in  the  values  for  k  .   Since  the 
theoretical  treatment  assnmes  that  the  pressure  dependence  arises 
from  the  similarity  parameter  E/p ,  if  k   is  the  same,  then  k 
should  be  the  same.   As  this  is  not  the  case,  and  the  analysis 
has  also  shown  that  the  temperature  gradient  should  not  affect 
this,  then  there  must  be  a  pressure  dependence  other  than  E/p. 

Von  Engle  C40j  points  out  that  molecular  gases,  such  as 
N„ ,  have  large  in-elastic  collision  cross  sections  which  causes 
their  mobility  to  be  constant  and  their  drift  velocity  to  vary 
linearly  with  E/p.   This  is  related  to  the  manner  in  which '>^,  the 
average  fraction  of  electron  energy  lost  in  a  collision,  increases 
with  E/p.   Until  inelastic  collisions  start  to  take  place,  ^  re- 
mains constant,  afterwards  it  increases  with  E/p.   The  relation 
between  the  dependence  of  JC  on  E/p  and  the  dependence  of  the 
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V^  (VOLTS) 


Figure  23.  .  Diode  Sensitivity  in  Dependence  of  d.c. 
Signal.   Argon  at  pressure  of  7.35  torr . 
Vertical  scale  is  arbitrary.   The  ordinate 
is  the  d.c.  component  of  the  signal  voltage 
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mobility-limited   diode   current    on   plate  potential    is    a    complicated 
one  and  would   be    difficult    to    establish    quantitatively.      However 
there    is    a    correlation   between    the    experimentally    determined   values 
of   k      and    the   slope   of   the    log  ^    -    log    E/p    curves    as    figure   63    of 
Von   Engle   [4ol.      For    simplicity,    k      is    defined  by 


^  -  (§1 


S 


which  is  used  merely  for  comparison.   The  values  of  k   determined 
from  Von  Engle  are  shown  in  Table  ivfoi^  comparison. 
4.   Conclusion 

The  mobility-limited  diode  nicrophone  should  serve  to  fill 
a  need  in  plasma-acoustic  studies.   It  is  simple  and  inexpensive 
to  construct  and  requires  no  complicated  circuitry.   By  virture  of 
changing  the  d.c.  signal,  sensitivities  of  different  diodes  can  be 
matched  for  direct  comparison  of  signals.   It  is  sensitive. 
Sensitivities  higher  than  .27  mV/microbar  (3^0  mV/torr)  at  5  torr 
ambient  pressure,  can  easily  be  attained.   It  is  transparent  to 
sound  and  offers  little  interference  to  light  in  the  axial  direction 
Calibration  at  zero  frequency  is  simple  and  easy  to  perform. 

It  is  susceptible  to  noise  from  contamination.   This  is  a 
drawback  which  conventional  microphones  do  not  have.   Its  use 
would  be  limited  to  inert  gases  or  those  where  chemical  activity 
with  the  hot  tungsten  filament  would  not  be  expected.   Even  with 
these  drawbacks,  a  diode  microphone  can  be  used  for  many  experi- 
ments where  other  techniques  can  not. 
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APPENDIX  B 

THE  USE  OF  BUTT  SEALS  FOR  FEEDTHROUGHS 
In  order  to  seal  electrode  supports  and  feedthroughs  into 
glass  with  a  minimum  disturbance  of  the  glass  surface,  a  special 
procedure  was  required.   Kovar  wire,  70  mil.  diameter,  was  used 
for  all  supports  and  feedthroughs.   Tantalum  sheet  5  mil.,  was 
used  for  the  electrodes.   Diode  filament  leads  were  wrapped  with 
fine  rhenium  wire  for  ease  in  spot  welding.   The  electrode  supports 
went  along  the  electrode  in  the  axial  direction  and  were  spot 
welded  in  several  places.   In  addition,  a  small  cut  in  the  elec- 
trode was  made  at  the  edge  and  the  support  wire  was  bent  back  and 
crimped  to  hold  the  electrode  in  a  mechanical  grip.   The  small 
cut  was  made  to  prevent  the  wire  from  extending  beyond  the  edge  of 
the  electrode.   The  materials  were  cleaned  prior  to  assembly  and 
the  components  were  cleaned  again  just  before  sealing  in  accordance 
with  procedures  recommended  by  Rosebury   L56J. 

For  final  assembly  and  sealing,  the  components  of  an  element 
were  assembled  on  a  special. jig.   The  jig  was  designed  to  hold  the 
support  wires  and  feedthroughs  in  a  plane  which  would  be  normal  to 
the  axis  of  the  glass  tubing.   After  assembly  of  the  components  the 
jig  was  attached  to  one  piece  of  the  glass  tubing,  earlier  mounted 
in  the  lathe.   The  jig  was  placed  so  that  the  feedthroughs  were 
flush  against  the  edge  as  shown  in  Figure  24.   The  edge  used  for 
the  seals  worked  best  if  it  was  diamond  cut.   Fire  cut  edges  were 
susceptible  to  leaks.   The  second  piece  of  tubing  was  matched  with 
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the  edge  of  the  first  as  shown  in  Figure  25.   A  low  pressure 
purge  of  nitrogen  was  used  to  prevent  the  flame  from  penetrating 
the  tube  and  contaminating  the  surfaces.   After  the  purge  was  on 
and  the  two  pieces  and  jib  were  turning  in  the  lathe,  heat  was 
applied  and  the  two  ends  pushed  toward  each  other.   When  suffi- 
ciently hot,  the  glass  vjould  flow  around  the  leads  and  the  two 
ends  would  butt  together,  forming  a  butt  seal.   A  small  ring  of 
7052  was  then  beaded  around  the  base  of  the  leads,  for  insurance. 
The  complete  tube  was  annealed  in  an  oven  following  the  recommended 
temperature  cycle.   For  the  tube  used  in  the  experiment,  due  to  the 
nature  of  the  main  discharge  filaments,  annealing  was  done  with 
the  tube  being  evacuated   by  a  diffusion  pump.   During  pre- 
liminary work   Several  seals  were  not  annealed  in  the  oven   but 
merely  flame  annealed  on  the  lathe  with  no  later  cracking. 
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Figure  24.   Diode  in  Process  of  Assembly, 
Showing  Placement  in  Jig. 


Figure  25.   Diode  in  Process  of  Assembly, 
Showing  Assembly  in  Lathe. 
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APPENDIX  C 


TAPE  RECORDER  SPECIFICATIONS 


FM  Record/Reproduce 
Input  Sensitivity:   (High  Level) 

(Low  Level) 

Input    Impedance:  (High   Level) 

(Low  Level) 


1    to  20    volts   peak-to-peak 
for   _+  ^0%   deviation 

0.030   to   1    volt  peak-to-peak 
foir   _+   40%   deviation 

20k   ohms    or    greater ,    unbalanced 
to   ground 

50k    ohms    or    greater ,    unbalanced 
to    ground 


Carrier    Frequency,    Response,    Signal-to-Noise   Ratio; 


High  Level  Center    Frequency 

.375    ips  500   hz 

3.75    ips  5        khz 

^7 .  "5      iDS  50      khz 


Response 
DC   to    100    hz    +    1    db 
DC   to    1    khz  jf    1    db 
DC   to    10   khz    +    1    db 


S/N 
35  db 
38  db 
42  db 


Low  Level 
.375  ips 
3.75  ips 
37.5   ips 

Gain  Drift 


500  hz 
5   khz 
50   khz 


Center  Frequency  Drift 

Linear  ity : 

Output  Level: 

Output  Impedance: 
Harmonic  Distortion: 


5  to  100  hz  +  1  db  35  db 
5  hz  to  1  khz  JH  1  db  38  db 
5  hz  to  10  khz  +  1  db  42  db 


2%   or  less  of  fall 
scale 

2%   or  less  of  full 
scale 


(temp,  range 
350  to  110°F, 
line  voltage 
range  110 
vac  to  125 
vac) 


1%   of  best  straight  line  through 
zer  o 

2.8  V  p-p  for  +_   40%  deviation 
into  1 ,000  ohms 

50  ohms  or  less 

1.5%  total  (signal  1  v  p-p, 
3,000  hz,  at  37.5  ips) 
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Output  Filter 
Cutoff  Frequency: 
Frequency  Response: 


Overshoot : 

Transport 
Tape   Speeds  : 
Speed  Accuracy: 


100   hz,     .1    khz,    10   khz,    100   khz 

Filter    set    for    flat    frequency 
response,    _+    1    db   --    Filter    set    for 
flat   phase   response   -4    to   -8    db 
at    cutoff. 

Filter    set    for    flat    frequency 
response,    15%  max.    --    Filter    set 
for    flat   phase   response,    2%  max. 


37.5,    3.75,    0.375  ips 
+   0.2%,    long-term 
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COMPUTER  PROGRAMS 


le  DATA  REDUCTION 

AFTER  THE  DATA  WERE  STORED  ON  TAPE  IN  32-BIT 
BINARY  FORM,  FIVE  PROGRAMS  WERE  USED  SEQUENTIALLY  TO 
PROCESS  THE  DATA.   DISK  STORAGE  WAS  USED  FOR  INPUT  AND 
OUTPUT  EXCEPT  TH^  FIRST  PROGRAM  USED  TAPE  FOR  INPUT  DATA 
EACH  PROGRAM  MAY  BE  THOUGHT  OF  AS  A  STEP  IN  THE  PROCESSING 

l.STEP  ONE    — AVERAGE  — 

C  THIS  PROGRAM  TAKES  THE  RAW  DATA  AND  BEGINS  THE 

C  PROCESSING.   THE  DATA  HAS  BEEN  PRESCREENEO  AND  INFORMATION 
C  STORED  IN  AN  ARRAY  INDICATING  WHICH  SETS  OF  WAVEFORMS  ARE 
C  ACCEPABLE  FOR  AVERAGING.   THIS  ARRAY  IS  CALLED  «KF'.   ALL 
C  ACCEPTABLE  WAVEFORMS  ARE  AVERAGED  AND  THE  STANDARD  DEVIA- 
C  TION  OF  EACH  DETERMINED  FOR  EACH  DIFFERENT  EXPERIMENTAL 
C  RUN.   THIS  INFORMATION  IS  STORED  FOR  USE  BY  THE  NEXT 
C  PROGRAM. 
C 

c 

DIMENSION  DATA (3600) tAVG{3600) tKSTAR(12) tKF{i00tl2) 

DIMENSION  STDEV(3600),DEV(3) 

REWIND  4 

REWIND  2 

DATA  KF,NF,KSTAR/ 1200*0, 1,3*2,3,1,1,3,1, 1,2/ 


C 

c 


c 
c 


READ(3)KF 

5  NR=0 

KG=0 
2    DO    7    I=1,360C 
STDFV(n=0. 

7  A^'G  (  I  *  ='^    (^ 

i   rjR=NR  +  i 

READ(2,ERR=1,END=50)    DATA 
IF(KF(NR,NF).Eq.O)     GO    TO    1 

DO    A     1^1,3600 

STDEV(  I  )  =  STDEV(  I  )  +  DATA  (  I  ) -^^^^2 
4    AVG( I)  =  AVG(I )+DATA(  I) 
KG  =  KG-»-l 
GO    TO    1 

50    DIV=KG 

DO    6    1=1,3600 

STDEV(  I  )==SQRT(  (  STDEV(  I  )-AVG(  I  )*«2/D  I V  )  /  (  DI V-1 .  )  ) 

6  AVGd  )=AVG(  I)/DI  V 

DO    8    J=l,3 
DEV{J)=0.0 
DO    9    1=1,1200 
II=I+( J-l)-1200 
9    DEV( J)=DEV( J)+STDEV{ II )^*2 

8  DEV( J)  =  SQRT(DEV( J) /1200o  ) 

WRITE (4)  AVG,DEV 

WRITE(6,20'J)  N!= 

WRITE (6,203)  OEV 

WRITE{6,2U1)  AVG 

WRITE(6,202)  (KF( I ,NF) , I=1,NR) 

NF=N-+1 

IF(NF.LE,10)  GO  TO  5 

W.RITE(13)  KF 

2C0  FORMAT('l  FILE  NO.  ',13) 

201  FORMAT  (60(  1  X,  1  :)F12  .  4,  /  )  ,/) 

202  FORMAT( '0» ,10( 1CX,10I5,/) ) 
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203  FORMAK'  AVERAGE  DEVIATIONS   'fSFlO.S) 

STOP 
END 


2. STEP  TWO    — SCREEN  AND  RE-AVERAGE — 

C  THIS  PROGRAM  USES  THE  AVERAGE  WAVEFORMS  AND  THE 

C  STyiNDARD  DEVIATION  C/\LCULATED  IN  THE  PREVIOUS  PROGRAM  TO 
C  SCREEN  OUT  DATA  WHICH  LIES  OUTSIDE  OF  TWO  STANDARD  DEVIA- 
C  TIONS  OF  THE  AVERAGE.  IT  RECOMPUTES  THE  AVERAGE   AND 
C  STANDARD  DEVIATION  WITH  THE  SCREENED  DATA. 


C 

c 


c 

C 


DIMENSION  DATA (3  600)  ,AVG(3600) ,KSTAR(12)  ,KF(10C,12) 
*KG(36C0) ,DEV(3) ,STDEV(3600) 
RFWIND  2 
REWIND  4 
NF=1 

READ(3)  KF 
REWIND  3 
STD=2. 

5  NR^O 

READ (4)  AVG»DEV 
DO  2  I=l,36ro 
STDFV(  I  )=AVG(I  )^'*2 

2  KG(I)=1 
1  NR=NR+1 

READ(2,ERR=1 ,END=50)  DATA 
IF{KF(NR,NF)«EQ.OO)  GO  TO  1 

KN=0 

DO  ?  1=1,3600 

J=T  n 200  +  1 

COM=AVG( i )/DIV 
TEST=ABS(CaM-DATA( I ) ) 
TEL  =  STD^*DEV{  J) 
IF(TeSTcGT,TFL)  GO  TO  3 
STOEV(  I  )=STCEV(  I  )+DATA  (I  )*-^2 
AVG(  n=AVG(  I)  +  OATA(  I  ) 
KN=KN+1 
KG(  n=KG(I  )  +  l 

3  CONTINUE 
K^{NR,NF)=KN/36 
GO  TO  1 

50  DO  6  I=1,36CC 
DIV  =  KG(  I  ) 

STnEV(  !)  =  (  STOEV(  I)-AVG(  I  )  *--''2/DI  V  )  /  (  D  IV- 1 .  ) 
1F(STDEV( I ).LT. 0.0)  WRITE (6,203)  I,STDEV(I) 
STDEV( I )=SQPT( ABS( STDEV( I ) ) ) 

6  AVG( I)=AVG( I )/DIV 

DO  8  J=l,3 

DEV( J)=0,0 

SUM=C. 

DO  *^  1  =  1,1200 

I  I  =  I  +  ( J-1)*120C 

SUM  =  SUM  +  KG(  I  I  ) 
9  DEV(J)=0-V(J)  +  STDEV(II  )=<=«2*KG(  II  ) 
8  DEV( J)=SQRT(DEV( J)/SUM) 

WRITE(9)  AVG,DEV 

WRITE (6,200)  DEV 

WRITE (6, 201)  AVG 

WRlTE(6,202)  (KF( I,NF) ,I=1,NR) 

NF  =  NF-H 

IF(NF.LE.IO)  GO  TO  5 
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WRITE(13)  KF 


200 
201 
202 
203 


FORMAK'l  DEVIATIONS    ARE 

FORMAT ( 60{lXrlOF 12. 4t/) »/) 


S3E12.3) 


FCRMAT( 
FORMAT ( 
STOP 
END 


'  n  •  »  1  n  (  1 C  X ,  1  01  5  ,  /  )  ) 


STDEV  SQUARED  IS  NEGATI VE • , 15 t E 12. 2 ) 


C 

c 
c 
c 
c 
c 
c 
c 


c 

c 

c 
c 


3.  STEP  THREE    — DELAY  ,  OVERLAY  AND  ZERO— 

THIS  PROGRAM  TAKES  THE  AVERAGED  WAVEFORMS  AND 
DELAYS  THEM,  REDUCING  THE  NUMBER  OF  POINTS  FROM  1200  TO 
600  IN  EACHo   IT  ALSO  ZEROES  THE  BASE  AND  OVERLAYS  A 
FALSE  SIGNAL  DUE  TO  AN  ELECTROSTATIC  INTERACTION.   THE 
RESULTS  ARE  STORED  FOR  THE  NEXT  STEP,  WHICH  IS  THE  FOURIER 
TRANSFORMATION. 


NSION  DA(12C0,3), IND(3,10) , SI  0(600,3)  , DE V( 3 ) , AVG (3 

IP0INT(2,3,  10) 
VALENCE  (DA( 1,1) ,SIG(1  ,1) ) 

I  PO I NT/C; ,  600 ,  0  ,  600  ,  0 ,  600 ,  240 ,600,120,270,0,288,24 
8 8, 0,2 8 8, 2 40, 600, 12 1,238, 0,28  8,240, 600,204,27 6,  0,2 
600,216,28  8,0,2  64,240,600,10  8,275,0,2  50,240,6  00,10 


DIME 
1 

EQUI 

DATA 
*96,2 

*240,    .    ^    .  .    . 
*0, 240,240, 600, 120, 220, 0,240, 240, 600,120, 240, 0,240/ 

KF=0 

NF  =  0 

RE 'a' I 

READ 


301 


311 
312 


302 
304 


322 

321 


324 


303 
323 


:nd  4 

(5,100) 


(  (IND( J,I) , J=l,3)  ,1=1,10) 


NF=Nr-  +  l 

READ(4,END=3C5)  DA,DEV 

KF=KF+1 

IF(NFoEQcl)  GO  TO  312 

A=URN(0) 

DO  311  1=100,2  40 

JI=130.*URN(l)+2  40 

DA(  I, 2)=DA( J  I, 2) 

CONTINUE 

DO  304  J=l,3 

KC=IND( J,KF  ) 

DO  302  1=1,60  0 

SIG( I, J)=0A(KC+I,J) 

CONTINUE 

DO  321  J=l,3 

AVG(J)=C.0 

IP1=IP0INT( 1,J,NF)- 

IP2=IP0IfJT(  2,  J,NF) 


DO    322    I=IP1,IP2 

AVG( J)=AVG( J)+SIG( I 

DIV=IP2-IP1 

AVG( J)=AVG( J)/DIV 


J) 


305 
IOC 


no    323    J=l,3 

DO    324    1=1,600 

SIG( I ,J)=SIG( I , J)-AVG(J) 

WRITE (6, 200)    KF,J 

DO    303    1=1,589,12 

WRITE (6, 202)     (SIG( I+KM-1,J) ,KM=1,12) 

CONTINUE 

WRITE(8)     SIG    ,DEV 

GO    TO    301 


STOP 

FORMAT (121 5) 
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200  FORMAK'l  R.TRIPP         FILE     ',13,' 

201  FORMAT (20A4) 

202  F0RMAT(1X,12F10.4) 
END 


4. STEP  FOUR   — TRANSFORM — 

C  THIS  PRCGRAM  DETERMINES  THE  SPECTRAL  DECOMPOSITION 

C  OF  THE  WAVEFORMS  AND  THEIR  STANDARD  DEVIATIONS,   GAIN  AND 
C  DIODE  SENSITIVITY  ARE  CORRECTED  FOR. 
C 

c 

C  INITIALIZE  VALUES 

C 

DIMENSION  SIG(600t3) t SPECT ( 600 , 3 ) TT(60r) ,FPEQ{60) , 

1  SHO(200) ,BIAS{3) ,GAIN(3),SENS(3) ,IP(3) , 

2  SDEV(4r, 3 ) tSREQ(60) , DEV(3) t INV(2CC)  ,S{20G) 
■  EQUlVALEfJCE  (  S  I  G  t  SPECT)  ,  (  SHO ,  SDEV  ) 

COMMON  SIGtSDEVtSt INV 

DATA  P I tGAI NvNPNTStNF, PERI 0D/3o 14159265, 3*50, ,600,0, 
1     0.00512/ 
DO  300  I=1,NPNTS 

300  T(I )=I 

DO    301     1=1,40 
FREQ(I)=(I-1 )/PERIOD 

301  SREQd  )=SQRT(FREQ(I  )) 
C 

C  READ  DATA  FOR  NEW  DATA  SET 

C 

17  CONTINUE 
NF=NF+1 

READ(8, END=9o)  SIG,DEV 
RFAD(5,101)  SENS 
WRITE(6,101)  SENS 
r 

C  CORRECT  FOR  AMPLIFIER  GAIN 

C  AND  DIODE  SENSITIVITY 

C 

DO  70  5  J =1,3 

FACT=GAIN( J)^SENS( J) 

DEV{ J)=DEV( J)/FACT 

DO  303  I=1,NPNTS 

303  SiGd,  J)=SIG(I  ,J)/FACT 
C 

C  DISPLAY  WAVEFORM  ON  PRINTER 

C 

DO  305  1=1,199 
305  SHn(  I)=(SIG(3>;^I-2,  J)+SIG(3*I,J)+SIG(3*I  +  2,  J)  )/3. 
SH0(2OC )=SlG(6e0, J) 
WRITE(6,200) 
705  CALL  PLOTP(T,SHO, 200,0) 
C 

c 

C  IF  THIS  IS  THE  FIRST  DATA  SET,  DETERMINE  CORRECT 

C      GAIN  FOR  REMAINDER  OF  DATA  SETS. 

C 

IF(Nf^.EQ.l)  GO  TO  313 
C 

C  DETERMINE  SPECTRUM  AND  STD.  DEVIATION 

C 

CALL  SPCTRA 

DO  304  J =1,3 

304  SDEVd,  J)  =  SQRT{SDEV(1,J)=<^*2  +  DEV(J)^'^*2) 

18  WRITE (10)( (SPECTd , J) ,SDEV(I, J) ,J=1,3  ) , FREQ ( I ) , I =1 , 25 
C 

C  DISPLAY  RESULTS  AND  PASS  THEM  TO  NEXT 

C  PROGRAM 

C 


DO  307  J=lj3 
WRITE(6,206) 
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c 
c 


DO    309    1=1,40 

SIG( 1+40, J)=SIG( I, J)-SDEV(Tt J) 

IF(SPECT(I+4Ct J l-LT.O.O)     SPECK  1+40, J)  =  0.0 

SIGCI+RO, J)=SIG( I, J)+SDEV( I, J) 
309  CONTINUE 

DO  308  KJ=1,3 

KiM=121-KJ*40 

CALL  PL0TP(FREQ,SIG(KM,J) ,40,KJ) 
3C8  CONTINUE 
307  CONTINUE 

1=0 

DO  302  K=l,27 

1  =  1+1 

WRITE (6,210)  K, (SPFCT(K,J) ,SDEV(K,J) , J=l , 3 ) , FREQ( K ) , SR 

RATI0=5PECT(K, 5 ) /SPECT ( K,  2  ) 

SPECT{I,2)=ALOG{SPECT(K,2) ) 

SPECT (I ,3)=RATI0 

SPECT (I ,1)=AL0G( RATIO) 
302  CONTINUE 
C 
C 

WRITE (6, 200) 

CALL  PLOT P{ PRE Q, SO ECT (1,3) ,27,0) 

CALL  PL OTP( SREQ, SPECT (1,1) ,27,0) 

GO  TO  17 
C 

C  CALIBRATE  GAIN 

C 

313  CALL  RHARM(SPECT(1,2) ,8,INV,S,IER) 

SPECT  (1,2)=      SgRT(SPECT(53,2)*^'2+SPECT(54,2)**^2) 

GAIN(2  )  =  GAIN(2  )=;'SPECT(1,2) 

CALL  RHAFM(SPECT  (1,3)  ,8,  IN\/,S,  lER) 

SPECT  (1,3)=      SQRT(SPECT(53,3)*=^2  +  SPECT(54,3)**2) 

G.ATN(3»=nAT(^(3  )  ASPECT  (1,3) 

'n'RITF(6y  10  3  )  GAIN;  S^ECTd  J?)  -  SPEC  T(  1.3) 

GO  TO  17 
C 

c 

C  FORMAT  STATEMENTS 

100  F0RMAT(20A4) 

101  FORMAT (3F1C« 5) 
200  FORMAT ( IHl  ) 

202  FORMATdHl,/,'      TIME  AND  SIGtJAL   •  ,  /  ,  50  (  1  2F  10o4  ,  /  )  ) 

203  FnRMAT(  IHl,/, TIC  SPECTRAL  ANALYSIS  OF  SIGNAL  FOR  ', 
I'DIOOE  » ,12,/, T6, ' INDEX   FREQUENCY     COSINE         ', 
2'SINE     AMPLITUDE    PHASE') 

204  FORMAT( 110, 1P3E 1 1. 3 , 0P2F 1 1 . 4) 

206  FORM AT {«  «,♦  UVL=',E10.2) 

207  FOPMAT( '      FILE  %  14) 

210  FORMAT( I5,8E11.3,I 5) 

211  F0R>^AT(20A4) 

500  FORMAT (IH  ,T10,'LESS  THAN  TEST') 
99  STOP 
END 
C 
C 

C      SUBROUTINE  SPCTRA 

C  THIS  SUBROUTINE  DETERMINES  THE  FOURIER  TRANSFORMS 

C  Of^  THE  WAVEFORMS  FOR  FIVE  START  POINTS,   IT  THEN  AVERAGES 
C  THEM  AND  CALCULATES  THE  STANDARD  DEVIATION  OF  EACH  OKE, 
C 
C 

SUBROUTINE  SPCTRA 
IMPLICIT  CGMPLEX'^8(0) 

DIMENSION   SPEC (516,3) ,S(200), I NV( 200 ) , M ( 3 ) , A (2, 2 58, 3) 
1  SIG(600,3),SDEV{40,3) ,AVG(40,3  ) 

COMMON  SICSOEV,  S,  INV 
EQUIVALENCE   (SPEC, A) 

155 


c 
c 


c 
c 


c 
c 


DO    201    I =1t^0 
DO    20  2    J=l»3 
AVG( I  5 J)=0.0 

202    SDEV( I, J)=r.O 

201    CONTINUE 


DO    706    KT=1,5 
NL=10*(KT-1 ) 
DO    705    J=lt3 
00    707     1  =  1 1 512 
SPEC(  I,  J)  =  SIG(f\iL  +  I,  J) 
7C7    CONTINUE 

CALL    RHAPfMSPEC(  l,J)t8tINV,S,IER) 
IFdER.NE.O)    WRITE(6,100)     lER 

705  CONTINUE 

DO  302  J=l ,3 

DO  301  1=1,40 

RSPEC=A(lv  I»J)*'«^2+A{2t  I,  J)**2 

SDEV( It J)=S^EV(  I,  J)+RSPEC 

301  AVGCI , J)=AVG(I t J)+SORT(RSPEC) 

302  CONTINUE 

706  CONTINUE 

DO  304  J=l,3 

DO  303  1=1,40 

SDEV( It J)=  ABS( (SOFV{ I, J)-AVG( I, J)**2*o2) ) 

SDEV(  I,  J)  =  SOPT(  SDEV(  I  T  J)  )'"^-.2  5 

303  SIG(  I, J)=AVG(I , J)*.2 

304  CONTINUE 
RETURN 

100  FORMAT ('    HARM  ERROR    IERRT='tI3) 
FKIQ 


5. STEP  FIVE   — TRANSFER  FUNCTIONS  — 

C  THIS  PROGRAM  TAKES  THE  FIRST  25  COMPONENTS  OF  THE 

C  SPECTRAL  OECOf^POSITIGN  OF  THE  WAVEFORMS  AND  USES  THEM  TO 
C  OBTAIN  THE  ACOUSTIC  TRANSFER  FUNCTION.   IT  ALSO  DETERMINES 
C  THE  SOURCE  TRANSFER  FUNCTION  AND  DISPLAYS  THE  DATA  ALONG 
C  WITH  A  CURVE  REPRESENTING  THE  PREDICTED  FORM  ,     WHICH  WAS 
C  FITTED  TO  PASS  THROUGH  ONE  POINT.   THE  STANDARD  DEVIATIONS 
C  OF  THE  RESULTS  ARE  ALSO  CALCULATED. 
C 

IMPLICIT  REAL*8  (  D  )  t  COMPLEX'^^e  (Q) 

REAL>:^8  TITLE(lO)/7*'  't'TRIPP    ',2-^'  '/ 

DIMENSION  DX( 25) tDY(25) ,DF(25) ,DW(25) tDELY(25) , 

1  DE(20) ,DSB(20) tX(3),SPECT(25,3) tFREQ(25), 

2  XX(2),YY(2),SDEV(25t3),Z(25,3),UL(2), 

3  SREQ( 80) tYZ{80) tUSRCE(25) 
C 

DATA  DW , X/ 2  5* 1 . DO  t  O . 0  t  98 . 7  t  251 . 3/  , 
1      XX,YY/C.0,4500.t10.»0o0/ 
K=25 
NF  =  0 
10  CONTINUE 
NF=NF+l 

READ(1^»END  =  99)  ( ( SPECT ( I , J )  t SOEV (I » J ) , J  =  l ,  3  )  , 
I  FREQ( I ) ,1=1,25) 

DO  300  1=1,25 
300  USRCE( I )=SPECT( I ,3)/SPECT(I ,2) 

WRITE (6, 20  3)  (USRCE( I ) ,FREQ( I ), I =1,25) 
C 

C  CONVERT    TO    DOUBLE    PRECISION 

C 

L=0 
DSUM=0.0 
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DO  301  I=2tK 

FSPQ=(SDEV(I,2)/SPECT (1,2) )**2 
1     +(SDEV{ I,3>/SPECT(I,3) )**2 
L  =  L  +  1 

DY(L)=ALCG(USRCE(  I)  ) 
DW(L»=1.D0/FSPQ 
DSUM=DSUM+DW(L ) 
nX(L)=SQRT(FREQ( I ) ) 
301    CONTIivJUE 

IF(L.LT.4)     GO    TO    10 

DO  311  1  =  1  J L 
311  DW( I )=DW(I i^DL/DSUM 
C 

C  DETERMINE  ALPHA 

C 

CALL  LSQPL2( L,l ,DX, DY, OW t OF , DELY, DB , DSB ,T ITLE ) 

ALPHA  =  -SN'GL(DB(  2))/(X(3)-X(2)) 

UALPHA  =  SNGL{0SB(2)  ) / ( X ( 3 ) -X ( 2  )  ) 

WRITF(6,205)  ALPHA, UALPHA 
C 

C  DISPLAY  LEAST  SQUARES  FIT 

C  DISPLAY  TUBE  TRANSFER  FUNCTION 

C 

ZINC=FRE0(K)/79. 

FA=DB{2) 

YZ(1)=DB(1) 

YZ(1)=FXP{YZ(1  )  ) 

SRFC( 1)=C, 

DO  30  4  I=2,8C 

SREQ{  n=(I-l  )*ZINC 
304    YZ( I ) =YZ(1)^EXP(F ASSORT (SREQC I) ) ) 
C 

CALL    PLOTP(SRE0TYZT&0,n 

CALL    PLCTP(FREQTUSRCEt25,3) 

GO    TO     ( 11 , 1 C , 1 C , 1 0 , 10 , 1 1 , 10 , 1 1 5  1 0  )     ,  NF 
1 1    rnMTj  wiiE 
C       ---•■■-  DETERMINE    SOURCE    AND    SOURCE    UNCERTAINTY 

C 

UL{  1)=('. 

UL(2)=G. 

WRITE(6,?0A ) 

Al  =  2=^l«3/152,6 

A2=-98.7/l?2o6 

DO    30  3     1=^1  ,K 

Z(I tl)=SPECT(I ,2)**A1*SPECT(I ,3)** A 2 

USRCEd  )  =  (  Al=f^SDEV(  Tt2)/SPECT(I,2)  )**2 
1  +(A2*SDFV(  It3)/SPECT(It3)  )*»2 

USRCF(  n=Z(  1,1  )-SQRT(USRCE(  I  )  ) 

Z(  I,2)  =  Z(I,n+USPCF{I  ) 

Z(I ,3)  =  Z{I ,1)-USRCE(I  ) 

UL(2)=AMAX1(Z{ 1,2) ,UL( 2)) 

WRITE( 6,201 )  I,Z(I ,1) ,USRCE( I) ,FREQ(  I) 
303  CONTINUE 
C 

C  DISPLAY  SOURCE  SPECTRA 

C 

WRITE(6,202) 

CALL  PLOTP( XX,UL,2, 1) 

CALL  PLOTP( FPEO,Z( 1  ,1)  ,K,2) 

C^LL  PLOTP(FREO,Z( 1,2) ,K,2) 

CALL  PL0TP(FREQ,Z(1,3) ,K,3) 
C 

C  DETERMINE  TRANSFER  FUNCTION  AND 

C  UNCERTAINTY 

C 

YY(1  )=:0.C 

DO  302  I=1,K 

TZ=Z(I,1) 

Z( I,1)  =  7(I,1 )/SPECT(T,  1  ) 

SDEV(  I  ,1  )=Z(  I,  1  )'i-SQRT(  (USPCE(  I  )/TZ)^=f^2 
1  -KSDEVt  I,1)/SPECT(  1,1)  )  ^- 2  ) 
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Z{ I,2)  =  Z(I, 1)+SDEV(  It 
Z(I,3)=Z(I,1)-S0EV(I, 
WRITF{ 6^201)  I,Z(I,1) 


302 


C 
C 

c 


YY(1)^AMAX1( YY( 1) 

CONTINUE 

YY( 1)=YY(1)-1.2 


Z(I 


1) 

1) 

»  SDEV(  1,1)  ,FREQ{ I) 

tl)) 


DISPLAY  TRANSFER  FUNCTION 


WRITE (6,200) 

CALL  PLOTP( XX,  YY,2, 
PLOTP(FREQ, Z( 1 
PLOTP (FRE0,Z(1 
PLCTP( FREQ,Z( 1 
7  5E-04 
DRAVv(K,FPPQ,Z( 


CALL 
CALL 
CALL 
YS=C 
CALL 


1) 
,1) 
,2) 
,3) 


K, 
K, 
K, 


2) 
2) 
3) 


C 

c 
c 


^,2, 5, 6,0, LAST) 
CALL  ORAW(K,FRFQ,Z( 1,2)  , 

2, 2, 5, 4,0, LAST) 
CALL  DRAW(l<,FRcQ,Z{  1,3)  , 

2,2, 5, 4,0, LAST) 


1,1)  ,1,1,T 


!TLE,TITLE,10C0.,YS,0,0, 
ITLE,TITLE,1D00, ,YS,0,0, 
2, 2, TITLE, TITLE, IC 00. ,YS, 0,0, 


,  c , 


DISPLAY  PREDICTED  CURVE  FITTED  TO  DATA 


102 


306 


READ{5,1^.2)  AA 

FORMAT (F 10. 5) 

FACT=ST 

FACT2-Z 

SREQ(l) 

DO  306 

FACT=SI 

Y  Z  (  I )  =  F 

CALL  DR 


305 


<?9 
100 
ICl 
200 
201 
202 
203 

2r'4' 
2C5 


DO  305 

SREO(l) 

SRFQ(2) 

SDEV(2, 

SDEVd, 

KK=I/!<  + 

CALL  DR 

GO  TO  1 
STOP 
FORMAT( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 

■  FORMAT ( 
FORMAT ( 
END 


N(  A 
(  14 

=  lc 

1=1 

n(A 

ACT 
AW{ 
2 
1=1 

=  FR 

=  FP 

D- 
1)  = 
2 

AW( 
0 
0 


A*  1  .  0  1  34E -04^^ F REQ (  1  4  )  )  /  FREQ  (  1  4  ) 
,1)/FACT 


A* 
vF 
80 
,2 
,K 
^0 
FO 
Z( 
7{ 


po 


1<.C134E-04=!=SREQ( 
ACT2 

,SREQ,YZ,2,C,TITL 
,5, 4,0, LAST) 


) 


I  )  )/SREQ(  I ) 
E,TITLE,1000.,YS,C,0, 


(I 
{  T  ) 
1,3) 
1,2) 


2,SREQ{1  ),SDEV,KK,0,  TITLE,  TITLE,  100').  ,YS,0, 
,2, 2, 5, 4,0, LAST) 


20A4) 

I^) 

•1      TRANSFER 

lOX, 13, 3E12»3) 

»  1     SOURCE  SPECTRA'  ) 

»!^  RATIO 

/,25{ 10X,2F12,5,/)  ) 

»'^-'  ,T15,  '  SOURCE 

'0» , '  ALPHA    IS' ,1 


FUNCTION  OF  SOURCE' ) 


FREO' , 

FREO') 

P1E1G.2,'  +/-'  ,Eir;.2,/) 
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II.  INGARD-SCHULZ  ANALYSIS 


C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


THE  «=IRST  PART  OF  THE  PROGRAM  INITIALIZES  VALUES. 
SYMBOLS  FOR  THE  CONSTANTS  OF  THE  GAS  ARE: 

CD    MOLECULAR  DIAMETER 

CMN   MASS  OF  ONE  MOLECULE 

CG    RATIO  OF  SPECIFIC  HEATS 
VALUES  OF  THE  INTEGRALS  USED  IN  EVALUATING  THE 
ELECTRON-NEUTRAL  COLLISION  FREQUENCY  AND  THE  AVERAGE  COSS 
SECTION  ARE  CONTAINED  AS  A  TABLE  OF  ENTRIES  AS  Yl  AND  Y2 
RESPECTIVELY^   EN  IS  THE  VALUE  OF  ELECTRON  TEMPERATURE 
USED  AS  ENTRY  INTO  THE  TABLE.   FOR  THE  COMPUT AT  I ONS »  THE 
MKSA  SYSTEM  OF  UNITS  ARE  USED  EXCEPT  MASSES  ARE  IN  GRAMS 
AND  LENGTHS  IN  CENTIMETERS. 


C 

c 
c 


PEAL 
*NAME 
DIME 
DIME 
CCMM 
CO  MM 
PI  =  3 
CK=1 
CME  = 
CA=„ 
CQE  = 
CE  =  8 


^^ 
(6) 

NSI 

NSI 

ON 

DM 

.1^ 

.38 

9.1 

5  2*^ 

1.6 

.85 


X(5) ,XI (2) tSTARPT(5) , LOWL I M( 5 ) ,UPL IM ( 5 ) t XMN ( 5 ) / 

ON  XXdOO)  ,YY(  100)  ,EE(40  )t  SIG(40) 

ON  P0(5) ,0S(5) 

/ INNER/ C ( 2C ), EN ( ICC) » Yl( 100) ,Y2( 10  0) ,ARG( 10) ,VA 

/BLOCl/OUT{2C)/BLOC2/MO 

159265 

0622E-16 

09558E-28 

17715E-8 

C21917E-19' 

42E-21 

EVALUATION  OF  CONSTANTS 


301 


C 

c 
c 

c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


READ(5, 101)  C D, CMN tCGtCPtNE, NAME 
DO  30  1  I  =1,  PC- 
OUT  (  I  )=^0.0 
C(  I  )=0.0 

C(] )=SORT(CG^CK/CMN) 
'~'21-1'/M-^''  A;^^'iP  T^r.n*''':?  ) 
C(3)=^pi=^-.5*CA=(^'^2 
C{7)^Cf-iE/CMN 
C(9)-2c=^=CK/CME/PI 
C  ( -^  )  -  (  C  G- 1  )  -  C  (  7  )  ^^*2  =5^0  (  9  )  '1==^  1 .  5 
C( 5)-S0RT(8,/PI/CG) 
C(6)-2.*C(3)';^S0RT(C(9)  ) 
C(8)=CG 
C(10)=CR 

C(  11)=CQE/SQRT(CME-^CE) 
C(12)=CQE--^2/4./?I/CE/CK 
C(13)=SQRT(CK/CME) 
C{  1A)  =  SQRT(  1.667=!^CK/CME) 
C(  15)  =  12.''^PI 
WRITE (6 t 201)  C 

.ESTABLISH  INTEGRAL  TABLE 

RE AD (5, 104)  (EN( I ) tYK  I)»Y2(I ) ,  I  =  lt 100) 


ACCE 

EVAL 

THE 

THE 

THE 

A  LA 

A  VA 

GIVI 

OTHE 

ATTE 

THE 


AFT 
SS  TO  T 
UATES  T 
ARRAY  ' 
A^RAY  ' 
EVALUAT 
PGE  QUA 
LUE  OF 
NG  THE 
R  VALUE 
NUATION 
UNCERTA 


ER  I 
H^  S 
HE  U 
OUT' 
OUT  ' 
lONc 
NT  IT 
KCC= 
VALU 
,  TH 
FAC 
INTY 


hJ  I  T  I  A 

UBROU 

ISPEP 

CO  NT 

CGNT 

THE 

Y  or 

1  CAU 
ES  OF 
E  SEC 
TOR  F 
OF  T 


LIZATI 

TINE  • 
SION  R 
A  INS  I 
AINS  V 
CONTRO 
INFORM 
SES  SE 
^LL  I 
OND  PA 
OR  A  S 
HE  RES 


ONt  the  MAIN  PROGRAM  CONTROLS 
INGAPD'.   the  subroutine 

elation  of  ingard  and  schulz, 

ALUES  OF  interest  CALCULATED  IN 
L  VARIABLE  'KCC  IS  USED  TO  GIVE 
ATION  IN  A  SINGLE  CALL  TO  INGARD 
VERAL  NAMELISTS  TO  BE  PRINTFDt 
NTERMEDIATE  QUANTITIES.   FOR  ANY 
RT  IS  USED  TO  CALCULATE  THE  EXPE 
ET  OF  PARAMETERS  AND  DETERMINE 
ULT. 


KCC  =  2 
J  =  3 
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REA0(5,103)     STARPT 

CALL     INGARO    (STARPT) 

WRITE(6,209)    0UT(19) 
C 

C  DETERMINE    ATTENUATION    CONSTANT    FOR 

C  COMPARISON    WITH    EXPERIMENT 


C 


c 
c 


c(  10=2.5 

LR=1 
800  READ(5,103,END=999)  STARPT 
REA0(5,103)  OS 
KCC=I 

DO  310  LQ=1 v2 

IF(LQoGT,i)  STARPT( 1)=6C0. 
WPITE(6,2CB)  STARPT  ,DS 
DO  304  M0=15t17 
CALL  INGARDiSTARPT) 
KCC  =  2 

TO  =  0'JT(MO) 
DO  3?2  L=1t5 
TX=STARPT{L) 
STARPT(L)=TX'!'l.l 
CALL  INGARD(STARPT) 
PG{ L)=(OUT(MG)-TO)/ (.1*TX) 

302  STARPT(L)=TX 
VAR=O.C 

DO  303  L=lt5 

303  VAR  =  VAR+{PG{L)  >"-DS(  L))*«2 
VAR=SORT(VAR) 
WRITE(6t206)  MO,TO,VAR,PO 

304  CONTINUE 
31C  CONTINUE 

IF(LR.GTcl)  C(10)=2, 

LP=LR+1 

GO  TO  SCO 

101  F0RMAT(4E10.0, I2t6Al) 

102  FORMAT (2F10.0) 

103  FORMAT(  SFl^'^^O) 

104  FORMAT (3L16c7) 
1C5  F0RMAT(315) 

201  FORMATdH  ,  / 1  4  (  IP  5E 12.  3  ,  /  )  ) 

202  FORMAT (IHl) 

203  FCRMATdHl, /TIO, 'MINIMUM  FOR  0UT(»tI2T')  IS  •,1PE9.1, 
-3Xt'  AT  'tSElOol) 

204  PORMATdHOtTlO,  •  ST^RT  POINT  WAS  •,lP5ElC.n 

205  FORMAT( IHO, TIC, 'OUTPUT  VALUES  ARE :  '  , / ,4( 1 P5F1 2 .3 t / )  ) 

206  FORMATC  0UT(',I2,M  IS  •  t  E 12.  3t  •  WITH  STD  DEV'fE12.3, 
1  /,'     PARTIALS  ARE  •,5E12.3) 

208  FORMATdH  t1P5E15c3t/) 

209  FORMAT('     SOURCE  TERM   'tlPlE12.3) 
214  FORMAT( 1H1,T3D, • PLOT  FOR  •,6A1) 

<^99  STOP 

END 


C  SUBROUTINE  INGARO(X) 
C 

c 

C  THIS  SUBROUTINE  SOLVES  FOR  THE  REAL  AND  IMAGINARY 

C  COMPONENTS  OF  THE  PROPAGATION  CONSTANT  'K'  OF  AN 

C  ACOUSTIC  WAVE  TP/^VELING  IN  A  CIRCULAR  TUP.E  THROUGH  AN 

C  IONIZED  GAS,  ACCORDING  TO  THE  EQUATIONS  DEVELOPED  BY 

C  INGARD  AND  SCHULZ. 

C 

C  VARIABLES  AND  THEIR  SYMBOLS  ARE: 

C 

C  ANGULAR  FREQUENCY        XCl) 

C  NEUTRAL  DENSITY          X(2) 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


ELECTRON  DENSITY  X(3) 
NEUTRAL  TEMPERATURE  X(4) 
ELECTRON  TEMPERATURE     X(5) 

VALUES  CF  CONSTANTS  ARE  COMPUTED  IN  THE  MAIN 
PROGRA-l  AND  PASSED  IN  COMMON.  LIKEWISE  THE  PRODUCTS  OF 
THE  INTEGRATION  SUBROJTINE  ARE  PASSED  FROM  THE  MAIN 
PROGRAM  IN  CCMMONo   ALSO  THE  RESULTS  OF  THIS  SUBROUTINE 
ARE  PASSED  BACK  TO  THE  CALLING  PROGRAM  IN  COMMOM.  THEY 
ARE  STORED  IN  THE  ARRAY  'OUT'.   THE  INDEX  NUMBERS  OF 
THIS  ARRAY  AND  THE  QUANTITIES  STORED  THERE  ARE: 

-SQUARED 

OF  N-SQUARED 
ROPAGATION  CONSTANT 

OF  PROPAGATION  CONSTANT 
ON  VELOCITY 
ELOCUY 

LECTRON  VELOCITY 
RON  VELOCITY 

)  TO  THAT  WITHOUT  IONIZATION 
ONENT  OF  'K'  WITHOUT  lONIZ. 
-SQUARE  (APRROX.) 

OF  N-SQUARED         (APPROX.) 
ROPAGATION  CONSTANT  (APPROX.) 

OF  PROP.  CONSTANT   (APPROX.) 
NSTANT  WITH  IONIZATION 
NSTANT  WITHOUT  IONIZATION 
NUATION  DUE  TO  IONIZATION  ' 
E  STRENGTH 


1 

REAL  PART  OF  N 

2 

IMAGINARY  PART 

3 

REAL  PART  OF  P 

4 

IMAGINARY  PART 

5 

MAGNITUDE  OF  I 

6 

PHASE  OF  ION  V 

7 

MAGNITUDE  OF  E 

8 

PHASE  OF  ELECT 

9 

RATIO  OF  OUT (4 

10 

IMAGINARY  COMP 

11 

REAL  PART  OF  N 

12 

IMAGINARY  PART 

13 

REAL  PART  OF  p 

14 

IMAGINARY  PART 

15 

ATTENUATION  CO 

16 

ATTENUATION  CO 

17 

CHANGE  IN  ATTE 

19 

ACOUSTIC  SOURC 

c 
c 
c 
c 

c 
c 


C 

c 
c 


SUBROUTINE  INGARD(X) 
IMPLICIT  REAL*8  (DJ 
IMPLICIT  COMPLEX'i^S  (Q) 
RFAi.^4  X(  SI  ,  n  »  17,1  OSS  tDtDI 
INTEGER*4  KLIC/0000/ 
EQUIVALENCE  (0Y1,RQY1) 

CCMMON  /IN^IER/C(  20)  tEN(  100)  »  Yl(  ICO)  ,Y2(100)  tARG(IC)  , 
*VAL(10)  ,KCC, J 
CCMMON  /BL0C1/0UT(20) 

NAMELIST  /LISTA/CN»CL,BULK,WALL,LOSS 

NAMELIST  /LI  STB/SI GMA,TAU»W in, WNI,WEN,WN£tGN, 11,12 

NAMELIST  /LISTC/GN,DtDI ,AA,B1»B2,B 

NAMELIST  /LISTD/WE, WI ,CCP,CLII , CL I , WE  I , W 1 E , CE , TAU I 

NAMEL 1ST  /L I  ST  F/OA  T  QB , QC  t  QD  T  OE , QF  t  QG  T  QH , Q I 

NAMELIST  /LISTF/gYl tQY2,QY3TQY 

NAMELIST  /LISTG/QVI ,QVE 

NAMELIST  /LISTH/RMtRDtTHETA 

CN=C(1)>^SQRT(X(4)  ) 

CL=C(2)/X{2) 

BULK=1,  54919*CL/CN*X{1  )*=*=2 

WALL  =  2o  07444?* SQRT(X(1)*CL''.^CN)/C(1C) 

L0SS=BULK+WALL 

IF(KCC.EQ.l)  WRITE(6,LISTA) 

ENTER  TABLE  OF  ENTRIES  WITH  VALUE  OF  ELECTRON  TEMP 
AND  OBTAIN  VALUE  OF  Y2 .  USE  TO  CALCULATE  TAU. 


CALL  ATSM(X(5)  tEN,Y2aOO,l,ARG,VAL,5) 
CALL  ALI(X( 5)tARGtVAL, I2,5,.1E-2,IER) 
SIGMA=C(3)^I2 
TAU=C(4)'^X(5)*^-=1.5i=X(3)*SIGMA/CN*^2 


ENTER  TABLE  OF  ENTRIES  WITH  VALUE  OF  ELECTr^ON  TE! 
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c 
c 

c 
c 


c 

c 


AND    DBT^IN    VALUF    OF    Yl.    USE 
CALL    ATS.M(X(5)  ,  EN, Yl»  100,1, 
CALL    ALI (X(5), ARG,VAL, 11,5, 
WEN  =  C(  6)'-i=X(2)^^SQRT(  X(5)  )=*!  1 
WNE=X(? )/X(2)*C(7)*WcN 
WIN=C(5)*CN/CL 
WNI=WIN'^X(3)/X(2) 
GN  =  V/NE  +  WNI 
IF(KCC.EQ.l)    WRITE(6,LISTB) 


IF(J.EQ.2)  GO  TO  401 


TO  CALCULATE 
ARG,VAL, 5) 
.lE-2, lER) 


WE.N 


500 


C 

c 


c 
c 


c 
c 


c 
c 


c 
c 


WE=C(11 )*SQRT( X(3) ) 
WI=WE*SQRT(C(7) ) 


RM=C(12)/X(4) 

RD=C(1?)/WE*SQRT(X(5) ) 

THETA=RD/RM 

IF(THETAcGE.U  )  GO  TO  5C0 

I F ( K  L I C . L  E . 2 )  WR I T  E ( 6 , L I STH ) 

THETA=lol 

KLIC=KLIC+1 

CONTINUE 

IF(KCC.EQ,1)  WRITt(6,LlSTH) 


CLII- 
CLI  =  I 
IF(CL 
TAUI- 
*/C( 10 
CCP  =  C 


Ul  P  T  =  ^, 
WIE  =  C 
CE=C( 
IF  (KG 


1./{X(3)';<RM=^^*2*'^L0G{THETA)  ) 

. /(K/CLII+I./CL) 

IcLE.O.O)     WRITE(6,LlSTn) 

1.5492^'^CLI/CN*X(  1  )  ^*2+2.  0  7444*SQRT  (  X  {  1  )  *CL  I*CN  ) 

) 

{ 15)*RD**3*X( 3) 

s  6  2  E- 1 ?-X ( 3 ) *  AL  OG ( CCP ) /X ( 5 ) 

14')*SQRT(  X{5)  ) 

C.EO.l)    WRITE(6,LISTD) 


=  ='^1.5 


CQ=(0, 
CA=1,+ 
QP=-QQ 
QC=(TA 
QD=-QQ 
QE=-(W 
QF=(WI 
QG=-OQ 
QH=(WE 
CI=-(W 
IF(KCC 


5J 

•-Ms' 

u- 
1/ 

/X 

/x 

E/ 

.E 


--MLOSS  +  GN)  /X(l  ) 

NI/X{1  ) 

QQ*VJME  )/X(l) 

IN/X( 1 ) 

X(  1  )  )**2  +  Q0*(TAUl+WIN  +  VnE)/X(  1) 

(1)  )--2--QQ*WlE/X(  1) 

EN/X(1) 

(1  )  )-=^?-QQ*WEI/X(l) 

X(l  )  )--2  +  l,/C  (8)^^(CE/CN)':^*2  +  00-(WEN+V,'EI  )/X(l) 

q.l)  WRITE(6,LISTE) 


QY1  =  QE-QI-QH-~^QF 

ROYl  =  -(  (Wl=fCE/CN)^*2-0.6+{TAUI-»-WIN)*(  WEN  +  WEI  )-WlE^NWEN) 
*/X{  I)'i^--!=2 
0Y2  =  QC^MQD'--0H-gE-QG) 
QY3=QB^MQG--^0F-QD*0I  ) 
QY=0A+(QY2+QY3 )/QYl 
IF(KCC,E0.1)  WRITE (6,LISTF) 

0UT(1)=REAL(QY) 
0UT(2)=AIMAG(QY) 
QS1=CSQRT{QY) 
0UT(3)=RF^L(QS1 )*X(1) /CN 
0UT(4)  =  ALMAG(QS1)=^X(1)/CN 


IF{J.E0.1)  GO  TO  999 
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c 
c 


c 
c 

c 
c 


AOl  CONTIMUF 

D  =  X(3  )':^X(3)>!=X(5)/C(8)/X{2)/X(4)/GN 

Di  =  i./(  io+n=<-D) 

AA=(GN+TAU)-0*OI/X( 1) 

B2  =  TAU^!=DT 

B=(L0SS+B1-B2) /X(l ) 
IFCKCC.EQoi)  VIRITE(6,LISTC) 
OUT( 11)=1,-AA 
0UT(]2)=B 

QX=CMPLX(QUT{11)  tOUT(12)  ) 
QS1=CSQRT{QX) 
0UT(13)=REAL(0S1)*X(1)/CN 
0UT(14)=AIMAG(QS1)*X(1 ) /CN 


IF(J„NE.3)  GO  TO  999 


C 
C 


c 
c 
c 


QVI=-(QD*0I-0G^^QF)/QY1 
OUT{'5)=CABS(0VI  ) 
T1=RFAL(QVI) 
T2=AIMAG(QVI) 
0UT(6)=ATAN2(T2,T1) 

QVE=-(OF^:<QG-QD-^QH)  /QYl 
GUT(7)=CABS{QVF) 
IF(KCC.EQ,1)  WRITE(6,LISTG) 
Tl  =  Rf:AL(OVE) 
T2=AIMAG(QVE) 
0UT(8)=ATAN2(T2,T1 ) 

Fl=le03284'CL/CN*X(  1 )  **2+l.  24467- SORT  (  X(  1  )*CL*CN)/C(  10) 

F2==1.2^1=^GL/rM*X(  1  )**2  +  2.0  744=^  SORT  (  X(  1  )'-^CL-M:N  ) /C  (IJ  J 

QN1=1,+Q0^VF1/X(  1) 

CN2  =  l«-o4*0Q*F2/X(  1  )/(  1»  +  QQ^^F2/X(  1  )  ) 

QN=QN1/QM2 

QSl=CSQRT(qM) 

OUT{10)=A1MAG(QS1)*X(1 ) /CN 


D0UT=0UT(4)/0UT( 10) 

D0UT2=nGUT-loDC 

0UT{9)=D0UT2 

0UT(  15)=OUT(4) /SQRT(X( 1 )/6«263183) 

0UT{16)=n!jT(  l'^  )/SORT(X{  l)/6,  28  3183) 

OUT(  17)=rjUT(4)-0UT(  IC  ) 

0UT(19)=CN**2--TAU 


999  RETURN 
END 


163 


BIBLIOGRAPHY 


1.  Part low,  James  G.,  Measurement  of  Acoustic  Pressures 

Associated  with  Moving  Strlations  in  a  Neon  Glow 
Discharge,  M.S.  Thesis,  Naval  Postgraduate  School 
Monterey,  I963. 

2.  Carretta,  A. A.  Jr.,  and  Moore,  W.N.,  Acoustical  Effects  in  a 

Neon  Glow  Discharge,  M.S.  Thesis,  Naval  Postgraduate 
School,  Monterey,  I965. 

3.  Wilson,  O.B.  Jr.,  Observation  of  Acoustic  Waves  Generated 

in  a  Neon  Glow  Discharge,  paper  presented  at  meeting  of 
Acoustical  Society  of  America,  Session  6h6 ,  Boston, 
3  June  1966. 

4.  Crandall ,  J.L.,  Jr.,  Acoustical  Perturbation  of  a  Neon  Glow 

Discharge,  Ph.D.  Thesis,  Naval  Postgraduate  School, 
Monterey,  I967. 

5.  Melville,  R.D.S.  Jr.,  Acoustic  Interaction  with  the  Positive 

Column  of  a  Neon  Glow  Discharge,  M.S.  Thesis,  Naval 
Postgraduate  School,  Monterey,  19^7 • 

6.  Ingard,  U.  ,  "Acoustic  Wave  Generation  and  Amplification  in  d, 

Plasma,"  Phys .  Rev. ,  v.  l45,  4l  (I966). 

7.  Ingard,  U.,  and  Schulz,  M.,  "Acoustic  Wave  Mode  in  a  VJeakly 

Ionized  Gas,"  Phys.  Rev.,  v.  158,  IO6  (I967). 

8.  Yatsui,  K.,  Kobayashi,  T.,  and  Inuishi,  Y. ,  "Ultrasonic  Wave 

Amplification  in  Gaseous  Plasma,"  Tech.  Reports  Osaka  U. , 
V.  18,  No.  812  (1968). 

9.  Yatsui,  K.,  Kobayashi,  T . ,  and  Inuishi,  Y. ,  "Ultrasonic  Wave 

Amplification  in  Gaseous  Plasmas,"  J.  Phys.  Soc.  Japan, 
V.  24,  1186  (I968). 

10.  Schulz,  H.M.  Ill,  Interactions  of  Electrons  and  Acoustic  Waves 

in  a  Plasma ,  Ph.D.  Thesis,  Massachusetts  Institute  of 
Technology,  Cambridge,  19^7. 

11.  Watanabe,  M.S.  and  Oleson,  N.L.,  "Traveling  Density  Waves  in 

Positive  Columns,"  Phys .  Rev. ,  v.  99,  1701  (1955). 

12.  Wojaczeck,  Karl,  "Schall  im  Nieder druckplasma" ,  Beitr . 

Plasmaphys . ,  v.  3,  127  (I96I). 


164 


13.  Strickler,  S.D.,  and  Stewart,  A.B.,  "Radial  and  Azimuthal 

Standing  Sound  Waves  in  a  Glow  Discharge,"  Phys .  Rev. 
Ltrs.  ,  V.  11,  527  (I963). 

14.  Riley,  J.R.,  and  Hall,  D.S. ,  "Neutral  Particle  Density 

Modulation  in  the  Positive  Colume  of  a  Mercury  Discharge," 
Proceedings,  VII  International  Conf.  on  Ionization  Phen. 
in  Gases .  Panel  3.1-^. 

15'       Schulz,  M.  and  Ingard,  U.,  "Acoustic  Kink  Instability  in  an 
Argon  Discharge,"  Phys.  Fluids,  v.  10,  IO31  (1967). 

16.  Fitaire,  M.  and  Mantei,  T.,  "Acoustic  Wave  Generation  by 

Temperature  Modulation  of  a  Plasma,"  Phys .  Rev .  Ltrs . , 
v.  29A,  84  (1969). 

17.  Turner,  R.,  Hochberg,  A.K.  and  Poehler ,  T.O.,  "Multiple  Pulse 

Emission  From  a  HCN  Laser,"  Appl.  Phys.  Ltrs.,  v.  12, 
104  (I968)  . 

18.  Born,  G.K.,  "Acoustic  Perturbations  in  Pulsed  Gas  Lasers," 

J.  Appl.  Phys.,  V.  39,  4479  (1968). 

19.  Hayess,  Erhard,  "Die  Dampfung  von  Schallwellen  in  der 

Entladungsalilo  des  Miede'*' druckplasmas  ,  "  Bei'tr.  Plasmaphvs. 
v.  4,  211  (1964). 

20.  Hayess,  Erhard,  "Schallimpulse  im  Nieder druckplasma ,  I. 

Erzeugung  und  Nachweis  der  Schallimpulse,"  Beitr . 
Plasmaphys. ,  v.  6,  377  (I966). 

21.  Nygaard,  K.J.  and  Meltz,  G. ,  "Sound  Production  in  Glow 

Discharges",  Proceedings,  VIII  Int.  Conf.  on  Phenomena  in 
Ionized  Gases ,  Vienna  (I967). 

22.  Nygaard,  K.J.  and  Meltz,  G.,  "Generation  and  Detection  of 

Acoustic  Waves  in  Pulsed  Electrical  Discharges,"  J.  Acoust 
Soc .  Am . ,  V.  44,  1566  (1968). 

23.  Goldstein,  L.  Roux,  M.R.  and  Dayton,  J. A.  Jr.,  "Production  of 

Acoustic  Waves  by  R.F.  Breakdown  in  Low  Pressure  Gases," 
Proceedings,  VI  Inter.  Conf.  on  Ionization  Phenomena  In 
Gases,  v.  Ill,  115  (I963). 

24.  Berlande,  J.,  Goldan,  P.D.  and  Goldstein,  L ., "Format ion  and 

Propagation  of  a  Pressure  Wave  in  a  Weakly  Ionized  Gas," 
Appl.  Phys  .  Ltrs. ,  v.  5,  51  (1964). 


165 


25.  Born,  G.K.  and  Buser ,  R .G. , "Determination  of  the  Gas 

Temperature  in  the  Afterglow  of  Pulsed  Discharges  by 
Microwave  Probing  of  Standing  Acoustic  Waves,"  J .  App 1 . 
Phys . ,  V.  37,  4918  (1966) . 

26.  Macon,  J.L.,  Generation  and  Amplification  of  Longitudinal 

Waves  in  a  Plasma,  M.S.  Thesis,  Massachusetts  Institute 
of  Technology,  Cambridge,  (I966). 

27.  Tanenbaum,  B.S.,  "Wave  Propagation  in  a  Partly  Ionized  Gas," 

Phys.  Fluids,  v.  5,  1226  (1962). 

28.  Sessler,  G.M.,  "Propagation  of  Longitudinal  Waves  in  a 

Weakly  Ionized  Gas,"  Phys.  Fluids,  v.  7,  90  (1964). 

29.  Veselovskiy,  I.I.,  "Sound  Propagation  in  Partly  Ionized 

Plasma,"  Geomagn,  and  Aeronomy  7,  1,  35-42  (I967). 

30.  Ishii,  T.,  and  Inuishi,  Y.,  "Amplification  of  Neutral  Waves 

in  Gaseous  Plasmas,"  Jap.  J.  Appl.  Phys.,  v.  8,  1531 
(I969). 

31.  Aubrecht,  L.,  "Theoretical  Study  of  the  Interaction  of  Acoustic 

Waves  with  a  Plasma,"  Phys .  Ltrs . ,v.  27A,  526  (1968). 

32.  Aubrecht,  L.,  "General  Solution  of  a  System  of  Equations  for 

Acoustic  Wave  Propagation  in  a  Partially  Ionized  Gas  in  an 
External  Electric  Field",  Czech.  J.  Phys. ,  B19,  1322 
(I969). 

33.  Aubrecht,  L.  ,  "Low  Frequency  Waves  in  a  Prjsrtive  Column  of 

Discharge,"   Czech.    J.    Phys.,    B20 ,    12    (1970). 

34.  Browne,    E.T.,     Introduction    to   the   Theory    of    Determinants    and 

Matr  ices ,    Univ.    of  N.C.    Press.,    1958. 

35.  Morse,    P.    and   Ingard,    U. ,    Theoretical    Acoustics,    McGraw  Hill, 

1968. 

36.  Present,  R.D.,  Kinetic  Theory  of  Gases,  McGraw  Hill,  1958. 

39.  McDaniel,  E.W.,  Collision  Phenomena  in  Ionized  Gases, 

McGraw  Hill,  1964. 

40.  Von  Engel,  A.,  Ionized  Gases  ,  2  ed.  ,  Oxford  University  Press, 

1968. 

41.  Tanenbaum,  B.S.,  Plasma  Physics,  348,  McGraw  Hill,  I967. 


166 


42.  Wu,  Ta-you,  Kinetic  Equations  of  Gases  and  Plasmas,  Addison- 

Wesley  Pub.  Co.,  I966. 

43.  Cooper,  A.W.M.,  Moving  Striations  in  Inert  Gas  Glow  Discharges, 

Ph.D.  Thesis,  Queens  University,  Belfast,  I96I. 

44.  Brown,  S. ,  Basic  Data  of  Plasma  Physics,  (MIT  -  Wiley  1959). 

45.  Taylor,  B.N.,  Parker,  W.H.,  and  Langenberg,  D.N. , "Determination 

of  e/h ,  using  Macroscopic  Quantum  Phase  Coherence  in 
Superconductors,"  Rev.  Mod.  Phys.,  v.  4l,  375  (1969). 

46.  Cook,  G.A.,  ed. ,  Argon,  Helium  and  the  Rare  Gases,  v.  1., 

Interscience  Publishers,  I96I. 

47.  Kinsler,  L.E.,  and  Fr ey ,  A.R.,  Fundamentals  of  Acoustics, 

2  ed. ,  226,  John  Wiley  and  Sons,  I962. 

48.  Bracewell,  R.N.,  The  Fourier  Transform  and  Its  Applications, 

McGraw  Hill,  I965. 

49.  Hayess,  E. ,  Venzke,  P.  and  Wojoczek,  K.,"Die  Neutragas- 

temperatur  in  der  Argon-Mitteldruckent ladung , "  Beitr 

J-clS  tUapuy -a  .,   v./,   -fwj.  y  1.  ^  O  /    J  . 

50.  Hayess,  Erhard,  "Schallimpuls e  im  Nieder dr uckplasma  II. 

Die  Schallausbr eitung , "  Beitr  Plasmaphys.,  v.  7,    327 
(1967). 

51.  Venzke,  D.,  Hayess,  E.  and  Wojoczek,  K,  "Similarity  Relations 

for  Discharge  Columns  in  Noble  Gases  at  Medium  Pressures," 
Beitr.  Plasmaphys,  v.  6,  (I966). 

52.  Chalupnik,  J.D. ,  Rule,  E.,  and  Suellentrop,  F. J. , "Pr essure 

Response  of  Condenser  Microphones  at  Low  Ambient  Pressures," 
J.  Acoust .  Soc.  Am. ,  v.  33,  p  177  (I96I). 

53-  Dayton,  J. A.,  Verdeyen,  J.T.,  and  Virob5k,  P. F., "Method  for 

Detecting  Weak  Sound  Waves  in  a  Low  Pressure  Gas,"  Rev . 
Sci.  Ins.,  V.  34,  1451  (I963). 

54-  Geppert,  D.V.,  Basic  Electron  Tubes,  p.  81,  McGraw  Hill,  1951. 

55.  Cobine,  J.D. ,  Gaseous  Conductors,  p.  128,  McGraw  Hi]l,  194l. 

56.  Rosebury,  F.,  Handbook  of  Electron  Tube  and  Vacuum  Techniques, 

Addison  VJesley,  I965. 


167 


Unclassified 


Setuntv  Classification 


DOCUMENT  CONTROL  DATA  -R&D 


Security  clas  ii  I  ic  alion  of  title,    body  of  absfrarf  and  indexing  annotation  must  be   entered  when   //le  overall  report  is   classified) 


1     ORIGINATING    ACTIVITY   (Corporate  author) 

Naval  Postgraduate  School 
Monterey,  California  93940 


2«.   REPORT    SECURITY     CLASSIFICATION 

Unclass  if ied 


2b.     GROUP 


3      REPOR  T     Tl  T  L  E 

Creation   and  Modification    of  Acoustic    Signals    by    Discharge  Plasma 


4     DESCRIPTIVE   NOT  ES  (Type  of  report  end,inctijsive  dales) 

Ph.D.  Thesis.  June  1971 


5     Au  TmORiSi  (First  name,  middle   initial,    last  name) 

Richard   Willis    Tripp,    Jr. 


6      REPOR  T    D A  TE 

June    1971 


7».     TOTAL    NO.    OF    PAGES 

170 


76.    NO.    OF    RE  FS 

56 


ea.    CONTRACT    OR    GRANT    NO. 


b.    PROJEC  T    NO 


9a.    ORIGINATOR'S    REPORT    NUMBER(S) 


9b.    OTHER    REPORT   NO(S)  (Any  other  numbers   that  may   be  at  signed 
this  report) 


10      DISTRIBUTION    STATEMENT 


Approved  for  public  release;  distribution  unlimited, 


\1      SUPPLEMENTARY    NOTES 


12.    SPONSORING    MILITARY    ACTIVITY 

Naval   Postgraduate    School 
Monterey,    California       939^0 


13      ABSTRACT 


Sound  production  and  modification  in  an  Argon  discharge  plasma  was  investi- 
gated theoretically  and  experimentally.   Sound  production  by  ions  in  the 
cathode  region  was  shown  to  be  feasible  and  consistent  with  ear  lier  ,  e>rpGr  i- 
mental  results.   A  numerical  analysis  of  the  theoretical  dispersion  relation 
of  Ingard  and  Schulz  showed  that  plasma  amplification  effects  will  not 
compensate  for  acoustic  losses.   The  theory  of  the  mobility-limited  ther- 
mionic diode  operated  as  a'  microphone  was  expanded.   Such  microphones  were 
used  in  an  experimental  investigation  of  sound  production  and  modification 
by  a  discharge  plasma.   Possible  plasma  effects  on  sound  propagation  were 
observed.   The  characteristics  of  sound  production  by  the  plasma  were  found 
to  be  in  reasonable  agreement  with  theory,   A  departure  from  theory  at  low 
frequencies  was  interpreted  as  arising  from  the  effects  of  particle  dif- 
fusion, an  effect  not  previously  observed  in  acoustics. 


FORM 
1  NOV  e 

S/N    01 01 -807-68) 1 


,1473     ^'^^^^  '> 


169 


Unclassified 


Security  Classification 


i-31<CS 


Unclassified 


Secuntv  Classifiration 


KEY     WO  R  OS 


Acoustic  interaction 

Plasma  acoustics 

Acoustic 

Plasma 

Diode  microphone 

Microphone 


)D/:o?„1473     fBACK) 


/'.        101-%C7-69Z1 


170  


Unclassified 


Security  Classification 


•  a  2  7  G 


Thesfs 

T799 
c.l 


4  NO/  V  J 


)  0  2  7  (J 


Thesis  128(^52 

T799    Tripp 

c.l        Creation  and  modifi- 
cation of  acoustic 
signals  by  discharge 
plasma. 


thesT799 


^L!,f."°"  ^"d  modification  of  acoustic  si     ' 


3  2768  002  03645  1 

DUDLEY  KNOX  LIBRARY 


